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CHAPTER  I 


INTRODUCTION 

Enclosed  waveguiding  structures  have  long  been  of  great 
importance  in  the  transmission  of  high-frequency  electro¬ 
magnetic  energy.  The  electromagnetic  properties  of  uniform, 
axial  waveguides  are  generally  well  understood,  but  practical 
applications  introduce  discontinuities  and  transition  regions 
into  waveguides  which  often  significantly  affect  their  re¬ 
sponse.  In  the  particular  case  of  discontinuities  which  are 
limited  to  the  plane  transverse  to  the  direction  of  propa¬ 
gation,  an  impedance  mismatch  and  a  region  of  reactive  energy 
storage  are  normally  introduced  into  the  waveguide.  Such 
effects  are  usually  undesirable,  but  properly  designed  dis¬ 
continuities  can  provide  very  desirable  impedance-matching 
devices  and  resonating  cavities. 

Analyses  of  planar  waveguide  discontinuities  quickly  fol¬ 
lowed  the  beginning  of  efforts  to  develop  useful  radio-fre¬ 
quency  and  microwave  devices  utilizing  waveguiding  structures. 
Low-frequency  techniques  such  as  field  mapping  yielded  useful 
results  for  transverse  electromagnetic  (TEM)  mode  propagation 
in  coaxial  and  parallel-plate  waveguides.  The  effective 
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treatment  of  discontinuities  when  higher-order  modes  exist 
in  coaxial  lines  and  for  waveguides  supporting  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  modes  was  much  more 
difficult.  A  discontinuity  requires  the  local  contributions 
of  an  infinite  set  of  modes  to  satisfy  the  electromagnetic 
boundary  conditions.  Mathematical  complexity  prevented  early 
analytical  treatments  from  considering  more  than  one  or  two 
of  the  higher-order  modes.  The  most  accurate  of  the  early 
treatments  was  based  upon  variational  solutions  to  susceptance 
expressions  which  had  been  developed  from  integral  equations. 
Very  useful  expressions  for  the  susceptances  of  a  large  variety 
of  waveguide  discontinuities  were  derived  by  Julian  Schwinger 
and  a  capable  group  of  associates  in  the  1940s.  Many  of  the 
results  of  this  work  were  compiled  by  Nathan  Marcuvitz  into 
the  well-known  Waveguide  Handbook,  Ref.[l],  This  work  is  of 
such  quality  that  even  today  it  remains  the  treatment  of 
choice  for  most  applications.  Very  good  descriptions  of  the 
techniques  employed  to  derive  these  expressions  have  been 
presented  in  books  by  Schwinger  and  Saxon,  Ref. [2],  Collin, 

Ref. [3],  and  Lewin,  Refs. [4],  [5],  and  [6].  The  earliest 
published  work  on  discontinuities  in  coaxial  waveguides  was 
by  J.  R.  Whinnery  and  H.  W.  Jamieson,  Ref. [7],  and  Whinnery, 
Jamieson,  and  Robbins,  Ref. [8].  These  investigators  used  an 
approach  suggested  by  W.  C.  Hahn  in  1941  which  involved  com- 
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putation  of  the  terms  of  a  Fourier  expansion  of  the  field  in 
the  discontinuity  to  calculate  the  capacitance  of  coaxial  dis 
continuities.  Computation  of  the  terms  of  the  series  was 
sufficiently  difficult  that  Whinnery,  et  al.,  used  only  the 
first  four  higher-order  modes  to  obtain  results  that  were 
later  shown  to  be  accurate  to  within  five  per  cent. 

The  variational  technique  employed  by  the  above  authors 
is  included  in  the  collection  of  mathematical  solution  pro¬ 
cedures  which  R.  F.  Harrington  has  identified  in  his  Field 
Computation  by  Moment  Methods ,  Ref. [9].  This  variational 
approach  is  identical  to  a  one-term  Galerkin's  solution. 

Ref. [9],  where  the  assumed  value  of  the  unknown  field  dis¬ 
tribution  in  the  discontinuity  serves  as  both  the  basis  and 
testing  function.  As  indicated  earlier,  this  approach  is 
very  effective  for  dominant  mode  operation  of  waveguides  when 
frequencies  are  low  enough  to  prevent  propagation  of  higher- 
order  modes.  Under  this  frequency  constraint,  susceptance  re 
suits  obtained  from  both  analyses  and  experimental  measure¬ 
ments  during  the  course  of  this  investigation  have  shown  re¬ 
markable  agreement  with  values  calculated  from  expressions  in 
the  Waveguide  Handbook.  These  latter  calculations  benefit 
from  the  characteristic  of  such  variational  solutions  in  that 
a  first-order  change  in  the  assumed  field  distribution  in  the 
discontinuity  produces  only  a  second-order  variation  in  the 
calculated  susceptance.  The  field  in  the  vicinity  of  the  dis 
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continuity  can  be  calculated  from  the  same  assumed  field  dis¬ 
tribution,  but  the  field  expressions  are  not  variational  and 
thus ,  do  not  share  the  property  whereby  a  first-order  error  in 
the  assumed  field  distribution  produces  a  second-order  error 
in  the  calculated  result. 

The  treatments  of  waveguide  discontinuity  problems  which 
employed  variational  techniques  required  a  liberal  amount  of 
ingenuity  because  such  problems  yielded  solutions  very 
grudgingly.  The  increasing  capability  of  the  digital  computer 
opened  new  avenues  of  approach  and  investigators  began  to  de¬ 
velop  techniques  which  took  advantage  of  this  tool.  P.  I. 
Somlo  used  the  computer  to  calculate  more  terms  in  the  Hahn 
series  approach  used  by  Whinnery,  et- al. ,  Refs. [10]  and  [11]. 
D.  Woods,  Ref. [12],  used  interpolation  on  Somlo's  results  to 
obtain  a  further  improvement  in  accuracy  and  A.  Jurkus  pro¬ 
vided  improvements  to  Somlo's  code.  Ref. [13].  E.  W.  Risley, 
Refs. [14]  and  [15],  adapted  the  Rayleigh-Ritz  procedure  used 
by  Schwinger,  et  al . ,  to  solution  by  a  computer  and  obtained 
both  upper-  and  lower-bound  solutions  for  the  capacitance  of 
the  coaxial- to-circular  waveguide  junction. 

A  number  of  more  recent  investigations  has  employed 
least-squares  techniques  to  minimize  the  error  in  satisfying 
boundary-conditions  at  the  discontinuities.  The  error  results 
from  the  fact  that  only  a  finite  number  of  higher-order  modes 
is  used  in  the  field  expansions  at  the  discontinuities. 
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Investigations  using  such  methods  have  been  reported  by  J.  B. 
Davies,  Ref. [16],  Davies  and  M.  Razaz,  Ref. [17],  H.  Oraizi 
and  J.  Perini,  Ref. [18],  R.  Jansen,  Ref. [19],  and  H.  J.  A. 
LaRiviere  and  Davies,  Ref. [20].  The  results  reported  in  Refs. 
[16]  -  [20]  are  largely  for  general  waveguide  discontinuities 
with  occasional  specialization  of  the  technique  to  coaxial 
waveguides  as  in  Ref. [17]. 

All  of  the  investigations  mentioned  thus  far  were  direc¬ 
ted  toward  the  computation  of  discontinuity  capacitance  and 
some  have  obtained  very  good  levels  of  accuracy.  Razaz  and 
Davies  claim  an  error  of  ±0.1  femtofarads  in  their  calcula¬ 
tions  for  the  capacitance  of  a  coaxial-to-circular  waveguide 
junction.  Ref. [21].  None  of  these  approaches  is  suitable 
for  treating  situations  where  higher-order  modes  may  propa¬ 
gate.  Alvin  Wexler,  Ref. [22],  has  developed  a  modal  analysis 
technique  for  treating  waveguide  discontinuities  which  in¬ 
volves  selecting  the  amplitudes  of  normal  modes  so  as  to 
satisfy  boundary  conditions  at  the  discontinuity.  The  method 
is  elegantly  developed  and  allows  the  treatment  of  as  many 
modes  as  desired.  It  has  the  computational  disadvantage  of 
requiring  the  solution  of  a  set  of  N+l  linear  equations  to 
treat  N  modes.  The  method  can  also  treat  coupled  discontinu¬ 
ities,  however,  each  discontinuity  must  be  treated  individually 
and  the  procedure  must  begin  with  the  discontinuity  most 
distant  from  the  excitation. 
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All  the  solution  techniques  mentioned  thus  far  begin 
with  the  same  expressions  for  the  field  in  the  waveguide.  It 
is  the  approach  to  approximating  or  obtaining  solutions  for 
the  transverse  field  in  the  apertures  of  the  discontinuities 
that  provides  the  differences  among  the  methods.  The  primary 
objective  of  this  investigation  is  the  development  and  vali¬ 
dation  of  a  frequency-independent  capability  of  calculating 
the  electromagnetic  properties  of  coaxial  waveguide  disconti¬ 
nuities  which  is  both  computationally  efficient  and  accurate. 
Such  a  method  necessarily  includes  the  treatment  of  higher- 
order  modes,  but  it  is  done  in  such  a  way  that  it  allows  for 
general  excitation  of  the  waveguides.  The  method  presented 
possesses  a  significant  advantage  over  earlier  approaches  in 
that  a  series-convergence  technique  is  employed  which  allows 
the  recovery  of  essentially  all  the  higher-order  mode  con¬ 
tributions  at  a  given  discontinuity  even  though  a  finite 
number  of  series  terms  is  computed.  The  accuracy  of  the 
method  can  be  controlled  (within  computational  errors)  by 
the  number  of  basis  expansion  pulses  and  higher-order  mode 
terms  considered  in  the  solution.  A  very  useful  aspect  of 
this  approach  is  the  capability  of  combining  two  or  more 
discontinuities  into  simple  or  complex,  one-  or  two-port 
cavities  as  illustrated  in  Figure  1. 

The  formation  and  analytical  treatment  of  simple  and 
partitioned  cavities  are  considered  extensively  in  the  body 
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of  this  report.  Also  presented  are  the  results  of  a  sub¬ 
stantial  experimental  program  designed  to  provide  validation 
for  the  analysis.  Two  sets  of  coaxial  cavities  with  differing 
dimensions  were  constructed  and  the  internal  field  was  meas¬ 
ured  for  one-  and  two-port  cavities  of  simple  and  parti¬ 
tioned  configurations.  The  input  admittance  of  the  two-port 
cavity  was  measured.  Extensive  comparisons  of  measured 
data  with  analytical  results  are  presented  which,  in  con¬ 
junction  with  a  discussion  of  possible  experimental  error, 
lead  to  quantitative  evaluations  of  the  accuracy  and  useful¬ 
ness  of  the  analytical  approach. 

Analytic  formulations  for  the  field  or  electromagnetic 
potentials  within  a  cavity  frequently  include  expressions 
that  contain  three  infinite  series.  For  coaxial  cavities 
considered  in  this  effort,  the  original  excitation  is  always 
assumed  to  consist  of  TEM  or  TMQn  modes  only.  This  constraint 
eliminates  any  variation  in  the  azimuthal  coordinate  and 
simplifies  the  appropriate  expressions  to  two  infinite  series. 
It  is  mathematically  (and  computationally)  expedient  to  treat 
these  cavities  from  a  waveguide,  i.e.,  traveling  wave,  point 
of  view  and  to  consider  the  front  and  rear  boundaries  of  the 
cavities  as  planar  waveguide  discontinuities.  The  axial  and 
radial  eigenvalues  in  this  approach  are  related  in  such  a 
manner  that  the  formulation  requires  only  one  infinite  series 
for  a  complete  set  of  eigenfunctions.  The  derivation  of  the 


magnetic  vector  potential  and  the  resulting  field  expressions 
is  presented  in  Chapter  II. 

If  the  waveguide  approach  is  used,  it  is  very  effective 
to  think  of  the  traveling  waves  as  penetrating  apertures  in 
the  discontinuities.  This  concept  allows  utilization  of 
existing  techniques  for  calculating  aperture  penetration,  in 
particular,  treating  the  aperture  field  as  the  unknown 
quantity  of  interest  since  it  is  common  to  the  adjoining  re¬ 
gions.  All  formulations  in  this  report  treat  the  transverse 
electric  field  in  the  aperture  as  the  unknown  quantity.  It 
is  the  objective,  therefore,  to  obtain  solutions  for  the 
aperture  field  for  each  discontinuity  or  partition  used  to 
form  a  cavity  structure.  When  this  approach  is  used,  it  is 
imperative  to  recognize  any  singular  behavior  of  the  aperture 
field  and  then  to  develop  a  numerical  solution  technique 
accordingly.  In  the  following  chapters  integral  equations 
are  developed  for  progressively  more  complex  configurations 
of  discontinuities  and  details  are  given  for  appropriate 
numerical  solution  techniques.  The  results  of  the  comple¬ 
mentary  experimental  measurements  program  serve  to  corrob¬ 
orate  the  analytical  results. 
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Figure  1.  Examples  of  Cavities  Formed  by  Placing  Planar  Discontinuities  in 
Coaxial  Waveguides 


CHAPTER  II 

ROT AT I ONALL Y  SYMMETRIC  JUNCTION  OF 
TWO  DISSIMILAR  COAXIAL  WAVEGUIDES 

A.  Introduction 

The  effects  on  the  propagating  modes  produced  by  a  dis¬ 
continuity  in  a  coaxial  waveguide  are  a  function  of  the  re¬ 
sulting  field  distribution  in  the  vicinity  of  the  discontin¬ 
uity.  The  local  field  distribution  is  a  modification  of  the 
propagating  modes  and,  in  general,  an  infinite  number  of 
higher-order  modes  is  required  to  satisfy  boundary  conditions 
at  such  a  discontinuity.  For  a  discontinuity  confined  to  a 
two-dimensional  plane,  a  knowledge  of  the  field  in  this 
"aperture  plane"  is  sufficient  for  determining  the  field 
everywhere.  From  the  knowledge  of  the  aperture  field,  the 
resulting  effects  on  any  propagating  modes  can  be  determined 
and,  if  desired,  microwave  circuit  representations  of  the 
discontinuity  can  be  derived.  In  this  chapter,  an  integral 
equation  is  formulated  for  the  aperture  electric  field  of  the 
generalized  geometry  depicted  in  Figure  2.  This  integral 
equation  is  then  specialized  to  the  simpler  geometries  of 
Figures  3  and  4  and  is  subsequently  applied  to  a  one-port 
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coaxial  cavity. 


B.  Derivation  of  Expressions  for  Field  Components 

An  integral  equation  suitable  for  numerical  solution  is 
obtained  by  developing  expressions  for  each  of  the  field  com¬ 
ponents  in  the  regions  either  side  of  the  discontinuity  and 
then  enforcing  continuity  of  the  transverse  field  components 
across  the  aperture  in  the  discontinuity.  These  expressions 
are  derived  from  the  appropriate  operations  on  a  vector 
potential  of  the  magnetic  and/or  electric  type.  The  geometry 
under  consideration  is  azimuthally  symmetric  and,  with  the 
assumption  of  similar  symmetry  of  the  excitation,  all  magnetic 
field  components  are  transverse  to  the  axial  coordinate  and, 
therefore,  the  magnetic  vector  potential  suffices  for  the 
complete  description  of  the  field.  This  potential  is  the 
solution  to  the  vector  wave  equation 

(V2  +  k2)  X  =  ft,  (1) 

subject  to  the  appropriate  boundary  conditions  and  where  k  is 
the  wave  number  of  the  medium.  A  coaxial  geometry  with 
azimuthally- symmetric  excitation  has  only  axially-directed  cur¬ 
rents  and,  therefore,  only  the  axial  component  of  the  magnetic 
vector  potential  is  necessary  to  derive  the  complete  field. 

The  geometries  treated  in  this  investigation  have  many  surfaces 
that  support  p-directed  currents;  however,  it  is  shown  in 


Appendix  A  that,  for  azimuthally-symmetric  excitation,  the 
field  derived  from  the  z-directed  vector  potential  is  still 
complete.  The  vector  potential  equation  is  reduced  to  the 
scalar  equation 

(V2  +  k2)  A  (p  ,  <j> ,  z )  =  0. 

2 


Since  the  axial  component  of  vector  potential  is  independent 

of  <p ,  A  may  be  represented  by 
z 

Az(p,z)  =  R(p)Z(z) . 


The  axial  function  Z{z)  satisfies  the  harmonic  equation 


(2) 


which  has  the  solution 

Z (z)  *  A  e"jkzz  +  B  e 

where  A  and  B  are  arbitrary  constants 
R (p )  is  a  solution  of 


jkzz 

.  The  radial 


function 


2  d  R 


dp 


p  ^  +  k" 
dp  p 


(3) 


The  parameters  kz  and  kp  are  related  by 


(4) 
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Equation  (3)  is  a  zero-order  Bessel's  equation  with  the 
solution 


R(P)  =  C  JQ(kpp)  +  D  NQ(kpp) 

where  Jq  is  the  zero-order  Bessel  function  of  the  first  kind, 
Nq  is  the  zero-order  Neumann  function  and  C  and  D  are  arbi¬ 
trary  constants.  Enforcement  of  the  boundary  condition  that 
the  tangential  (or  z)  component  of  electric  field  be  zero  at 
the  surface  of  the  inner  conductor  yields  the  function 

R(p)  =  E  [N0(kpa)J0(kpp)  -  J0(kpa)N0(kpP)] 

where  a  is  the  radius  of  the  inner  conductor  and  E  is  an  ar¬ 
bitrary  constant.  The  second  boundary  condition  requires  the 
z  component  of  the  electric  field  to  be  zero  at  the  surface 
of  the  outer  conductor  and  results  in  the  expression 

R(P)  =  *n(p)  =  W^VV5*  "  J0(Yna)N0(Ynp)  (5) 

where  the  are  the  infinite  set  of  k  which  satisfy 

Vb)  =  N0(Yna)J0(Ynb)  “  J0(Yna,N0{Ynb)  *  °'  n=1'2'*-*  <6) 

where  b  is  the  radius  of  the  outer  conductor.  The  new  rep¬ 
resentation  for  k  is  a  where 

z  n 
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For  the  special  case  kp  «  0,  the  solution  to  Eq.  (3)  is 

R{p)  =  B'  lnp  +  C’ 

where  the  constant  C'  can  be  discarded.  In  this  case, 
k  =  k  which  yields  the  special  case 

A^EM(p,z)  =  b+  lnp  e"jkz  +  b“  lnp  ejkz  (8) 

mpu 

where,  as  indicated,  is  the  vector  potential  for  the 

TEM  mode.  The  complete  expression  for  the  vector  potential 
is  given  by 

Az(p,z)  =  bj  lnp  e"jkz  +  bg  lnp  ejkz 

+  l  (b*  e'jotnz  +  b~  ejotnz)  *  (p)  (9) 

I*  il 

n— l 

The  wave  number  an  for  the  higher-order  modes  may  be 
either  real  or  imaginary  depending  on  the  relative  magni¬ 
tudes  of  k  and  y  .  For  the  former  condition,  the  mode  would 
propagate  and,  for  the  latter,  it  would  be  cutoff.  Proper 
root  selection  is  required  for  each  case  to  assure  conser¬ 
vation  of  energy.  The  following  selections  are  made 

an  *  42  -  *n  '  *2  »  *n  <10a> 

an  *  -3  ^  -  k2  >  ^  (Wl 

In  Eqs.(10a&b),  k  is  assumed  to  be  a  real  number. 
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The  electric  field  is  calculated  from  the  expression 


E  =  -j  %  [k2  2Az  +  77 *  (2AZ)]  . 
k 


The  radial  component  of  the  electric  field  is  given  by 

32A 

Ep  3  jT  3p3z  ' 
p  .  U)  e^f.  ,  u,  h-  e3*2 

Ep  k  b0  p  k  b0  p 

,  oj  \  „  ,  y+  -jot  z  ,  -  jot  z.  d$  ( p ) 

+  — ^  ;  a_  (-b_  e  J  n  +  b_  e  n  )  n 

.  2  L,  n  n  n  — 3— — 

k  n=l  dp 


The  axial  component  of  the  electric  field  is  given  by 

.2  1 


e2  -  -j  % 
2 


k2  +  r 


3z 


Az  ' 


Ez  *  "3  72  l  ^  -  «n>  «£  +  b’n  ^  V» 

k  n=l 


The  magnetic  field  is  calculated  from 


H  -  i  »*(*»,> 


where  y  is  the  permeability  of  the  dielectric  medium, 
azimuthal  component  of  magnetic  field  is  given  by 


H , 


1  ^Z 

V  9P  ' 


(ID 


(12) 


The 
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„  -  I  „♦  -  1  b-  £1 

(j)  y  0  p  y  0  p 


i  I  <b+  e‘*anz  +  fa’  e^anz)  dV0) 

u  n=l  n  n  — dp 


(13) 


At  this  point,  the  constants  used  in  defining  the  field 
components  are  redefined  to  make  the  TEM  terms  agree  with 
physical  concepts  so  that  the  sign  conventions  on  the  TEM 
and  the  higher-order  mode  terms  agree,  e.g.. 


II 

+  O 

X) 

-  -  B+ 
u  B0  ' 

(14a) 

II 

1  o 

XI 

£  B0  ' 

(14b) 

*>;  = 

wa  n 

n 

(14c) 

bn  - 

sLb-  . 

u>a  n 

(14d) 

The  field  components  are  now  given  by 

E  -  e-3k*  +  b;  *  l  (B+  e'3“nz  +  ejonz)  , 

K  P  P  n=l  dp 


Ez  =  j  l 

2  n=l 

and 

H 


f.2  2 

k  -  a 

n 


<B£  e'30tnZ  "  Bn  e3“nZ)  *n  ' 


-  B+  e"3kz  -  B*  ejkz  +  -  7  — (B+  e"janZ  -  B"ejcinZl  d$n 

*  B0  B0  +  r,  A  a  (Bn  e  n  Bne  n  '  ^  ' 

HP  np  n=l  n  dp 


1  7 


(15) 


where  n  is  the  characteristic  impedance  of  the  dielectric 
medium, 


C.  Field  Expressions  for  Two  Dissimilar  Coaxial  Regions 

The  analysis  of  the  junction  of  two  coaxial  regions  with 
differing  inner  and/or  outer  radii  and/or  dielectric  media 
requires  separate  field  expressions  for  each  region.  Such  a 
junction  is  shown  in  Figure  2  which  also  displays  schematic 
representations  of  the  incident  and  reflected  waves  in  each 
region.  In  this  figure  and  others  in  the  report,  the  assump¬ 
tion  is  made  that  the  truncated  sections  of  waveguide  have 
matched  terminations  (or  equivalently,  extend  to  ±  infinity). 
Under  such  an  assumption,  all  incident  waves  represent  know 
sources  with  no  reflected  components.  The  junction  in  Fig¬ 
ure  2  is  arbitrarily  located  at  z  =  s  in  order  to  retain  gen¬ 
erality  in  the  field  expressions.  For  single  discontinuities, 
specifying  z  *  0  in  the  resulting  equations  will  simplify  the 
mathematics  somewhat. 

For  the  region  z  <  s,  with  media  characteristics 


and 


k 


(16a) 


(16b) 


the  field  components  are 


E  *  B„  - 

pup 


+  e-ik-z  ,  eJk-z 


+  B-  S-^  +  I  (B*  .'iV  ♦  .3°»z)«n  ,  ,17) 


n=l 


dp 
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n=l  n 


(18) 


and 

,  -jk  z  ik  z  k  «  , 
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where 


r  l  b<  -  b;  -3“"z>^a 

-  n=l  n  dp 


dp 

(19) 


Vp)  =  N0(Yna)J0(Ynp)  ‘  J0(Yna)N0(Ynp) 


(20) 


and 


a2  =  k2  -  Y2 


n 


(21) 


The  Yn  are  the  roots  of  the  transcendental  equation 

Vb)  -  N0(Yna)J0(Ynb)  -  J0(Yna)N0(Ynb)  =  0  '  {22) 

where  a  and  b  are  the  inner  and  outer  radii  of  the  coaxial 
region  for  z  <  s. 

For  the  region  z  >  s,  with  the  dielectric  media  char¬ 
acteristics 


k+  =  u>  /u+e+ 


(23a) 


and 


= 


(23b) 


the  field  components  are 
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where 


dp 

(26) 


Vp)  =  N0<Xnc)J0(Xnp)  -  J0(Xnc)N0(Xnp)  (27) 


and 


2  2  2 

6  =  k“  -  A  . 

n  +  n 


(28) 


The  An  are  the  roots  of  the  transcendental  equation 


4'n(d)  =  N0(Xnc)  J0(Xnd)  "  J0  (Xnc)  N0  (Xnd)  “  0  '  (29) 


where  c  and  d  are  the  inner  and  outer  radii  of  the  coaxial 
region  for  z  >  s.  Note  that  the  annular  disk  attached  to 
the  inner  conductor  has  no  effect  on  the  field  definitions. 
The  disk  merely  imposes  a  boundary  condition  to  be  enforced 
in  the  plane  of  the  aperture. 

D.  Evaluation  of  Constants  in  the  z  <  s  Region 

The  unknown  constants  of  Eqs . (17 ) - (19)  are  evaluated  in 
terms  of  integrals  over  the  unknown  aperture  electric  field 
by  application  of  Fourier  series  techniques.  The  first  step 
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1 


in  evaluating  Bq  and  Bq  is  to  integrate  both  sides  of  Eq. (17) 
over  the  z  =  s  plane  (illustrated  in  Figure  2)  with  respect 
to  p  which  yields 


b 

Ep 


B, 


jk  s  dp 


b"  Jk-S 
Bq  e 


rb 
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dp 
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b 

a  n=l 


l  (B*  e  -*an“  +  B^  eJ  "n~')  ’“"n  dp 


-  ja  s.d<l> 

l  r*  1  r 

n 


dp 


After  evaluation  of  the  integrals  on  the  right-hand  side  of 
the  equation,  one  obtains 


-  j  “jk  s  .  -  jk  s 

Ep  dp  =  (B0  e  J  -  +  Bq  eJ  - 


)  lnp  | 


+  l  (B+  e~ianS  +  Bn  e^anS)  5n(p) . 


n=l 


Note  that  $R(b)  =  0  from  Eq. (22)  and  that  <t>n(a)  =  0  from 
the  definition  in  Eq. (20) ;  therefore 


B0  e 


-jk-s  -  b:  e^-s  =  A  fbE* 


ln-  ^ 


dp 


(30) 


The  presence  of  the  perfectly  conducting  disk  in  the  aper¬ 
ture  affects  the  value  of  E~  over  the  aperture. 


j”  0  ,  pe (a,h) 
|^E°,  pe  (h,b) 


With  the  constraint  of  Eq.(31),  Eq. (30)  becomes 


(31) 
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B0  e 


"jk-s  +  B«  e^k-s  -  -4  E°  dp  . 

0  ln|  Jh  p 

d 


The  B 


+  e"jk-z 


term  in  Eq. (17)  represents  a  TEM  wave  traveling 


in  the  +2  direction.  This  wave  is  the  only  such  positive 
traveling  TEM  wave  for  z  <  s  and,  hence,  is  identified  as  the 
known  excitation.  The  subject  term  may  be  expressed  as 

„+  e'jlt-2  V*  a-ik-2 


The  quantity  V  is  the  voltage  (complex  in  general)  of  the 
excitation  applied  to  the  coaxial  waveguide  in  the  z  <  s 
region  and  is  defined  by 

rb 


-  -r 


Ep  dp  . 


The  value  of  V  is  normally  a  specified  parameter  of  the 
problem.  The  following  definitions  may  now  be  obtained: 


V+ 

ln- 

a 

e”3k_s 

i  b 

in-  [• 

h  P 

.  v+  *-:>k-s 


(34a) 


(34b) 


The  higher-order  mode  coefficients  are  determined  in  a  sim- 

»  cl$ 

ilar  manner.  The  first  step  is  to  multiply  Eq. (17)  by  p  _ m 

dp 

and  integrate  with  respect  to  p  over  the  z  =  s  plane. 
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B+  e^V3  +  B"  e^n3  = 


n 


n 


d<J>  . 
p  _ n  dp  . 

dp 


(37) 


As  indicated  earlier,  Bg  and  Bg  are  the  complex  coef¬ 
ficients  of  the  TEM  waves  in  Region  1  which  are  traveling 
in  the  positive  and  negative  z  directions,  respectively. 

Under  the  assumption  of  a  matched  waveguide  in  the  negative 
z  direction,  Bg  will  be  non-zero  only  in  the  case  of  an  ex¬ 
ternally  excited  wave  traveling  in  the  positive  z  direction. 
The  TEM  wave  traveling  in  the  negative  z  direction  repre¬ 
sents  the  reflected  component  of  the  incident  TEM  wave  and/or 
the  TEM  component  of  a  transmitted  wave  incident  on  the  aper¬ 
ture  from  Region  2.  The  same  definitions  apply  for  B*  and 
Bn  which  are  the  complex  coefficients  for  the  higher-order 
modes  in  Region  1 .  In  order  for  the  solutions  to  the  prob¬ 
lems  considered  in  the  present  investigation  to  be  tractable, 
the  Bq  and  the  Bn  must  be  either  zero  or  known  quantities. 

In  theory,  there  is  no  limit  to  the  number  of  incident  waves 
which  can  be  treated  as  long  as  the  values  of  the  coeffi¬ 
cients  are  specified.  Satisfaction  of  the  boundary  condi¬ 
tions  at  a  perfect  discontinuity  would  require  the  excita¬ 
tion  of  an  infinite  number  of  higher-order  modes  and,  thus, 
all  Bn  should  be  non-zero.  The  procedure  for  accurately 
evaluating  an  infinite  series  of  Bessel  functions  which  is 
introduced  in  Appendix  B  allows  the  solution  of  such  problems 
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with  a  relatively  small  number  of  terms  of  the  resulting 
series  (normally  less  than  20) . 


E.  Evaluation  of  Constants  in  the  z  >  s  Region 

The  procedures  for  evaluating  the  constants  in  the 
field  expressions  for  the  z  >  s  region  are  the  same  as  for 
the  z  <  s  region.  For  a  single  discontinuity,  C~  and  C~ 
are  the  complex  coefficients  of  waves  incident  from  the  nega 
tive  z  direction  and,  if  non-zero,  are  assumed  to  be  known 


F.  Enforce  continuity  of  Transverse  Fields  in  Aperture 

The  unknown  constants  in  the  field  component  expressions 
for  the  regions  on  both  sides  of  the  discontinuity  are  now 
defined  in  terms  of  weighted  integrals  over  the  unknown  aper¬ 
ture  electric  field.  Fundamental  boundary  conditions  require 
that  the  transverse  components  of  the  electromagnetic  field 
be  continuous  across  the  aperture.  This  constraint  has  al¬ 
ready  been  applied  to  the  electric  field  component  since 


E~(p,s)  *  E*(p,s)  ■  E°  ,  pe(h,b)  . 


Enforcing  continuity  of  the  transverse  magnetic  field  is  ac¬ 
complished  by  equating  Eqs. (19)  and  (27)  at  the  z  *  s  plane, 
i.e.  , 

h~(p,s)  *  h*(p,s)  ,  pe (h,b)  . 

After  the  evaluation  of  B~,  B~,  C*,  and  C*  through  the  use  of 
Eqs. (32),  (37),  (38),  and  (39),  this  equation  becomes 
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(41) 
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The  terms  involving  unknown  constants  are  placed  on  the 
left-hand  side  of  the  equation  to  obtain 


.  - jk  s  ik ,  s  2k  «  .  ,  „ 

2B+  2-LjL-  +  2C:  sii-  +  —  l  i_  B-  e'3ans  dtn 

0  1.0  0  n+p  n.  a„  n  3J- 


2k ,  » 


+  r  1  . j 8_s  df.  .  f  1  .ill  f  „ o 

— —  )  -3—  C_  eJ  n  _ n  =  - r-  +  - -r  —  EA  dp 

n+  n=l  6n  n  dp  n  lnr  n .  In-  p  ' h  p 


a  '+  cl 


k-  r  1  f°o  d$  .  d* 
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A  proper  integral  equation  (i.e.,  one  equation,  one  unknown) 
may  be  obtained  from  Eq. (41)  by  the  application  of  sufficient 
constraints  to  evaluate  the  unknown  constants.  A  commonly 
encountered  situation  is  TEM  excitation  only  and  from  only 
one  direction.  For  an  incident  TEM  wave  in  Region  1  travel¬ 
ing  in  the  positive  z  direction,  B*,  Cq  and  c”  would  be  zero. 
The  constant  Bq  would  then  be  evaluated  from  Eq. (34a) . 

G.  Application  to  Simpler  Geometries 

The  expression  presented  in  Eq.  (41)  pertains  to  the  ge¬ 
ometry  depicted  in  Figure  2  which  is  a  coaxial  waveguide  junc¬ 
tion  which  includes  both  an  annular  disk  on  one  conductor  and 
steps  in  the  radii  of  both  conductors.  Most  coaxial  waveguide 


discontinuities  are  not  this  complex  and  the  appropriate 

integral  equations  are  special  cases  of  Eq. (41) .  It  should 

be  noted  that  the  presence  of  an  annular  disk  on  one  or  both 

coaxial  conductors  affects  only  the  limits  of  integration  of 

the  three  integrals  in  Eq. (41) ;  consequently,  this  type  of 

geometry  variation  is  easily  treated.  For  a  discontinuity 

consisting  of  a  step  in  the  inner  or  outer  radius  only, 

Eq.  (41)  is  applicable  with  suitable  adjustments  of  radii 

2  2 

values  used  in  defining  $  ,  a„ ,  ,  ¥  ,  B  and  U  .  In  ad- 

n  n  n  n  n  n 

dition,  the  proper  changes  must  be  made  in  the  limits  of 
integration.  For  the  case  of  a  discontinuity  consisting 
solely  of  an  annular  disk  attached  to  the  inner  or  outer 
coaxial  conductor,  as  shown  in  Figure  3,  or  both  conductors, 
as  shown  in  Figure  1,  Eq. (41)  is  simplified  considerably. 
Since  the  coaxial  regions  on  both  sides  of  the  discontinuity 
are  now  the  same,  the  $n  of  Eq. (20)  and  the  4^  of  Eq.  (27) 
are  now  the  same.  An  important  and  common  simplification 
occurs  when  the  dielectric  media  in  the  two  regions  are  the 
same  $  i • 6 « / 


and 


=  e 


e 


y+  =  y_  =  y  . 


(42a) 

(42b) 


The  constraints  of  Eq.(42)  result  in  k+  =  k_,  on  =  Bn,  and 
n+  *  n_. 

An  integral  equation  appropriate  for  the  geometry  of 
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Figure  3a  under  the  assumption  that  the  excitation  consists 
only  of  a  TEM  wave  traveling  in  the  positive  z  direction  is 
given  by 
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where  and  a*  pertain  to  the  negative  and  positive  z 

regions  respectively.  If  the  constraints  of  Eq. (42)  are 
,  2 

applied  to  Eq. (43)  a  factor  of  —  may  be  divided  from  each 
term  in  the  equation  (noting  that  the  two  series  would  now 
be  identical) .  It  is  possible  to  change  Eq. (43)  into  an 
equation  appropriate  for  the  geometry  of  Figure  3b  by  chang¬ 
ing  the  limits  of  integration  for  all  integrals  from  (a,g) 
to  (hfb) . 

An  integral  equation  appropriate  for  the  geometry  of 

Figure  4a  with  the  same  excitation  employed  in  Eq.  (43)  and 

2 

with  the  parameters  V  ,  Bn  and  Un  properly  defined  is  given 
by 
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z  *  s 


Figure  4b,  Step  Discontinuity  in  the  Inner  Conductor 
of  a  Coaxial  Waveguide 
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n  n 


This  equation  may  be  changed  to  one  appropriate  for  Figure 

4b  by  changing  the  limits  of  integration  from  (a,b)  to  (c,b) 

and  by  making  the  appropriate  changes  to  the  definitions  of 

6  .  ¥  and  U2. 
n  n  n 


H.  Integral  Equation  for  a  One-Port  Coaxial  Cavity 

An  effective  way  to  create  a  one-port  coaxial  cavity  is 
to  place  a  perfectly  conducting,  transverse  plate  across  the 
coaxial  waveguide  at  a  point  on  the  positive  z  side  of  any 
of  the  discontinuities  depicted  in  Figures  3  and  4 .  The  cav¬ 
ity  treated  in  this  section  is  created  from  Figure  4a  and  is 
shown  in  Figure  5.  The  expressions  for  the  field  components 
in  the  coaxial  regions  are  those  given  by  Eqs. (17) - (19)  and 
(24)  —  (26)  with  the  accompanying  definitions.  Note  that,  for 
the  subject  geometry,  the  radii  a  and  c  are  equal. 

The  constants  in  Eqs. (17) - (19)  are  defined  by  Eqs. (34) 
and  (37)  with  the  assumption  that  the  only  excitation  is  a 
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Figure  5.  One-Port  Cavity  Formed  from  a  Coaxial  Waveguide 
with  a  Step  Discontinuity  in  the  Outer  Radius 
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TEM  wave  traveling  in  the  positive  z  direction  in  the  region 
z  <  s.  Because  of  the  perfectly  conducting  plate  at  z  =  t, 
the  following  boundary  condition  is  imposed : 


E*(P,t)  =  0  . 


An  inspection  of  Eq. (24)  reveals  that  this  condition  is  en¬ 
forced  by  requiring 


co  e 


•iV  ♦  c'  e3V 


and 


C*  e'3 V  +  C-  e3 
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=  0 
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The  resulting  relationships  between  constants  are 


and 


rj2k+ 
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(45a) 

(45b) 


Substitution  of  Eq. (45a)  into  Eq. (38)  and  Eq. (45b)  into 
Eq.  (39)  yields  the  following  expressions  for  the  constants: 
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(47) 


An  integral  equation  for  the  unknown  aperture  electric 
field  is  again  obtained  by  constraining  the  transverse  mag- 
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netic  field  to  be  continuous  across  the  aperture.  The  fol¬ 
lowing  expression  is  obtained  by  equating  Eqs.(19)  and  (26) 
with  the  constants  defined  by  Eqs.(46)  and  (47): 
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Collecting  terms  yields  the  equation: 
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I.  Microwave  Circuit  Concepts 


As  described  in  Chapter  I,  essentially  all  early  treat¬ 
ments  of  waveguide  discontinuities  were  directed  toward  cal¬ 
culating  the  resulting  change  in  admittance  of  the  waveguide. 
With  the  assumption  of  lossless  waveguide  walls  and  dielec¬ 
trics,  these  admittance  changes  are  purely  reactive  in  na¬ 
ture.  The  susceptance  expressions  found  in  the  Waveguide 
Handbook  provide  excellent  data  for  comparison  with  results 
obtained  from  solutions  to  the  previously  derived  integral 
equations.  Since  such  solutions  yield  only  the  radial  com¬ 
ponent  of  electric  field  in  the  aperture,  a  small  amount  of 
additional  analysis  is  required  to  compute  values  of  suscep¬ 
tance  due  to  the  discontinuities. 

Consider  a  uniform  transmission  line  extending  to  z  =  ±<=° 
which  has  a  single  discontinuity  at  z  =  s.  Such  a  transmis¬ 
sion  line  is  depicted  in  Figure  6.  The  admittance  at  z  =  s 
looking  in  the  positive  z  direction  is  given  by 

YL  =  Yo  +  jB 

where  Y*  is  the  transmission  line  TEM  characteristic  admit¬ 
tance  for  z  >  s  and  B  is  the  TEM  susceptance  of  the  disconti¬ 
nuity.  The  characteristic  admittance  of  the  line  for  z  <  s 
is  designated  Yq.  The  complex  TEM  (voltage)  reflection  coef¬ 
ficient  in  the  region  z  <  s  is  defined  by 
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In  the  transmission  line  (or  waveguide)  region  z  <  s, 
the  TEM  voltage  reflection  coefficient  is  the  ratio  of  the 
reflected  TEM  voltage  (or  electric  field)  wave  to  the  inci¬ 
dent  TEM  voltage  wave.  As  noted  in  Section  D  of  this  chap¬ 
ter,  Bq  and  Bq  are  the  complex  coefficients  of  the  positive 
and  negatively  traveling  TEM  waves  for  Egs .  (17) - (19)  which 
are  defined  for  the  z  <  s  region.  The  reflection  coefficient 


can  now  be  represented  by 


T(z)  »  T(s)  ej2k-(z“s) 
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The  reflection  coefficient  can  be  defined  in  terms  of  the 
aperture  electric  field  through  the  use  of  Eqs. (34a&b) 


dp 


(53 


The  susceptance  of  the  discontinuity  may  then  be  calculated 
through  use  of  Eq. (50)  or  Eq. (51) .  It  should  be  noted  that 
Eq.  (53)  is  a  normalized  expression  since  E°  is  linearly  re¬ 


lated  to  V+. 


CHAPTER  III 


NUMERICAL  SOLUTION  TECHNIQUE 


A.  Introduction 

A  solution  for  the  unknown  aperture  electric  field  dis¬ 
tribution  in  Eq. (43 )  may  be  obtained  through  use  of  the 
method  of  moments.  Ref. [9].  Application  of  such  a  method 
requires  intelligent  choices  of  a  set  of  expansion  or  basis 
functions.  For  the  waveguide  discontinuities  being  treated 
in  this  investigation,  once  a  solution  for  the  aperture 
field  distribution  is  obtained,  all  other  quantities  of  in¬ 
terest  can  be  calculated.  The  desired  field  and  susceptance 
values  are  calculated  from  integrals  of  the  aperture  electric 
field  and  weight  functions  over  the  aperture.  The  sample 
solution  technique  illustrated  in  this  chapter  is  applicable 
to  all  single-aperture  equations. 

B.  Selection  of  Basis  and  Testing  Functions 

An  inspection  of  the  integral  under  the  summation  sign 
in  Eq. (43)  reveals  that  the  integral  is  easily  evalusted  if 
the  value  of  pE°  is  constant  over  the  range  of  integration. 
This  may  be  effected  by  the  choice  of  the  following  sub¬ 
sectional,  constant-value  pulse  expansion  for  the  unknown 
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aperture  field.  The  expansion  is  illustrated  in  Figure  7. 


P  E® (o) 


X  V0’ 


where 


1,  pe(r -A/2,r  +A/2) 
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1  0,  otherwise 


:54) 
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an<^  rq  =  a  +  (<3“)s)^  where  A  =*  (g-a)/N  is  the  pulse  width. 
Substituting  Eq. (54)  into  Eq.  (43)  yields 
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The  first  integral  is  easily  evaluated: 


(9  ?  1  N  r  q'“'‘ 
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N 
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(57) 


The  last  two  integrals  have  the  same  value  which  is  given 
by 
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=  MATCH  POINTS 


=  PULSE  WIDTH 
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After  the  substitution  of  the  integral  values  of  Eqs. (57) 
and  (58)  into  Eq. (56) ,  the  following  equation  is  obtained 
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An  important  consideration  in  the  selection  of  a  set  of 
testing  functions  is  the  avoidance  of  "sampling”  the  integral 
equation  at  a  point  where  the  aperture  field  has  a  singular¬ 
ity.  All  planar  discontinuities  discussed  in  the  report 
possess  90  degree  or  180  degree  exterior  corners  at  which 
points  the  normal  component  of  electric  field  is  singular. 

The  integral  equation  for  the  aperture  field  of  a  planar  co¬ 
axial  waveguide  discontinuity  such  as  given  by  Eq. (43)  pos¬ 
sesses  the  same  important  features  as  does  the  integral  equa¬ 
tion  for  the  two-dimensional,  conducting  strip  illuminated  by 
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a  TM  wave  (and  the  infinite  slot  illuminated  by  a  TE  wave) . 
The  experiences  of  Butler,  et  al. ,  Refs. [23]  and  [24],  sup¬ 
port  the  selection  of  constant-value,  pulse  basis  functions 
with  the  testing  functions  consisting  of  Dirac  delta  func¬ 
tions  at  the  pulse  centers.  Delta  functions  used  in  this 
manner  are  commonly  called  "match  points."  The  set  of  R 
match  points,  which  is  illustrated  in  Figure  7,  may  be  ex¬ 
pressed  mathematically  by 

T(p)  =  (6 (p-r1) , . . . ,5  (p-rp) , . . . ,6 (p-rr)|  .  (60) 

The  testing  functions  force  the  subject  equation  to  be  valid 
at  each  point 

p  *  rp  »  a  +  (p-^)A  ,  P  —  1,2, ...,R  . 
Substitution  of  the  expressions 


+  A/2 

=  a  +  qA  =  r+ 

(61a) 

-  A/2 

=  a  +  (q-l)A  *  r~ 

(61b) 

for  the  pulse  limits  in  Eq. (59)  simplifies  the  notation 
somewhat.  By  sampling  Eq. (59)  at  a  set  of  R  =  N  match 
points,  a  set  of  N  algebraic  equations  for  N  unknown  Pulse 
amplitudes  is  obtained  which  can  be  solved  for  the  unknown 
electric  field  in  the  aperture.  This  set  of  equations  is 
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(62] 


The  final  expression  given  by  Eq. (62)  is  now  in  a  de¬ 
sirable  form  for  numerical  determination  of  E^.  Special 
attention  must  be  paid  to  the  convergence  of  the  infinite 
series  of  Bessel  functions.  In  general,  such  series  have 
poor  convergence  properties  and  operations  on  the  series  to 
improve  the  rate  of  convergence  are  warranted.  A  very  ef¬ 
fective  procedure  for  improving  the  convergence  of  series 
such  as  that  in  Eq. (62)  is  described  in  Appendices  B  and  C. 
Application  of  the  procedure  does  not  affect  the  form  of 
the  equation  since  the  operations  are  only  upon  the  series , 
but  the  additional  step  of  determining  the  asymptotic  form 
of  the  series  elements  is  required. 


C.  Numerical  Results  for  Single  Discontinuities 

A  digital  computer  has  been  employed  to  obtain  solu¬ 
tions  to  Eq.(62)  with  constant  ratios  of  a,  b  and  c  at  five 
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different  frequencies.  The  real  and  imaginary  parts  of  the 
aperture  electric  field  are  plotted  for  a  frequency  range 
that  extends  from  below  to  above  the  cutoff  frequency  for 
the  first  higher-order  mode  in  the  waveguide.  The  first 
three  frequencies  are  below  and  the  last  two  are  above  cut¬ 
off.  All  plotted  field  distributions  display  a  ^  variation 
near  the  center  conductor  and  a  singular  behavior  as  the 
edge  of  the  disk  is  approached.  For  frequencies  above  the 
cutoff,  the  imaginary  part  normally  exhibits  one  or  more 
zero  crossings.  The  aperture  field  distributions  for  an 
annular  disk  attached  to  the  center  conductor  of  a  coaxial 
guide  (Figure  3b)  are  shown  in  Figure  9.  In  this  set  of 
curves,  the  singular  behavior  is  again  displayed  near  the 
edge  of  the  disk,  and  the  ^  variation  becomes  apparent  near 
the  outer  conductor. 

Similar  results  for  the  steps  in  outer  and  inner  radii 
of  the  coaxial  guides  illustrated  in  Figures  4a  and  4b  are 
presented  in  Figures  10  and  11.  The  behavior  of  the  aper¬ 
ture  field  distributions  is  analogous  to  that  shown  in  Fig¬ 
ures  8  and  9  with  the  singular  behavior  being  less  pro¬ 
nounced  near  the  90  degree  edges  than  for  the  180  degree 
ones.  In  Figure  12  are  shown  the  aperture  field  distribu¬ 
tions  obtained  when  a  shorting  plate  is  added  to  the  geom¬ 
etry  of  Figure  10  to  create  a  one-port  cavity. 
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In  Figure  13,  the  effect  of  varying  the  number  of 
basis  function  pulses  across  the  aperture  is  presented  for 
the  goemetry  of  Figure  3a.  The  plotted  curves  show  that 
the  correct  ~  behavior  is  captured  by  even  a  few  pulses 
across  the  aperture.  The  nature  of  the  field  singularity 
becomes  much  more  evident  as  the  number  of  pulses  is  in¬ 
creased.  Because  most  of  the  quantities  of  interest  in 
this  investigation  are  calculated  from  an  integral  of  the 
aperture  field,  the  effect  of  missing  the  fine  detail  of 
the  singularity  is  substantially  reduced,  a  fact  that  is 
demonstrated  in  Table  6 . 

In  Tables  1  through  5  are  presented  values  of  suscep- 
tance  computed  from  expressions  in  the  Waveguide  Handbook 
for  the  discontinuities  depicted  in  Figures  3a,  3b,  4a,  4b 
and  5.  Also  in  each  table  are  the  corresponding  susceptance 
values  calculated  from  Eq. (50)  and  (51)  with  the  aperture 
field  approximated  by  means  of  solutions  to  integral  equa¬ 
tions.  The  use  of  Eqs. (50  and  (51)  requires  evaluation  of 
the  aperture  field  integral  in  Eq.  (53) .  The  value  of  this 
integral  is  obtained  by  evaluating  the  right-hand  side  of 

Eq. (57)  with  the  known  values  of  E_. 

q 

The  variational  approach  used  in  the  derivation  of  the 
susceptance  expressions  in  the  Waveguide  Handbook  also  in¬ 
volves  determining  an  approximation  for  the  aperture  field. 
For  assumed  aperture  electric  field  distributions  that  dif- 
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fer  from  the  true  distribution  by  a  small  amount,  the  re¬ 
sulting  susceptance  computations  are  only  slightly  larger 
than  those  that  would  be  computed  from  the  true  distribution. 
It  should  be  recalled  from  Chapter  I  that,  with  the  varia¬ 
tional  approach,  a  first-order  error  in  the  aperture  distri¬ 
bution  results  in  only  a  second-order  error  in  the  calculat¬ 
ed  susceptance.  This  discussion  tends  to  indicate  that  the 
smaller  of  two  calculated  susceptance  values  is  the  more  ac¬ 
curate;  however,  caution  must  be  exercised  since  the  exact 
answer  is  not  known.  The  concept  of  the  true  susceptance 
value  being  a  local  minimum  is  supported  by  the  results  pre¬ 
sented  in  Table  6  where  complex  admittance  calculations  are 
given  for  three  variations  of  the  discontinuity  in  Figure  3a 
as  the  number  of  basis  function  pulses  is  increased.  Both 
the  real  and  imaginary  parts  of  the  admittance  appear  to  be 
converging  to  lower  bounds.  For  the  two  cases  where  the 
frequency  is  below  cutoff,  the  real  part  rapidly  tends  toward 
zero  as  the  number  of  pulses  is  increased.  As  the  number  of 
pulse  functions  is  increased,  there  is  a  requirement  for  a 
larger  number  of  terms  in  the  Bessel  function  series  to  main¬ 
tain  accuracy  because  representing  a  narrow  pulse  with  a 
Fourier-Bessel  series  requires  a  greater  number  of  terms 
than  does  a  wider  pulse.  It  is  probable  that  at  some  point 
numerical  error  will  become  significant,  but,  even  with  this 
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caveat ,  the  concept  of  the  lower  bound  on  calculated  suscep- 
tance  being  the  correct  value  is  supported. 

In  Tables  7  and  8,  the  effects  of  increasing  the  num¬ 
ber  of  terms  in  the  modified  higher-order  mode  series  are 
demonstrated  for  two  different  size  pulse  expansions  in 
three  different  apertures.  The  convergence  behavior  shown 
in  Table  6  is  not  present.  After  a  relatively  few  terms  in 
the  series  have  been  computed,  an  oscillatory  behavior  is 
observed  which  appears  to  converge  toward  an  ultimate  correct 
value.  It  is  expected  that  numerical  error  becomes  a  factor 
before  the  final  limit  is  identified,  but,  because  this  con¬ 
vergence  behavior  manifests  itself  in  relatively  few  terms, 
the  magnitude  of  such  error  appears  to  be  small. 

The  calculated  integral  equation  data  presented  in 
Tables  1  through  5  contain  calculations  made  at  the  cutoff 
frequency  for  the  TMq^  mode  in  the  larger  section  of  wave¬ 
guide.  These  calculations  are  actually  made  in  the  limit 
as  the  frequency  approaches  the  cutoff  point  and  are  courtesy 
of  the  small  inaccuracies  introduced  by  the  computer.  The 
integral  equations  have  terms  which  are  undefined  at  these 
transverse  resonance  frequencies.  This  attempt  to  sidestep 
the  singular  nature  of  the  equations  is  rewarded  by  negative 
real  parts  for  the  admittance  in  three  of  the  tables.  The 
relatively  large  real  parts  (although  very  small  when  com¬ 
pared  to  the  imaginary  parts)  of  the  admittances  calculated 


at  the  cutoff  frequencies  in  Table  5  are  also  a  product  of 
computer  error. 

It  is  not  possible  to  draw  conclusions  on  the  relative 
accuracies  of  the  integral  equation  and  Waveguide  Handbook 
approaches  without  establishing  parameters.  It  appears  to 
be  always  possible  to  take  sufficient  expansion  pulses  and 
series  terms  to  calculate  an  admittance  correction  which  is 
smaller  than  the  corresponding  Waveguide  Handbook  figure. 

All  calculations  made  for  Tables  1  through  5  and  for  Tables 
10,  11,  12,  and  13  used  20  pulses  and  30  series  terms.  This 
pulse-series  combination  gives  very  acceptable  results  with 
minimal  computer  time  expenditure. 

It  is  possible  to  compare  integral  equation  results  with 
more  accurate  calculations  of  coaxial  discontinuity  suscep- 
tance.  A.  Jurkus,  Ref. [13],  has  produced  a  refinement  of 
P.  I.  Somlo's  computer  code.  Ref. [10],  which  employs  the 
Hahn  series  approach  developed  by  Whinnery,  et  al . ,  Ref. [8]. 
Jurkus  presents  results  for  the  capacitance  of  a  coaxial 
waveguide  with  a  step  in  the  inner  conductor.  His  calcula¬ 
tions  contain  10  digits  and  the  numerical  precision  is  claimed 
to  be  ±1  in  the  fifth  digit.  In  Table  9,  comparisons  are  pre¬ 
sented  for  the  results  obtained  for  the  Jurkus,  integral 
equation,  and  Waveguide  Handbook  approaches.  The  integral 
equation  calculations  were  made  with  50  pulses  and  50  series 
terms.  Jurkus1  results,  which  employed  1000  series  terms 
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compared  to  Somlo's  40,  are  uniformly  slightly  less  than  0.3% 
below  the  integral  equation  results.  The  approach  employed 
by  Jurkus,  et  al . ,  does  not  involve  the  solution  of  a  set  of 
simultaneous  linear  equations,  so  there  can  be  no  straight¬ 
forward  comparison  of  numerical  computation  expense.  Ex¬ 
perience  gained  during  this  effort  strongly  indicates  that 
the  magnitude  of  the  computed  integral  equation  capacitance 
could  be  reduced  below  that  of  Jurkus  by  using  a  sufficient 
number  of  expansion  pulses.  It  should  be  noted  that  the 
agreement  of  the  integral  equation  capacitance  for  both 
geometries  was  obtained  by  different  calculations.  In 
Chapter  IX,  Table  14  presents  a  number  of  comparisons  of  the 
capacitance  of  a  coaxial-to-circular  waveguide  junction  in 
which  the  integral  equation  results  are  shown  to  be  very 
accurate. 
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TABLE  1 .  CALCULATED  SUSCEPTANCE  OF  AN  ANNULAR  DISK  ATTACHED  TO  THE 
OUTER  CONDUCTOR  OF  A  COAXIAL  WAVEGUIDE  (NORMALIZED) 


c/X 

a  ■  c/6 

Wavacuide  Handbook 

b  *  c/2 

Integral  Equation 

O.OOIO 

0.0  +  jO. 00156520 

0.00058660 

♦ 

jO. 00157331 

0.0100 

0.0  +  jO. 01565200 

0.00058659 

♦ 

jO. 01565393 

0.0500 

0.0  ♦  jO. 07850436 

0.00058634 

jO. 07890602 

0.1000 

0.0  +  jO. 15850730 

0.00058555 

+ 

jO. 15928538 

0.2000 

0.0  +  jO. 33006177 

0.00058209 

+ 

jO. 33135521 

0.3000 

0.0  +  jO. 53492015 

0 . 000S7511 

jO. 53571828 

0.4000 

0.0  +  jO. 82068523 

0.00056125 

jO. 31659679 

0.5000 

0.0  +  jl. 38924230 

0.00052592 

+ 

jl. 34509956 

0.5600 

0.0  +  j2. 82512207 

0.00043400 

■f 

j2. 38070170 

0.5817 

Cutoff  Freq 

-0.00005884 

+ 

j7. 67352323 

0.6000 

Above  Cutoff 

3.11703197 

+ 

jl. 61050941 

0.7000 

Above  Cutoff 

1 . 6441C242 

+ 

jO. 38996023 

a  *  c/2  b  *  2c/3 


C/X 

Wavacuide  Handbook 

Integral  ) 

Equation 

0.0010 

0.0 

♦ 

jO. 00197451 

0.00003314 

♦ 

jO. 00199170 

0.0100 

0.0 

♦ 

jO. 01974593 

0.00003314 

+ 

jO. 01991780 

0.0500 

0.0 

+ 

jO. 09882843 

0.00003314 

♦ 

jO. 09968915 

0.1000 

0.0 

+ 

jO. 19827879 

0.00003313 

jO. 20000873 

0.2000 

0.0 

+ 

jO. 40167294 

0.00003311 

jO. 40520104 

0.4000 

0.0 

jO. 84941823 

0.00003302 

jO. 85702615 

0.5000 

0.0 

jl. 11344261 

0.00003293 

♦ 

jl. 12347096 

0.6000 

0.0 

jl. 42841388 

0.00003280 

+ 

jl. 44118492 

0.8000 

0.0 

♦ 

j2. 43356469 

0.00003227 

4* 

j2. 45173828 

0.9000 

0.0 

j3. 61737121 

0.00003151 

♦ 

j3. 62722902 

0.9945 

Cutoff  Freq 

0.00001057 

■f ; 

)30. 33650604 

1.1000 

Above  Cutoff 

4.15095650 

■f 

jO. 85150308 
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.  CALCULATED  SUSCEPTANCE  OF  AN  ANNULAR  DISK  ATTACHED  TO  THE 
INNER  CONOUCTOR  OF  A  COAXIAL  WAVEGUIDE  (NORMALIZED) 


T 


c/X 

a  ■  c/6 

Wavecuide  Handbook 

b  -  c/2 

Integral  Equation 

O.OOIO 

0.0  +  jO.  00362468 

0.00027253 

* 

jO. 00365967 

0.0100 

0.0  +  jO. 03625057 

0.00027254 

jO. 03660017 

0.0500 

0.0  +  jO. 18170876 

0.00027267 

+ 

jO. 18342761 

0.1000 

0.0  +  jO. 36631767 

0.00027309 

jO. 36956776 

0.2000 

0.0  +  jO. 75753951 

0.00027487 

jO. 76232629 

0.3000 

0.0  +  jl. 21013329 

0.00027839 

jl. 21116118 

0.4000 

0.0  *  jl. 80461695 

0.00028501 

+ 

jl. 78144529 

0.5000 

0.0  +  j2. 88524378 

0.00029977 

+ 

j2. 67602259 

0.S600 

0.0  +  jS. 72519939 

0 . Q0032S87 

4* 

33.85053889 

0.5817 

Cutoff  Freq 

0.00038165 

+ 

j5. 81638213 

0.6000 

Above  Cutoff 

2.15903638 

+ 

j5. 12827991 

0.7000 

Above  Cutoff 

3.62498163 

j4. 91082115 

c/X 

a  -  c/2 

Wavecrui.de  Handbook 

b  -  2c/3 

Integral  Equation 

0.0010 

0.0  +  jO. 00083319 

0.00007785 

* 

jO. 00085388 

0.0100 

0.0  +  jO. 00833212 

0.00007785 

♦ 

jO. 00853904 

0.0500 

0.0  ♦  jO. 04168796 

0.00007786 

♦ 

jO. 04272324 

0.1000 

0.0  +  jO. 08354763 

0.00007787 

♦ 

jO. 08562268 

0.2000 

0.0  +  jO. 16849223 

0,00007795 

♦ 

jO. 17268213 

0.4000 

0.0  +  jO. 34899602 

0.00007648 

♦ 

jO. 35781976 

0.5000 

0.0  «•  jO. 44883165 

0.00007856 

♦ 

jO. 46044428 

0.6000 

0.0  +  jO. 55951281 

0.00007895 

jO. 57451142 

0.9000 

0.0  +  jO. 34654291 

0.00008040 

♦ 

jO. 86911853 

0.9000 

0.0  +  jl. 08977905 

0.00008196 

+ 

jl. 09639242 

0.9945 

Cutoff  Freq 

0.00008834 

♦ 

jl. 70338315 

1.1000 

Above  Cutoff 

0.71330313 

jl. 61754223 
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TABLE  3.  CALCULATED  SUSCEPTANCE  OF  A  STEP  IN  THE  OUTER  RADIUS  OF  A 
COAXIAL  WAVEGUIDE  (NORMALIZED  BY  THE  LEFT-HAND  WAVEGUIDE) 


i 

a  »  c/6  b  -  c/2 


c/X 

Wavequide  Handbook 

Integral  Equation 

0.00X0 

0.0  ♦  jO. 00052259 

0.00051126 

jO. 00052505 

0.0100 

0.0  +  jO. 00522663 

0.00051126 

+ 

jO. 00525117 

o.osoo 

0.0  ♦  jO. 02621644 

0.00051116 

jO. 02633891 

0.1000 

0.0  +  jO. 05296478 

0.00051083 

jO. 05320808 

0.2000 

0.0  ♦  jO. 11056802 

0 .00050942 

+ 

jO. 11103353 

0.3000 

0.0  ♦  jO.  18006588 

0.00050655 

+ 

jO. 18064396 

0.4000 

0.0  +  jO. 27869199 

0.00050074 

jO. 27879113 

o.sooo 

0.0  +  jO. 47982823 

0.00048527 

♦ 

jO. 47401488 

0 . 5600 

0.0  +  jO. 98287410 

0.00044321 

-4* 

jO. 90923393 

0.5817 

Cutoff  Freq 

-0.00007031 

♦ 

j5. 68891810 

0.6000 

Above  Cutoff 

1.16333598 

♦ 

30.26601990 

0.7000 

Above  Cutoff 

0.53058243 

30.07136124 

a  *  c/2 

b  »  2c/3 

c/x 

wavequi.de  Handbook 

Integral 

Equation 

0.0010 

o.o  +  jO. 00044256 

0.00002425 

•4» 

30.00044315 

0.0100 

0.0  +  jO. 00442577 

0.00002425 

♦ 

30.00443165 

0.0500 

0.0  ♦  jO. 02215205 

0.00002425 

30.02218148 

0.1000 

0.0  +  jO. 04445010 

0.00002425 

♦ 

30.04450931 

0.2000 

0.0  +  jO. 09010117 

0.00002424 

30.09022265 

0.4000 

0.0  *  jO. 19102869 

0.00002420 

+ 

30.19129943 

0.5000 

0.0  ♦  jO. 25094968 

0.00002417 

♦ 

30.25131700 

0.6000 

0.0  jO. 32292430 

0 . 0Q002413 

♦ 

30.32340570 

0.3000 

0.0  +  jO. 55683513 

0.00002394 

+ 

30.55737816 

0.9000 

0.0  *  jO. 84136428 

0.00002366 

30.84018797 

0.9945 

Cutoff  Freq 

0.00000713 

♦314.87807544 

1.1000 

Above  Cutoff 

0.90926887 

♦ 

30.12537772 

TABLE  4 


CALCULATED  SUSCEPTANCE  OF  A  STEP  IN  THE  INNER  RADIUS  OF  A 
COAXIAL  WAVEGUIDE  (NORMALIZED  BY  THE  LEFT-HAND  WAVEGUIDE) 


c/X 

a  «  c/6 

wavacuida  Handbook 

b 

-  c/2 

Intaqral  Equation 

O.OOIO 

0.0  +  jO. 00209161 

0.00044132 

+  jO. 00210399 

0.0100 

0.0  +  jO. 02091867 

0.00044132 

+  jO. 02104239 

0.0500 

0.0  +  jO. 10490495 

0.00044144 

+•  jO. 10551809 

0.1000 

0.0  >  jO. 21179806 

0.00044182 

+  jO. 21298804 

0.2000 

0.0  +  jO. 44086145 

0.00044349 

+■  jO. 44287046 

0.3000 

0.0  +  jO. 71368241 

0.00044685 

+  jO. 71510837 

0.4000 

0.0  +  jl. 09143147 

0 . 000453S6 

♦  jl. 08610510 

O.SOOO 

0.0  *  jl. 82817276 

0.00047050 

+  jl. 76514962 

0.5600 

0.0  +  j3. 63777780 

0.00050965 

+  j2. 99413580 

0.5817 

Cutoff  Fraq 

0.00067809 

+  j7. 50582984 

0.6000 

Abova  Cutoff 

3.70562650 

♦  j2. 80386573 

0.7000 

Abova  Cutoff 

2.35142230 

+  jO. 69252237 

c/X 

a  ■  c/2 

wavaouida  Handbook 

b 

-  2c/3 

Intaqral  Equation 

0.0010 

0.0  *  jO. 00050401 

0.00010060 

+  jO. 00050544 

0.0100 

0.0  +  jO. 00504031 

0.00010060 

-  jO. 00505455 

0.0500 

0.0  +  jO.  02522438 

0.00010060 

+  jO. 02529554 

0.1000 

0.0  +  jO. 05059232 

0.00010062 

+  jO. 05073454 

0.2000 

0.0  +•  jO. 10236195 

0.00010067 

+  jO. 10264547 

0.4000 

0.0  «•  jO. 21519440 

0.00010092 

+  jO. 21574369 

0.5000 

0.0  +  jO. 28051462 

0.00010113 

♦  jO. 28115931 

0.6000 

0.0  +  jO. 35676104 

0.00010145 

♦  jO. 35739963 

0.8000 

0.0  +  jO. 58432744 

0.00010275 

♦  jO. 58288265 

0.9000 

0.0  ♦  jO. 82080969 

0.00010446 

■f  jO. 80792331 

0.9945 

Cutoff  Fraq 

0.00011939 

+  j2. 24993063 

1.1000 

Abova  Cutoff 

1.06489802 

+  jO. 63059216 
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TABLE 


5.  CALCULATED  SUSCEPTANCE  OF  A  ONE-PORT  COAXIAL  CAVITY 
(NORMALIZED  BY  THE  INPUT  WAVEGUIDE) 


J 


c/X 

a  -  c/6 

Waveguide  Handbook 

b 

-  c/2  d 

Integral 

■  3c 

.  Equation 

O.OOIO 

0.0  -  j32. 52409365 

0.00031693 

- 

j  32 . 53439935 

0.0100 

0.0  -  j  3.20900351 

0.00031693 

- 

j  3.20999767 

0.0500 

0.0  -  j  0.41926107 

0.00031687 

- 

j  0.41927976 

0.1000 

0.0  +  j  0.25213838 

0.00031667 

+ 

j  0.25249477 

0.2000 

0.0  -  j  0.73335669 

0.00031579 

- 

j  0.73315769 

0.3000 

0.0  +  j  1.02399059 

0.00031401 

+ 

j  1.02483368 

0.4000 

0.0  +  j  0.07946839 

0.00031041 

j  0.07950574 

0.5000 

0.0  -  j  4.9826x10* 

881.54673050 

- 

j  4.9841x10* 

0 . 5600 

0.0  +  j  0.69434893 

0.00027447 

j  0.62456391 

0.5817 

Cutoff  Freq 

-0.00004352 

♦ 

j  5.67004190 

0.6000 

Above  Cutoff 

0.00020467 

♦ 

3  2.26837281 

0.7000 

Above  Cutoff 

0.00035314 

- 

3  1.12708972 

c/X 

a  «  c/2 

Waveguide  Handbook 

b 

»  2e/3 

Integral  Equation 

0.0010 

0.0  -  j22. 01537281 

0.00001714 

- 

322.01574949 

0.0100 

0.0  -  j  2.17127699 

0.00001714 

- 

3  2.17130840 

0.0500 

0.0  -  j  0.27939034 

0.00001714 

- 

3  0.27936609 

0.1000 

0.0  +  j  0.17930395 

0.00001713 

♦ 

3  0.17936609 

0.2000 

0.0  -  j  0.48114894 

0.00001713 

- 

3  0.48103724 

0.4000 

0.0  ♦  j  0.05617483 

0.00001711 

3  0.05644327 

0.5000 

0.0  -  j  3.3727x10* 

403.14432740 

- 

3  3.3728x10* 

0.6000 

0.0  +  j  0.45777816 

0.00001705 

♦ 

3  0.45826186 

0.3000 

0.0  +  j  1.12803524 

0.00001692 

4* 

3  1.12863794 

0. 9000 

0.0  ♦  j  0.70651042 

0.00001672 

■4 

3  0.70533203 

0.9945 

Cutoff  Freq 

0.00000501 

♦ 

318.90491948 

1.1000 

Above  Cutoff  ■ 

0.00001645 

3  1.53150570 
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TABLE  6 


CALCULATED  SUSCEPTANCE  OF  AN  ANNULAR  DISK  ON  THE  OUTER 
CONDUCTOR  OF  A  COAXIAL  WAVEGUIDE  AS  THE  NUMBER  OF 
EXPANSION  PULSES  IN  THE  APERTURE  IS  VARIED  (NORMALIZED) 


a  -  c/5  b  ■  4e/5  X  ■  2c  x cutoff  “  1>65c 


Waveguide 

Handbook  Susceptance  »  0.0  +  iO. 16470255 

No.  of 

A 

Integral  Equation 

Pulses 

X 

Susceptance 

5 

0.0600 

0.01727536  +  jO. 18367030 

10 

0.0300 

0.00446745  +  jO. 16502277 

IS 

0.0200 

0.00199920  +  jO. 15899635 

20 

0.0150 

0.00112726  +  jO. 15601510 

25 

0.0120 

0.00072224  +  jO. 15423484 

30 

0.0100 

0.00050185  +  jO. 15305148 

35 

0.0086 

0.00036883  +  jO. 15220771 

40 

0.0075 

0.00028244  +  jO. 15157575 

50 

0.0060 

0.00018081  +  jO. 15069203 

a  »  c/5 

b 

-  2c/5 

X  "  2c  X cutoff  "  l-65c 

Waveguide 

Handbook  Susceptance  -  0.0  +  12.27553677 

NO.  of 

A 

Integral  Equation 

Pulses 

X 

Susceptance 

5 

0.0200 

0.00307077  ♦  32.32601738 

10 

0.0100 

0.00076432  f  j2. 28872639 

15 

0.0067 

0.00033898  +  32.27678771 

20 

0.0050 

0.00019045  +  32.27091264 

25 

0.0040 

0.00012180  +  j2. 26741777 

30 

0.0033 

0.00008454  +  32.26510059 

35 

0.0029 

0.00006209  +  j2.26345157 

40 

0.0025 

0.00004752  +  j2.26221827 

50 

0.0020 

0.00003040  +  32.26049647 

a  -  c/5 

b 

-  2c/5 

X  -  X cutoff  -  1*65c 

Waveguide  Handbook  Susceptance  »  Above  Cutoff  Freq 


NO.  Of 

A 

Integral  Equation 

Pulses 

X 

Susceptance 

5 

0.2000 

2 . 63616915 

30.18763705 

10 

0.1000 

2.59431266 

♦ 

30.17881963 

15 

0.0667 

2.58133988 

+ 

30.17437867 

20 

0.0500 

2.57510034 

+ 

30.17181552 

25 

0.0400 

2.57144336 

♦ 

30.17015885 

30 

0.0333 

2.56904400 

♦ 

30.16900265 

35 

0.0286 

2.56734976 

♦ 

30.16815072 

40 

0.0250 

2.56609013 

♦ 

30.16749720 

50 

0.0200 

2.56434298 

♦ 

30.16656094 

«  width  of  Puis*  30  terms  were  calculated  in  the  series 
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TABLE  7.  CALCULATED  SUSCEPTANCE  OF  AN  ANNULAR  DISK  ATTACHED  TO  THE  OUTER 
CONDUCTOR  OF  A  COAXIAL  WAVEGUIDE  AS  THE  NUMBER  OF  TERMS  IN  THE 
HIGHER-ORDER  MODE  SERIES  IS  VARIED  (NORMALIZED) 

APERTURE  FIELD  EXPANDED  IN  20  PULSES 


4  •  c/S 
20c  X 


4c/5 


cutoff 


1.65c 


4  ■  e/5 

20c  X 


cutoff 


2e/S 
•  1.65c 


Wavacuida  Handbook  Suacaptanca  Ha would#  Handbook  Suacaptanca 


0.0  *  JO. 01280232 


0.0  ♦  JO. 14230349 


a  •  c/5  b  ■  2e/5 
X  "  C/S  X cutoff  ‘  1-S5c 

Hawculda  Handbook  Suacaptanca 
Abow  Cutoff  Fraquancy 


do.  Tana  Intaqral  Equation 
In  Saslaa  Suacaptanca 


No.  Ta sms  Intaqral  Equation 

In  Saslaa  Suacaptanca 


No.  Tarn*  Intaqral  Equation 
In  Saslaa  Suacaptanca 


1 

0.00118614 

* 

JO. 01351677 

1 

0.00020059 

♦ 

JO. 15513573 

1 

1.37815686 

♦ 

J4. 44734382 

5 

0.00116467 

JO. 01319154 

5 

0.00019970 

+ 

JO. 14343800 

5 

2.43647412 

♦ 

JO. 37385966 

10 

0.00116489 

JO. 01314581 

10 

0.00019960 

♦ 

JO. 14332278 

10 

2.57052629 

♦ 

JO. 16619757 

15 

0.00116502 

+ 

JO. 01314133 

IS 

0.00019952 

♦ 

JO. 14334945 

15 

2.57297198 

♦ 

JO. 17056063 

20 

0.00116507 

♦ 

JO. 01314024 

20 

0.00019956 

JO. 14334071 

20 

2.57426588 

♦ 

JO. 17093339 

25 

0.00116508 

♦ 

JO. 01313959 

25 

0.00019955 

♦ 

JO. 14334229 

25 

2.57486588 

♦ 

JO. 17160180 

30 

0.00116507 

♦ 

JO. 01313960 

30 

0.00019955 

♦ 

J0.14334152 

30 

2.57510034 

♦ 

JO. 17181552 

40 

0.00116506 

♦ 

JO. 01313975 

40 

0.00019955 

■4* 

JO. 14334104 

40 

2.57524851 

-f 

JO. 17171150 

SO 

0.00116506 

* 

JO. 01313969 

SO 

0.00019955 

♦ 

JO. 14334051 

SO 

2.57532391 

♦ 

JO. 17195233 

TABLE  8.  CALCULATED  SUSCEPTANCE  OF  AN  ANNULAR  DISK  ATTACHED  TO  THE  OUTER 
CONDUCTOR  OF  A  COAXIAL  WAVEGUIDE  AS  THE  NUMBER  OF  TERMS  IN  THE 
HIGHER-ORDER  MOOE  SERIES  IS  VARIED  (NORMALIZED) 

APERTURE  FIELD  EXPANDED  IN  50  PULSES 


a  »  e/5 
20c  X 


4c/ 5 


cutoff 


1.65c 


a  •  c/5 
X  -  20c  X 


b  • 

cutoff 


2c/ 5 
•  1.65c 


a  -  c/5 
c/5  X 


b  •  2c/5 


cutoff 


•  1.65c 


wawqulda  Handbook  Suacaptanca 
0.0  *  JO. 01280232 

No.  Tasaa  Intaqral  Equation 
In  Saslaa  Suacaptanca 

1  0.00018983  *  JO. 01310735 

5  0.00018667  *  JO. 01278999 

10  0.00018669  ♦  JO. 01274783 

IS  0.00018671  *  JO. 01274350 

20  0.00018671  *  JO. 01274270 

25  0.00018672  ♦  JO. 01274224 

30  0.00018672  *  JO. 01274206 

40  0.00018671  *  JO. 01274218 

50  0.00018671  *  JO. 01274213 


Hawquida  Handbook  Suacaptanca 
0.0  +  J  .14230349 

No.  Tasaa  Intaqral  Equation 
in  Saslaa  Suacaptanca 

1  0.00003203  *  JO. 15424329 

5  0.000031B9  *  JO. 14273483 

10  0.00003186  *  JO. 14261662 

15  0.00003186  *  JO. 14264501 

20  0.00003187  *  JO. 14263675 

25  0.00003187  *  JO. 14263748 

30  0.00003187  *  JO. 14263699 

40  0.00003187  *  JO. 14263666 

50  0.00003187  *  JO. 14263606 


Hawquida  Handbook  Suacaptanca 
Abow  Cutoff  Fraquancy 


No.  Tana  Intaqral  Equation 
in  Saslaa  Suacaptanca 


1 

1.37399736 

♦ 

J4. 39469120 

S 

2.42168653 

♦ 

JO. 367 07202 

10 

2.55995962 

♦ 

JO. 16163260 

IS 

2.56233020 

♦ 

JO. 16542236 

20 

2.56360487 

♦ 

JO. 16577117 

25 

2.56413201 

♦ 

JO. 16636831 

30 

2.56434298 

♦ 

JO. 16656094 

40 

2.56447746 

♦ 

JO. 16663587 

50 

2.56453807 

♦ 

JO. 16666803 
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TABLE  9.  COMPARISONS  OF  CALCULATED  CAPACITANCE  OF  COAXIAL 
STEP  DISCONTINUITIES  WITH  RESULTS  OF  A.  JURKUS 


2b 

r 

1  1  : 

2c 

_ 

2*  I 

tz 

i 

PT _ : 

a  ■  0.209672  cm 

b  -  0.714375 

cm  c  *  0.310205  cm 

Freq 

GHz 

JurJcus 

Capacitance 

femtoFarads 

Intaqral  Equation 

Complex  Capacitance 
femtoFarads 

Waveguide  Handbook 
Capacitance 
femtoFarads 

0.000 

18.87112905 

0.0006691302 

+  J18. 92582501* 

•• 

1.000 

18.87742237 

0.0000000008 

+  jl8. 93211996 

18.96409053 

2.000 

18.89635437 

0.0000000003 

+•  jlB.  95105772 

18.98333401 

3.000 

18.92808243 

0.0000000002 

+  jl8. 98279505 

19.01559548 

4.000 

18 . 97287338 

0.0000000002 

+  j!9. 02759828 

19.06116396 

—  ——A*  1 

2a 

* 

::  2c 

_ 3 

’ 

a  -  0.310205  cm 

b  -  0.714375 

cm  c  *  0.209672  cm 

Freq 

GHz 

JurJcus 

Capacitance 

femtoFarads 

Integral  Equation 
Complex  Capacitance 
femtoFarads 

Waveguide  Handbook 
Capacitance 
femtoFarads 

0.000 

18.87101768 

0.0006691302 

+  jl8. 92582501* 

** 

1.000 

18.87731097 

0.0000000006 

+  jl8. 93211996 

18.96409053 

2.000 

18.89624288 

0.0000000004 

+  JIB. 95105772 

18.98333401 

3.000 

18.92797080 

0.0000000002 

+  J18. 98279505 

19.01559548 

4.000 

18.9727-6154 

0.0000000001 

+  J19. 02759828 

19.06116396 

Integral  Equation  calculation  at  1000  Hz 
**Ho  Raault 


CHAPTER  IV 


INTEGRAL  EQUATIONS  FOR  A  TWO-PORT  COAXIAL  CAVITY 

A.  Introduction 

In  Chapter  II,  integral  equations  are  derived  for  the 
transverse  electric  field  in  the  apertures  of  the  coaxial 
waveguide  discontinuities  depicted  in  Figures  3  and  4 .  Sub¬ 
sequently,  it  is  demonstrated  that  a  one-port  coaxial  cavity 
can  be  formed  by  adding  a  shorting  plate  to  a  waveguide  con¬ 
taining  a  discontinuity,  and  the  integral  equation  appropri¬ 
ate  for  such  a  geometry  is  derived.  The  analysis  is  now  un¬ 
dertaken  for  the  treatment  of  a  two-port,  coaxial  cavity 
formed  Ly  placing  two  discontinuities  in  a  coaxial  waveguide. 
Such  a  geometry  is  illustrated  in  Figure  14  where  three  re¬ 
gions  of  interest  are  defined: 

Region  1:  z  <  s 
Region  2:  s  <  z  <  t 
Region  3 :  z  >  t  . 

B.  Expressions  for  Field  Components 

In  order  to  describe  the  field  in  each  of  the  three  re¬ 
gions,  three  sets  of  expressions  equivalent  to  those  of 
Eqs.  (17) - (22)  are  required.  The  following  set  of  equations 
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represents  the  field  in  each  of  the  three  regions  where 
individual  regions  are  specified  by  i  =  1,  2,  or  3. 


i  +  e“3k; 2  _ 

E1  =  B  - - —  +  B._  - — — 

p  lO  p  lO  p 


+  l  <bJ 

n=l 


-ict.  z  ,  ]a.  z,  d<&  . 

e  J  in  +  B.  eJ  m  )  m 
in  m  - 


dp 


2 

oo  v  , 

y  in 

(B+  e“ 

jctin2  - 

,»1  ain 

Bin  e 

xn 

e'^i2 

'  -  BIo 

ejv 

i  niP 

nip 

k .  00 

—  T 

1  ,n  + 

2 

l  n— l 

ain  in  = 

B'n  e^in2)  *in 
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where 


*i.nlp)  *  -  V^in^'W 

and 

2  2  2 
in  *  kin  -  ^in 

The  are  the  roots  of  the  transcendental  equations: 

4lnlb>  *  N0(1'lna>J0(lrlnb)  *  J0(Ylna)  VYlnb)  "  0  ' 
«2„<c)  =  VY2na,VY2nc)  -  VY2na)VY2nc>  -  0  , 

*3n<d)  -  N0(Y3nal J0(,f3ndl  *  J0(Y3na)  VY3ndl  =  0 


the 

(63) 

(64) 

(65) 

(66) 

(67) 

(68a) 

(68b) 

(68c) 


where  a,  b,  c  and  d  are  radial  dimensions  defined  in 


Figure  14.  The  dielectric  medium  in  each  of  the  three  re¬ 
gions  is  characterized  by  the  parameters: 


and 


(69a) 


(69b) 


C.  Evaluation  of  Constants  in  Field  Expressions 

The  procedure  for  evaluating  the  expansion  constants 
in  Eqs.(63),  (64)  and  (65)  for  i  =  1,  2  and  3  is  the  same 
as  that  employed  in  Chapter  II.  The  unknown  aperture  field 
in  each  of  the  two  discontinuities  is  designated 

Ep  (p  ,  s)  =  Ep(p,s)  *  E*  (70a) 

and 

E*(p,t)  *  E^(p,t)  =  zl  .  (70b) 

The  above  notation  allows  the  unknown  constants  in  all 
three  regions  to  be  expressed  in  terms  of  weighted  inte¬ 
grals  over  the  appropriate  aperture  fields.  A  notation  is 
adopted  for  representing  the  individual  integrals  which 
will  allow  the  resulting  equations  to  be  written  more  com¬ 
pactly.  Each  integral  will  be  denoted  by  1^  where  i  is 

Sr 

the  number  of  the  aperture,  p  =  0  represents  a  weight  func¬ 
tion  of  unity,  and  p  =  mn  indicates  the  use  of  the  proper 
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Bessel  function  weighting  factor.  The  following  expressions 
are  obtained  for  the  constants  in  each  region  by  application 
of  the  Fourier-Bessel  series  techniques  described  in  Chapter 
II. 


Region  1: 


B*  Q  e'jkls  +  8^  Q  ejkls 


dp  ■  3 
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Region  2 : 
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(72a) 


t-4 

2n 


(73a) 


»Jn  e’»2n* 


i  r 

m2  J, 


E l  o  d$2n  dp  =  -i-  (73b) 

T  «2„  2" 


*  The  subscript  notation  1,0  is  used  where  there  is  a  chance 
of  misinterpretation;  otherwise,  the  comma  is  omitted  in  the 
subscripts  of  the  unknown  constants  and  integral  symbols. 
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Region  3 : 


B3,0  6 
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.2 


El  P  d$3n  dp  =  -i-  1^  .  (74b) 
M-n  1  a  —  M3n  3" 


In  Region  1,  the  constants  for  the  waves  traveling  in 
the  negative  z  direction  may  be  expressed  in  terms  of  the 
constants  for  the  positive-traveling  waves: 


e-3kls 


e-%s 


e'^ln3 
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In 


(75a) 


(75b) 


In  Region  2  there  are  four  unknown  constants  and  also  four 
equations  given  by  Eqs. (69)  and  (70) .  It  is  therefore  pos¬ 
sible  to  express  each  constant  in  terms  of  aperture- field 
integrals : 


2 , 0 


j2sink-(t-s)  ln^- 

4  d 


tii  e±jk: 


+  ij  e±jk2J 


(76a) 
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2n 


j2sina2n(t-s)  M2n 


2n  •****'  ?  4 


(76b) 


In  Region  3,  the  constants  for  the  waves  traveling  in  the 
positive  z  direction  are  expressed  in  terms  of  the  constants 
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for  the  negative-traveling  waves: 
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D.  Enforcing  Continuity  of  Transverse  Fields  in  Apertures 
The  next  step  in  formulating  an  integral  equation  for 
the  transverse  electric  field  in  each  of  the  two  apertures 
is  to  require  that  the  transverse  component  of  the  magnetic 
field  be  continuous  across  each  aperture.  The  continuity 
requirement  on  the  transverse  electric  field  already  is  en¬ 
forced  by  Eqs . (70a&b) .  An  equation  is  now  obtained  by  re- 
1  2 

quiring  that  H^(p,s)  »  H^fPrS): 
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Collecting  the  terms  with  unknown  constants  on  the  left-hand 
side  yields  the  equation 


2R+  ^kis 

2B1,0  e  1 

n1P 


'1  n=l 


-2-  bT  e-J“ins 
In  ln  “3o 


In 


.  b 
nlp  lnI 


k.  <*> 

+  —  l 

nl  n=l 


a,  M? 

ln  ln 


rl  d$ 
In 


ln 


dp 


.  jsink2 (t-s)  n2P 


In* 


cosk2 (t-s) 


I 


2) 

0 


3  i  _ i _ 

n2  n-1  jsina2n(t-s)  a2nM2n 


*2n  <=°sa2n(t-s)  -  I 


2n 


d$ 


2n 


dp 


(78) 


Inspection  of  Eq. (78)  reveals  that  this  equation  would 

•  <  «f* 

be  identical  in  form  to  Eq.  (48)  if  all  the  constants  Bln  and 
.  2  2 

the  integrals  Ig  and  I2n  were  equal  to  zero.  The  reason  for 
this  similarity  is  that  Eq. (78)  is  the  general  form  of  the 
coupled  integral  equations  obtained  when  two  planar  discon¬ 
tinuities  are  placed  in  a  coaxial  waveguide.  A  shorting 
plate  is  a  special  case  of  a  planar  discontinuity  and,  there¬ 
fore,  Eq.  (48)  is  representative  of  a  two-discontinuity  geom¬ 
etry.  It  follows  that  the  integral  equation  obtained  from 
constraining  H^(p,t)  to  be  equal  to  H^(p,t)  has  the  same 
form  as  Eq. (78)  and  can  be  written  as 
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If  the  constants  Q,  B^n,  B3  Q  and  B~n  are  known, 


Eqs . (78 )  and  (79)  represent  two  integral  equations  for  the 
two  unknown  aperture  field  distributions.  If  the  output 
waveguide  is  matched  and  there  is  no  excitation  incident 
from  the  negative  z  direction,  B3  g  and  all  B3n  will  be 
zero.  The  numerical  solution  technique  described  in  Chapter 
III  may  be  applied  to  each  equation  in  order  to  obtain  ex¬ 
pansions  of  the  form  of  Eq.(54)  for  each  aperture  electric 
field. 

Aperture  electric  field  distributions  have  been  calcu¬ 
lated  for  two  different  size  two-port  cavities  at  two  fre¬ 
quencies  for  each  cavity  and  the  results  are  presented  in 
Figures  15  and  16.  The  two  cavities  that  are  treated  are 
symmetric  in  that  the  radii  b  and  d  as  depicted  in  Figure  14 
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are  equal.  The  two  frequencies  are  chosen  such  that  one  is 
below  and  the  other  is  above  the  cutoff  frequency  for  the 
first  higher-order  mode  in  the  larger  section  of  wave¬ 
guide.  The  curves  in  Figures  15  and  16  display  more  variety 
than  the  aperture  distributions  shown  in  Figures  11  and  12 
which  are  for  similar-size  step  discontinuities.  The  distri¬ 
butions  in  Figures  15  and  16  do,  however,  possess  the  same 
characteristics  of  decreasing  in  magnitude  as  i-  at  points 
near  the  center  conductor  and  then  increasing  without  bound 
as  the  discontinuity  is  approached.  In  contrast  to  the 
aperture  distributions  for  the  one-port  cavity,  there  are 
examples  of  real  and  imaginary  components  crossing  the  zero 
axis. 

The  approximation  for  the  aperture  field  at  z  =  s 
allows  the  calculation  of  the  input  admittance  of  the  cavity 
at  that  point  through  the  use  of  Eqs. (50)  and  (53)  .  A 
similar  expression  may  be  obtained  for  the  input  admittance 
looking  into  the  cavity  at  z  =  t  by  accounting  for  the  direc¬ 
tions  of  wave  propagation.  Calculated  values  of  input  admit¬ 
tance  are  presented  in  Table  10  for  two  different  size 
cavities  as  a  function  of  frequency.  The  admittances  are 
computed  both  from  solutions  to  Eqs.  (78)  and  (79)  and  from 
transmission  line  theory  with  susceptance  corrections  from 
an  expression  in  the  Waveguide  Handbook .  The  latter  cal¬ 
culations  are  valid  only  for  frequencies  below  the  cutoff 
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frequency  of  the  first  higher-order  mode  in  the  larger  cav¬ 
ity.  Comparisons  between  calculated  admittance  values  and 
measured  ones  for  a  two-port  coaxial  cavity  are  presented  in 
Chapter  VII  where  the  experimental  procedures  are  also  de¬ 
scribed. 

The  approximations  for  the  aperture  fields  may  be  used 
to  calculate  the  field  anywhere  in  the  three  regions  by  using 
the  appropriate  form  of  Eq. (63) ,  (64)  or  (65)  with  the  cor¬ 
responding  expressions  for  the  constants  in  these  equations. 
In  Chapter  VIII,  calculated  and  measured  internal  field  com¬ 
ponents  at  the  inner  surfaces  of  the  end  plates  of  a  two-port 
coaxial  cavity  are  presented.  Comparisons  between  calculated 
and  measured  impedance  and  field  data  make  possible  quanti¬ 
tative  assessments  of  the  accuracy  obtained  with  integrals 
of  the  aperture  field  solutions.  A  more  direct  validation 
of  the  accuracy  of  the  aperture  field  solutions  is  not  with¬ 
in  the  scope  of  this  effort.  Calculated  values  of  the  impe¬ 
dance  and  distributions  of  field  components  do  provide  a  de¬ 
finitive  estimate  of  the  usefulness  of  the  integral  equation 
solutions  since  these  quantities  are  of  particular  interest 
when  characterizing  a  waveguide  section  or  cavity. 
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TABLE  10.  CALCULATED  ADMITTANCE  OF  A  TWO-PORT  COAXIAL  CAVITY 
(NORMALIZED  BY  THE  INPUT  WAVEGUIDE) 


e/\ 

a  *  e/6  b  ■  c/2 

waveguide  Handbook 

d  »  3c 

Integral  Equation 

0.0010 

0.99944272  -  jO. 01812500 

1.00196814 

- 

jO. 01813840 

0.0100 

0.94773463  -  jO. 16771975 

0.95022001 

- 

jO. 16784703 

0.0500 

0.48852657  -  jO. 21444090 

0.49048806 

- 

jO. 21455535 

0.1000 

0.39078594  +  jO. 15363671 

0.39245002 

+ 

jO. 15414453 

0.2000 

0.70760169  -  jO. 21443225 

0.71040425 

- 

jO. 21394403 

0.3000 

0.53118295  +  jO. 48006425 

0.53239471 

jO.  48203314 

0.4000 

0.41303114  +  jO.  04664547 

0.41447072 

jO. 04752152 

0.5000 

1.00000000  +  jO.  95965646 

1.00239613 

jO. 95982762 

0.S6OO 

0.30982384  ♦  jO. 47922307 

0.40973071 

+ 

jO. 37257491 

0.5817 

Cutoff  Freq 

0.62882691 

+ 

jl. 31400852 

0.6000 

Above  Cutoff 

0.28442809 

♦ 

jl. 19896375 

0.7000 

Above  Cutoff 

1.22676713 

+• 

jO. 26613012 
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a  ■  c/2  b  *  2c/3 

waveguide  Handbook 

d  *  3c 

Integral  Equation 

0.0010 

0.99833581  -  jO. 03663773 

0.99847306 

- 

jO. 03663551 

0.0100 

0.85851552  -  jO. 30720201 

0.85865845 

- 

jO. 30718608 

0.0500 

0.24412759  -  jO. 21118345 

0.24422795 

- 

jO. 21110388 

0.1000 

0.18450972  +  jO. 14622063 

0.18457563 

jO. 14642800 

0.2000 

0.40485674  -  jO. 28635255 

0.40510632 

- 

jO. 28601009 

0.4000 

0.13984262  +  jO. 045510-6 

0.18990183 

jO. 04635422 

0.5000 

1.00000000  jO. 50189935 

1.00013674 

♦ 

jO. 50467629 

0.6000 

0.15744697  +  jO. 38570238 

0.15731578 

♦ 

jO. 38731123 

0.8000 

0.21939101  ♦  jO. 88059348 

0.21886714 

+ 

jO. 88348416 

0.9000 

0.12703252  +  jO.  61676062 

0.12713311 

♦ 

jO. 61875045 

0.9945 

Cutoff  Freq 

0.60087728 

■4* 

j5. 25158940 

1.1000 

Above  Cutoff 

1.21393128 

- 

jO. 29987103 
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CHAPTER  V 


INTEGRAL  EQUATIONS  FOR  PARTITIONED  CAVITIES 

A.  Introduction 

In  the  preceding  chapter,  it  is  demonstrated  that  the 
same  general  form  applies  for  both  of  the  coupled  integral 
equations  for  the  aperture  electric  fields  which  are  present 
when  two  planar  discontinuities  are  placed  in  a  coaxial 
waveguide.  This  general  form  can  be  used  to  formulate  the 
appropriate  coupled  integral  equations  for  any  combination 
of  two  of  the  discontinuities  depicted  in  Figures  3  and  4 . 

In  this  chapter,  it  is  shown  that  with  one  additional 
general  form  of  an  integral  equation,  it  is  possible  to 
treat  any  combination  of  planar  discontinuities  with  no 
limit  on  the  number.  The  use  of  such  a  procedure  is  demon¬ 
strated  by  placing  annular  disks  within  the  geometries  de¬ 
picted  in  Figures  5  and  14  such  that  the  cavities  are  further 
partitioned  into  two  regions  which  are  coupled  through  the 
aperture  in  the  partition. 

B.  Field  Components  for  a  Partitioned  One-Port  Cavity 

The  one-port,  partitioned  cavity  depicted  in  Figure  17 
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is  formed  by  placing  a  shorting  plate  across  Region  3  of  the 

two-port  cavity  of  Figure  14 .  An  annular  disk  is  attached 

to  the  outer  conductor  at  the  junction  of  Regions  2  and  3  to 

make  the  problem  more  general.  The  field  components  in  the 

three  regions  are  described  by  Egs. (63) - (65)  with  the 

accompanying  definitions  in  Eqs. (66)- (69) .  The  unknown 

aperture  electric  fields  in  each  of  the  two  discontinuities 
1  2 

are  designated  E  and  E  with  the  same  definitions  as  given 

by  Eqs. (70a&b) .  The  expressions  for  the  unknown  constants 

in  the  field  expansion  for  Region  1  are  given  by  Eqs. (75a&b) . 

The  expressions  in  Eqs. (76a-d)  are  applicable  to  Region  2  if 

2  2 

the  limits  of  integration  for  IQ  and  I2n  are  changed  from 
(a,d)  to  (a,g)  to  reflect  the  presence  of  the  annular  disk. 
The  constants  for  the  field  components  in  Region  3  are  de¬ 
rived  in  the  same  manner  as  for  the  interior  region  of  the 
one-port  cavity  in  Chapter  II.  The  constants  have  the  same 
form  as  those  defined  in  Eqs. (47)  and  (48)  and  are  given  by 
the  following  expressions: 


j2 

j2sink3(u-t)  ln*^  0 

±e±^a3nu  z2 
j2sina3n(u‘t)  M3n  ^ 


(80a) 


(80b) 
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C.  Integral  Equations  for  the  Aperture  Fields 

Two  coupled  integral  equations  for  the  transverse  elec¬ 
tric  field  components  in  the  two  apertures  of  the  geometry 
depicted  in  Figure  17  are  obtained  by  constraining  the  trans¬ 
verse  component  of  magnetic  field  to  be  continuous  across 
both  apertures.  The  first  integral  equation,  which  is  ob¬ 
tained  for  the  aperture  connecting  Regions  1  and  2,  is 

given  by  Eq.(78)  with  the  previously  mentioned  changes  to 

,  2  2 
the  limits  of  the  integrals  IQ  and  I2n*  This  is  tiie  appro- 

priate  form  of  the  integral  equation  for  an  aperture  which 

faces  another  discontinuity  in  one  axial  direction  and  an 

infinite  waveguide  in  the  other  direction.  The  aperture 

connecting  Regions  2  and  3  faces  discontinuities  in  both 

axial  directions,  thus  a  slightly  different  general  form  is 

required  for  the  integral  equation.  This  equation  is  ob- 

2  3 

tained  by  requiring  H^(p,t)  =  H^(p,t). 
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As  can  be  seen,  this  equation  contains  no  unknown  constants. 
Collecting  terms  yields 
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Because  the  discontinuity  at  z  =  u  has  no  aperture,  Eq. (81) 
does  not  represent  the  most  general  form  of  the  integral 
equation  for  an  aperture  which  faces  discontinuities  in  both 
axial  directions.  The  most  general  form  will  be  derived  in 
the  next  section  where  a  partitioned,  two-port  cavity  is 
treated . 

If  the  constants  q  and  B^n  are  known  quantities, 

Eqs. (78)  and  (81)  represent  two  integral  equations  for  the 
two  unknown  aperture  field  distributions.  Solutions  for  the 
aperture  fields  may  be  obtained  by  application  of  the  numer- 
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ical  procedure  described  in  Chapter  III.  Aperture  electric 
field  distributions  which  have  been  calculated  for  two 
different  size  partitioned,  one-port  cavities  at  two  fre¬ 
quencies  each  are  presented  in  Figures  18  and  19.  The  two 
radii  c  and  d  in  Figure  17  are  equal  for  each  of  the  cavities 
but  the  apertures  are  of  different  size.  One  of  the  fre¬ 
quencies  is  above  and  the  other  is  below  the  cutoff  frequency 
or  the  first  higher-order  mode  in  the  larger  section  of 
waveguide.  The  aperture  distributions  are  very  similar  in 
form  to  those  shown  in  Figure  8  for  the  annular  disk  and  in 
Figure  12  for  the  simple  one-port  cavity.  The  singular  be¬ 
havior  of  the  distributions  differs  in  the  two  apertures  be¬ 
cause  of  the  different  geometry  of  the  discontinuities.  No 
zero  crossing  of  the  real  or  imaginary  parts  is  observed, 
as  was  the  case  in  Figure  12  for  the  simple  one-port  cavity. 

The  input  admittance  of  the  cavity  can  be  computed  from 
the  approximate  solution  for  the  aperture  electric  field  at 
z  =  s.  Calculated  values  of  input  admittance  are  presented 
in  Table  11  for  two  different  size  cavities  as  a  function  of 
frequency.  As  in  previous  examples,  the  admittance  values 
are  computed  both  from  solutions  to  integral  equations  and 
from  transmission  line  theory  with  susceptance  corrections 
from  the  Waveguide  Handbook.  The  admittance  data  appear 
to  support  earlier  observations  that  the  transmission  line 
theory  with  corrections  is  more  accurate  than  the  integral 
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equation  approach  for  geometries  where  the  aperture  is  a 
significant  percentage  of  the  radius  of  the  larger  section 
of  waveguide.  The  integral  equation  solutions  appear  to 
give  better  results  than  the  transmission  line  theory  when 
the  aperture  is  a  comparatively  small  part  of  the  waveguide 
radius.  There  is  a  pronounced  peak  in  the  input  admittance 
at  frequencies  where  there  is  an  axial  resonance.  The  cal¬ 
culated  admittance  at  these  frequencies  has  a  real  part  that 
is  much  greater  than  zero  but  is  still  many  orders  of  magni¬ 
tude  smaller  than  the  imaginary  part.  These  calculations 
actually  give  only  a  qualitative  description  of  the  admit¬ 
tance  at  axial  resonance  because  the  value  actually  approaches 
infinity.  The  numerical  error  introduced  by  the  digital  com¬ 
puter  allows  a  calculation  to  be  made. 

D.  Field  Components  for  a  Partitioned  Two-Port  Cavity 

A  two-port  partitioned  cavity  may  be  created  from  the 
geometry  depicted  in  Figure  17  by  removing  the  shorting  plate 
at  z  *  u  and  connecting  Region  3  to  an  infinitely- long  section 
of  waveguide  of  a  different  size.  Such  a  cavity  is  illustra¬ 
ted  in  Figure  20.  The  field  components  in  the  four  regions 
are  described  by  Eqs . ( 63) - (65)  with  i  now  equal  to  1,  2,  3  and 
4.  The  unknown  aperture  electric  field  at  z  =  u  is  defined 

Ep(p,u)  =  Ep(p,u)  *  E^  .  (82) 
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TABLE  11.  CALCULATED  ADMITTANCE  OF  A  PARTITIONED,  ONE-PORT 
COAXIAL  CAVITY  (NORMALIZED  BY  THE  INPUT  WABEGUIDE) 


O.OOIO 

0.0  +  j  0.04542289 

0.00000014 

3 

0.04541979 

0.0100 

0.0  +  j  0.46065099 

0.00001455 

+ 

j 

0.46062760 

0.0500 

0.0  +  j  3.61464770 

0.00089707 

+ 

j 

3.61712020 

0.1000 

0.0  -  j  5.20242693 

0.00184342 

- 

j 

5.18534685 

0.2000 

0.0  ♦  j  2.49915937 

0.00043327 

+ 

3 

2.51429549 

0.4000 

0.0  -  jll. 66961875 

0.00889162 

- 

311.39779337 

0.5000 

0.0  +  j  0.62501141 

0.00002848 

3 

0.64550426 

0.6000 

0.0  -  j  2.00806112 

0.00027188 

- 

j 

1.99606828 

0.3000 

0.0  -  j  0.41045183 

0.00001040 

- 

3 

0.39189766 

0.9000 

0.0  -  j  1.21202495 

0.00009623 

j 

1.19927743 

0.9945 

Cutoff  Fr«q 

0.00000000 

+ 

3 

0.00008060 

1.1000 

Abova  Cutoff 

0.00022378 

- 

3 

1.81561836 

■  c/2 

b  «  2c/3  d  *  5c/6 

x  -  2c 

X  »  c 

C/\ 
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Intagral  Equation 

0.0010 

0.0  +  j  0.04542285 

0.00000014 

+ 

3 

0.04541974 

0.0100 

0.0  ♦  j  0.46059965 

0.00001454 

3 

0.46057495 

0.0500 

0.0  +  j  3.59553783 

0.00088736 

+ 

3 

3.59749163 

0.1000 

0.0  -  j  5.37370058 

0.00196982 

- 

3 

5.36018579 

0.2000 

0.0  *  j  2.52256821 

0.00044187 

+ 

3 

2.53908775 

0.4000 

0.0  -  jl2.30740830 

0.00992517 

- 

312.04200640 

0.5000 

0.0  *  3  0.00000000 

0.00000000 

+ 

3 

0.00000000 

0.6000 

0.0  -  3  1.95779087 

0.00025766 

- 

3 

1.94321006 

0.8000 

0.0  -  3  0.52793367 

0.00001805 

- 

3 

0.51626468 

0.9000 

0.0  -  3  1.29588364 

0.00011173 

- 

3 

1.29207960 

0.9945 

Cutoff  Fraq 

0.00000000 

- 

3 

0.00020133 

1.1000 

Abova  Cutoff 

0.00001291 

- 

3 

0.45044112 

Figure  20  Partitioned,  Two-Port,  Coaxial  Cavity 
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The  expressions  for  the  unknown  constants  in  the  field  ex¬ 
pansions  for  Regions  1  and  2  are  the  same  as  for  the  one-port, 
partitioned  cavity.  The  constants  for  Region  3  are  defined 
similarly  to  those  of  Region  2  since  the  region  is  bounded  on 
both  ends  by  apertures. 
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The  expressions  for  the  constants  in  Region  4  are  similar 
to  those  defined  in  Eq.  (74)  *. 
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E.  Integral  Equations  for  the  Aperture  Fields 

The  structure  depicted  in  Figure  20  has  three  apertures 
with  unknown  fields,  and,  by  following  previously  established 
procedures,  one  can  derive  three  coupled  integral  equations 
for  these  unknown  fields.  The  equation  for  the  aperture  at 
z  *  s  is  again  given  by  Eq.(78).  It  is  apparent  at  this 
point  that  the  form  of  Eq.  (78)  is  not  influenced  by  any 
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geometrical  features  beyond  the  discontinuity  at  z  =  t,  al¬ 
though  the  final  solutions  for  the  aperture  fields  certainly 
are.  The  aperture  at  z  =  t  faces  discontinuities  containing 
apertures  in  both  axial  directions.  The  integral  equation 
appropriate  for  this  aperture  is  the  most  general  of  the 
equations  to  be  developed  for  treatment  of  such  apertures 
in  discontinuities.  The  equation  is  obtained  by  requiring 
that  H^(p,t)  =  H^(p,t) : 
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The  fact  that  weighted  integrals  of  all  three  aperture 
fields  appear  in  Eq. (85)  clearly  demonstrates  the  interdepen¬ 
dence  of  these  fields.  If  more  apertures  were  added,  the 
interdependence  would  still  be  present  although  integrals  of 
all  aperture  fields  would  not  appear  in  the  same  equation. 

The  integral  equation  obtained  by  constraining  continuity  of 
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the  transverse  magnetic  field  across  the  aperture  at  z  =  u 
is  of  the  same  form  as  Eq. (79)  and  is  given  by 
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With  the  assumption  that  the  constants  Ql  B^n,  B~  Q 
and  B4n  are  known  quantities,  Eq.  (78) ,  (85)  and  (86)  form  the 
set  of  equations  necessary  to  obtain  solutions  for  the  three 
unknown  aperture  fields.  The  numerical  technique  of  Chapter 
III  is  again  suitable  for  obtaining  solutions.  There  is  no 
additional  difficulty  involved  with  treating  three  simultane¬ 
ous  equations  beyond  the  "bookkeeping"  problems.  Simple 
logic  advises  that  the  optimum  solutions  for  a  given  total 
number  of  expansion  pulses  would  be  obtained  by  dividing 
each  of  the  apertures  into  pulses  that  are  of  approximately 
equal  width.  This  effect  actually  is  diminished  significantly 
because  the  equations  involve  integrals  of  the  aperture 


fields.  It  has  been  demonstrated  in  earlier  tables  of  admit¬ 
tance  values  that  the  integrals  of  the  aperture  fields  are 
not  very  sensitive  to  the  number  of  expansion  pulses  (above 
a  usually  small  minimum  number) .  The  result  of  this  situa¬ 
tion  is  that  very  good  solutions  can  be  obtained  with  pulse 
distributions  that  vary  widely  for  different  apertures  of 
the  same  cavity  configuration.  If  there  were  reason  to  con¬ 
centrate  attention  on  one  cavity  region,  the  non-adjoining 
aperture  could  be  adequately  treated  with  comparatively  few 
expansion  pulses. 

Aperture  electric  field  distributions  which  have  been 
calculated  for  two  different  size,  partitioned,  two-port 
cavities  are  presented  in  Figures  21  and  22.  These  cavities 
have  been  created  by  adding  annular  disks  to  the  outer  con¬ 
ductors  of  the  two-port  cavities  of  Figures  15  and  16  so 
that  the  effects  of  partitioning  such  cavities  can  be  ob¬ 
served.  The  resulting  aperture  distributions  are  very  fre¬ 
quency  dependent  as  the  first  transverse  resonance  fre¬ 
quency  is  approached  and  for  frequencies  beyond.  A  quali¬ 
tative  idea  of  the  aperture  field  variation  with  frequency 
can  be  interpreted  from  the  comparisons  of  measured  and  cal¬ 
culated  field  components  in  Chapter  VIII  which  are  given  for 
a  number  of  frequencies. 

The  input  admittance  at  either  port  can  be  computed  from 
the  approximation  to  the  aperture  electric  field  at  that  port. 
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In  Table  12  comparisons  are  presented  for  admittance  values 
calculated  from  solutions  to  integral  equations  and  computed 
from  transmission  line  theory  with  susceptance  corrections 
from  the  Waveguide  Handbook.  This  data  again  appears  to 
demonstrate  that  the  integral  equation  approach  is  superior 
to  the  transmission  line  approach  when  the  aperture  is  small 
enough  that  the  singular  behavior  of  the  aperture  field 
caused  by  the  discontinuity  is  seen  over  most  of  the  aperture. 


TABLE  12.  CALCULATED  ADMITTANCE  OF  A  PARTITIONED,  TWO-PART 
COAXIAL  CAVITY  (NORMALIZED  BY  THE  INPUT  WAVEGUIDE) 


a  »  c/2 

b  »  2c/3 

d  -  Sc/6 

x  -  3c/2 

y  -  3c/2 

c/\ 
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Integral  Equation 

O.OOIO 

0.99834364  - 

jO. 03646526 

0.99889305 

- 

jO. 03645737 

0.0100 

0.35913835  - 

jO. 30578823 

0.85970885 

- 

jO. 30573056 

0.0500 

0.24708957  - 

jO. 20829428 

0.24755160 

- 

30.20801032 

0.1000 

0.19468650  + 

jO. 15497574 

0.19523825 

♦ 

jO. 15561873 

0.2000 

0.33681097  - 

jO. 26536484 

0.33590023 

- 

jO. 26378041 

0.4000 

0.21549346  + 

jO. 06722641 

0.21648982 

+ 

jO. 07016078 

0.5000 

0.96235858  - 

jO. 06022518 

0.95726980 

* 

jO. 06694721 

0.6000 

0.12576944  + 

jO. 43536054 

0.12489168 

+ 

jO. 44122015 

0.8000 

0.58503247  + 

j 1.01100177 

0.59984037 

4* 

jl. 02307863 

0.9000 

0.06010651  + 

jO. 77547372 

0.05651798 

4* 

jO. 79681264 

0.9945 

Cutoff 

Freq 

0.52758674 

4- 

j6. 67688686 

1.1000 

Above  Cutoff 

0.94391508 

4* 

jO. 03856458 
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- 
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0.0100 

0.85932100  - 
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0.85989129 

- 

30.30587505 

0.0500 

0.24701063  - 

jO. 21024318 

0.24745830 

- 
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0.1000 

0.19266498  + 

jO. 14812225 

0.19314166 

4* 

30.14879100 

0.2000 

0.35572846  - 

jO. 23677548 

0.35517085 

- 

30.23466776 

0.4000 

0.15438273  + 

jO. 06244383 

0.15359618 

+ 

30.06533893 

0.5000 

1.00000000 

jO. 59496446 

1.00054745 

+ 

jO. 60669547 

0.6000 

0.20108709  + 

jO. 39082365 

0.20211057 

■f 

30.39476634 

0.8000 

0.46097070  + 

jl. 23927333 

0.47402267 

+ 

31.24904843 

0.9000 

0.06190716  + 

jO. 71737474 

0.06739482 

+ 
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0.9945 

Cutoff 

?req 

0.60245035 

35.48157554 

1.1000 

Above  Cutoff 

0.23160684 

♦ 

30.48464206 
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CHAPTER  VI 


INTEGRAL  EQUATION  FOR  A  TRUNCATED  COAXIAL  INNER  CONDUCTOR 

A.  Introduction 

The  discontinuities  which  are  the  subjects  of  the  anay- 
ses  in  the  preceding  chapters  are  all  coaxial  in  nature.  It 
is  demonstrated  that,  in  theory,  any  combination  of  planar, 
coaxial  discontinuities  can  be  treated  by  proper  application 
of  one  or  more  of  the  general  equation  forms  developed  in 
Chapters  II,  IV  and  V.  An  additional  type  of  discontinuity 
is  formed  by  truncating  the  center  conductor  of  a  coaxial 
waveguide  and  it  can  be  treated  by  a  similar  approach  to 
that  used  for  the  coaxial  discontinuities.  The  analysis 
for  such  a  structure  is  presented  in  this  chapter,  and  it 
is  subsequently  demonstrated  that  the  expressions  which  are 
derived  for  the  circular  waveguide  region  are  analogous  to 
those  derived  for  the  coaxial  regions,  noting  that  there 
are  no  corresponding  TEM  terms.  This  fact  facilitates  the 
derivation  of  integral  equations  for  waveguides  where  the 
previously-treated  discontinuities  and  the  truncated  center 
conductor  are  mixed. 


B.  Expressions  for  Field  Components 

A  coaxial  waveguide  with  a  truncated  inner  conductor  is 
depicted  in  Figure  23.  The  outer  radius  of  the  waveguide  is 
also  reduced  in  the  plane  of  the  truncation  to  demonstrate 
that  other  discontinuities  can  be  treated  simultaneously. 

For  Region  1  (z  <  s)  ,  the  field  component  definitions  of 
Eqs. (17) - (22)  are  appropriate.  For  the  circular  waveguide 
region,  (z  >  s) ,  the  field  components  must  be  derived  from 
an  appropriate  magnetic  vector  potential  which  is  a  solution 
to  Eq. (1) .  The  axial  component  of  magnetic  vector  potential 
is  sufficient  to  provide  a  complete  field  definition,  as  is 
demonstrated  for  the  coaxial  case  in  Appendix  A.  The  mag¬ 
netic  vector  potential  component  is  defined  as  the  product 
(no  variation  with  4> )  : 

A°(p,z)  =  RC(P)  Z(z)  (87) 

where  the  axial  function,  Z(z),  which  satisfies  Eq. (2) ,  has 
the  form 

Z (z )  =  cx  e’jkzz  +  c 2  ejkzz  (88) 

and  the  radial  function,  Rc(p),  is  a  solution  of  Eq. (3)  and 
has  the  form 

Rc(p)  =  c3  JQ(kpp)  +  c4  N0(kpp)  (89) 

The  boundary  condition  that  the  field  have  a  finite 
value  at  p  *  0  eliminates  NQ(kpp)  as  a  possible  solution. 


Figure  23.  Discontinuity  Formed  in  Coaxial  Waveguide  by  Truncating  Center  Conductor 


A  second  boundary  condition  requires  that  the  tangential 


component  (in  this  case,  E  )  of  the  electric  field  be  zero 

z 

at  the  inner  surface  of  the  waveguide  (at  p  =  c) .  The  re¬ 
sulting  function  is 


R  (p)  =  c 3  J0UnP) 


(90) 


where  the  are  the  infinite  set  of  kp  which  satisfy 


Jq(X^c)  0  ,  n  ~  1,  2 ,  •• 


(91) 


The  eigenvalue  representation  for  kz  is  now  given  by 

kz  "  8n  =  k2  -  xn  * 


The  following  definitions  for  Bn  are  used 


8n  = 


42  _  ,2  J2  2 

k2  An  '  k2  An 
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6n  =  ' 


1  'i2  -1? 

3  An  *2 


2  2 
n  >  k2 


(92b) 


The  resulting  expression  for  the  magnetic  vector  potential  is 

00 

Az(P,z)  »  l  (c+  e"36nz  +  cn  e3Snz)  JQ(XnP)  •  (93) 

n=l 


The  field  components  for  Region  2  can  now  be  computed  from 
Eqs.(ll),  (12),  and  (13)  operating  on  Eq. (93) .  The  resulting 
expressions  are 
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The  constants  c  and  c  are  now  redefined  to  give  the  field 
n  n  3 

expressions  a  more  standard  form: 
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and 
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The  field  components  are  now  given  by 
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C.  Definitions  of  Constants  in  Field  Expressions 

The  procedure  to  obtain  useful  definitions  for  the  un¬ 
known  constants  of  Eqs. (95)- (97)  employs  similar  orthogo- 
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nality  properties  to  those  of  the  method  used  in  Chapter  II . 
The  first  step  is  to  multiply  both  sides  of  Eq.(95)  by 
P  and  then  integrate  both  sides  over  the  interval 

(0,c).  This  yields 


fc  2 

J  A  • 


(X  o)  dp  = 


#c  00 

Jo  P  (Cn  e’*lSnS  +  cn  e:,6nS)  Ji^mp)  Jl(Xnp)  dp  ' 


Because  of  the  orthogonality  properties  of  the  integral  on 
the  right-hand  side,  Eq.(98)  has  zero  value  except  when  m  =  n. 
The  constants  may  now  be  defined  by 

<  e"j6nS  +  C"  ej6nS  =  - - - i -  Pe®  p  J.  (X  p)  dp 

n  n  re  -Jain 

~X  I  p[J. (X  p) ] ap  a 

nJ  0  A  n  fQQ) 


where  the  radial  component  of  the  electric  field  in  the 


aperture  is  defined  by 


EpCp) 


E  ,  a  <  p  <  c 

cl 

0  ,  0  <  p  <  a 


(100) 


The  normalization  integral  in  Eq. (99)  has  the  value: 
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D.  Enforce  Continuity  of  Transverse  Fields  in  the  Aperture 
An  equation  containing  the  integral  of  the  unknown  aper¬ 
ture  electric  field  is  obtained  by  constraining  the  trans¬ 
verse  component  of  magnetic  field  to  be  continuous  across 

1  2 

the  aperture,  i.e.,  H^fP/S)  =  H^(p,s).  This  yields 
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An  integral  equation  appropriate  for  the  geometry  of  Figure 
23  is  easily  obtained  from  Eq. (102)  if  the  waves  incident  on 
the  aperture  from  both  directions  are  fully  defined.  As  an 
example,  assume  that  the  excitation  traveling  in  the  +z 
direction  in  Region  1  consists  of  a  TEM  wave  with  amplitude 
Bq.  The  existence  of  only  a  TEM  mode  implies  that  all  the 
B+  =  0.  It  is  also  assumed  that  the  waveguide  on  the  right- 
hand  side  of  the  aperture  is  either  infinitely  long  or  termina¬ 
ted  in  a  matched  load,  either  of  which  implies  that  all  the 
C~  *  0.  The  resulting  integral  equation  for  the  coaxial  wave¬ 
guide  with  a  truncated  inner  conductor  is,  from  Eq.(99), 
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E.  Application  to  a  One-Port  Cavity 

The  geometry  depicted  in  Figure  24  possesses  two  discon¬ 
tinuities:  the  first,  at  z  =  s,  is  the  now-familiar  step- 
discontinuity  in  the  outer  conductor  of  a  coaxial  waveguide: 
and  the  second,  at  z  =  t,  is  the  truncated  center  conductor 
of  a  coaxial  line.  The  two  unknown  aperture  field  distribu¬ 
tions  are  obtained  from  solutions  to  a  pair  of  coupled  inte¬ 
gral  equations  which  are  derived  according  to  previously  de¬ 
termined  procedures.  It  is  convenient  to  define  the  field 
components  in  Regions  1  and  2  with  the  expressions  in  Eqs.(63)- 
(65)  with  i  =  1  and  2.  The  constants  in  these  expressions  are 
defined  in  Eqs. (71) - (73)  where  the  notation  for  the  weighted 
integrals  of  the  aperture  fields  retains  the  same  meaning. 

The  field  components  for  Region  3  are  given  by  Eqs (95) -(97) 
where  the  subscripts  and  superscripts  which  denote  parameters 
corresponding  to  a  particular  region  are  changed  from  2  to  3. 
The  constants  in  Eqs. ( 95) - ( 97)  are  evaluated  by  application 
of  the  boundary  condition  requiring  zero  tangential  electric 
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Figure  24.  One-Port,  Coaxial  Cavity  with  the  Coaxial  Center  Conductor  Truncated 


field  on  the  plate  at  z  =  u.  These  constants,  which  are 
analogous  in  form  to  those  defined  in  Eq.(48),  are  given  by 


te^n® 


n  -jX  d2  sin8„(u-t)  J?(X  d) 
n  n  in 


rd  2 

I  p  Jl(Xnp)  dp  (104) 


The  first  of  the  required  pair  of  integral  equations  is  given 
by  Eq. (78)  with  z  =  s.  The  second  integral  equation,  which 
results  from  the  substitution  of  the  proper  expressions  for 
the  constants  into  Eq. (102)  with  z  =  t,  is  quite  similar  to 
Eq. (81) : 
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where  the  X  are  the  solutions  to 
n 


Jq ( X ^d )  =  0  ,  n  =  1,  2,  ... 

and  the  Bn  are  defined  by  Eq. (92) . 

It  is  a  straightforward  procedure  to  combine  the  dis¬ 
continuity  depicted  in  Figure  23  with  any  or  all  of  the 
coaxial  discontinuities  illustrated  in  Figures  3  and  4  to 
create  an  arbitrarily  complex  cavity  structure.  The  equa- 


tion  for  the  truncated  center  conductor  would  be  included 


where  necessary  to  obtain  a  tractable  set  of  integral  equa¬ 
tions.  A  similar  treatment  of  the  junction  of  two  circular 
waveguides  of  differing  radii  can  be  accomplished  with  a 
simple  extension  of  the  material  presented  in  this  chapter. 
Such  a  configuration  was  not  of  interest  during  this  inves¬ 
tigation  primarily  because  circular  waveguides  do  not  sup¬ 
port  the  TEM  mode. 


F.  Numerical  Solution  Considerations 

The  numerical  procedure  employed  to  obtain  solutions  to 
Eqs. (103)  and  (104)  is,  in  principle,  the  same  technique  de¬ 
scribed  in  Chapter  III.  The  expansion  for  the  aperture  elec 
trie  field  defined  by  Eqs. (55)  and  (56)  is  appropriate.  The 
resulting  integral  is  easily  evaluated  to  yield 
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The  testing  function  set  given  by  Eq. (61)  is  used  to 


generate  the  N  algebraic  equations  necessary  to  obtain  a 
solution  for  the  aperture  field.  The  two  sets  of  algebraic 
equations  obtained  from  Eqs.(78)  and  (105)  are  sufficient 
to  obtain  solutions  for  both  aperture  electric  fields  when 
the  constants  of  Eq. (78)  are  defined.  In  Figures  25  and  26 
are  presented  calculated  aperture  field  distributions  for  the 
coaxial-to-circular  waveguide  junction  of  Figure  23  without 
the  step  in  the  outer  conductor  and  the  one-port  cavity  of 
Figure  24. 

G.  Susceptance  of  a  Truncated  Center  Conductor 

The  lumped  TEM  susceptance  introduced  into  a  transmission 
line  when  the  coaxial  center  conductor  is  truncated  can  be 
calculated  by  following  the  procedure  described  in  Chapter  II. 
Once  the  aperture  distribution  is  obtained,  its  integral  can 
be  used  to  evaluate  the  voltage  reflection  coefficient  through 
Eq.  (54)  and  the  susceptance  can  then  be  computed  with  the  aid 
of  Eq. (51)  where  the  characteristic  admittance  (in  the  TEM 
sense  of  the  circular  waveguide  section,  Yg,  is  equal  to  zero. 
The  resulting  expression  is 

il  .  i  -  r<«)  .  no?, 

y0  i  +  ns) 

The  expression  on  the  right-hand  side  of  Eq. (107)  is  the 
definition  of  the  normalized  admittance  at  z  *  s  expressed  in 
terms  of  the  voltage  reflection  coefficient.  For  frequencies 
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below  the  cutoff  frequency  for  the  first  propagating  mode  in 
the  circular  waveguide,  the  TEM  admittance  is  imaginary.  At 
frequencies  above  cutoff  for  one  or  more  circular  waveguide 
modes,  a  portion  of  the  incident  energy  is  converted  to  a 
propagating  mode(s)  in  the  circular  waveguide.  Under  the 
latter  circumstance,  the  admittance  is  complex.  Whereas 
previous  susceptance  expressions  are  generally  applicable  to 
TEM  waves  incident  from  either  direction,  the  expression  in 
Eq.  (107)  is  only  valid  for  a  TEM  wave  incident  from  the 
coaxial  waveguide  side  of  a  truncated  center  conductor  dis¬ 
continuity.  In  Table  13  are  presented  TEM  admittance  values 
calculated  from  solutions  to  integral  equations  for  four 
different  ratios  of  outer  coaxial  radius  to  inner  radius. 
There  is  no  corresponding  expression  in  the  Waveguide  Hand¬ 
book  for  comparison,  although  Marcuvitz  does  provide  a  sus¬ 
ceptance  expression  for  a  coaxial-to-circular  waveguide 
junction  where  the  inner  conductor  is  also  a  circular  wave¬ 
guide.  The  b/X  =  0.3827  frequency  for  each  geometry  is  the 
cutoff  frequency  for  the  TM^  mode  in  the  circular  waveguide. 
As  for  the  other  planar  discontinuities,  the  admittance  has  a 
non-zero  real  part  when  power  is  converted  to  propagating 
higher-order  modes. 

As  mentioned  in  Chapter  I,  the  coaxial-to-circular  wave¬ 
guide  junction  has  been  treated  by  a  number  of  other  investi¬ 
gators.  In  Table  14  comparisons  of  calculated  capacitance 
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are  presented  for  six  junctions  of  varying  size  along  with 
measured  values  for  three  of  the  junctions.  The  integral 
equation  results  are  computed  with  20  pulses  across  the 
aperture  and  30  series  terms  and  with  50  pulses  and  50 
series  terms.  The  two  calculations  give  a  qualitative  feel 
for  the  convergence  of  the  integral  equation  approach.  The 
integral  equation  results  compare  extremely  well  with  the 
other  values  in  the  table  and  indicate  that  this  approach 
would  be  very  useful  when  extreme  accuracy  is  desired. 


TABLE  13.  CALCULATED  ADMITTANCE  OF  COAXIAL-TO-CIRCULAR  WAVEGUIDE  JUNCTIONS 
(NORMALIZED  BY  THE  COAXIAL  WAVEGUIDE) 


b/X 

a  •  10  a 

b  *  60  a 

b/X 

4  »  10  si 

0.0100 

0.00191089  + 

j  0 . 017443  52 

0.0050 

0.00040563  * 

0.0500 

0.00191390  + 

jO. 08771871 

0.0250 

0.00040576  - 

0.1000 

0.00192371  * 

j  0 . 1787  08  52 

0.0500 

0.00040616  * 

0.2000 

0.00197098  ♦ 

jO. 38887531 

0.1000 

0.00040787  - 

0.3000 

0.00210462  ♦ 

jO. 71558373 

0.2000 

0.00041627  * 

0.3500 

0.00231081  ♦ 

jl. 06790249 

0.2500 

0.00042511  * 

0.3600 

0.00240087  + 

jl. 20129638 

0.3000 

0.00044171  * 

0.3300 

0.00291225  + 

jl. 86998604 

0.3500 

0.00048546  - 

0.3827 

0.00348897  + 

j2. 55084881 

0.3750 

0.00057774  - 

0.4000 

1.31789521  * 

3 1 .45373287 

0.3800 

0.00065154 

0.5000 

1.09293073  * 

jO. 44577374 

0.3827 

0.00090171  " 

0.6000 

0.77851531  * 

jO. 31616103 

0.4000 

2.47353062 

b/X 

a  »  0.001521 

b  -  0.0035 

b/x 

a  *  3  si 

0.0053 

0.00019634  + 
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0.0050 

0.00514131 

0.0292 

0.00019640  + 
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0.0250 

0.00514312 

0.0583 
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0.0500 

0.00514884 

0.1167 
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30.26241962 

0.1000 

0.00517258 

0.1750 

0.00019907  + 

jO. 41843924 

0.2000 

0.00528542 

0.2333 
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jO. 61781907 
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0.00539763 

0.2917 

0 .  Q00208Q6  + 
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0.3000 

0.00559469 

0.3500 
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0.3500 
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0.3617 
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0.3750 

0.00632352 

0.3733 
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j2. 78448222 

0.3800 
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0.4000 
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b  »  30  si 

jO. 0073853C 
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TABLE  14.  COMPARISONS  OF  CALCULATED  AND  MEASURED  CAPACITANCE 


CHAPTER  VII 


ADMITTANCE  MEASUREMENTS  FOR  TWO-PORT  COAXIAL  CAVITIES 

A.  Introduction 

During  the  course  of  the  present  investigation,  two  sets 
of  coaxial  cavities  were  constructed  and  a  number  of  measure¬ 
ments  of  the  electromagnetic  properties  of  the  cavities  were 
performed.  The  term  "sets  of  cavities"  is  used  because  a 
number  of  interchangeable  parts  were  constructed  for  each 
cavity  which  allows  it  to  be  configured  as  a  simple  or 
partitioned  cavity  and,  for  one  set,  as  either  a  one-  or  two- 
port  cavity.  Descriptions  of  the  cavities  and  interesting  de¬ 
tails  of  their  construction  are  given  in  Appendix  E.  It  was 
found,  as  anticipated,  that  the  voltage  standing  wave  ratio 
(VSWR)  on  the  waveguide  feeding  a  one-port  cavity  was  far  too 
large  to  measure  accurately.  Consequently,  the  admittance 
measurements  presented  in  this  chapter  are  only  for  the  two- 
port  configuration  with  one  port  being  excited  by  an  RF 
generator  and  a  matched  load  placed  at  the  other  port. 

B.  Measurement  Procedure  and  Cavity  Configuration 

A  schematic  representation  of  the  apparatus  used  in  per¬ 
forming  the  admittance  measurements  is  given  in  Figure  27 . 
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Frequency  Counter 


Figure  27.  Schematic  Representation  of  Instrumentation  For  Admittance 
Measurements  for  Two-Port  Coaxial  Cavities 


The  primary  instrument  used  is  a  General  Radio  (now  GenRad) 

50  ohm  precision  slotted  line.  Three  quantities  were 
measured  on  the  slotted  line  which  together  provided  the 
necessary  data  for  calculating  the  admittance  of  the  cavity 
at  the  discontinuity  where  the  input  waveguide  joined  the 
cavity  proper  (z  =  s  in  Figure  14).  The  first  quantity  was 
the  VSWR  on  the  transmission  line  between  the  cavity  and 
the  source  generator.  The  second  quantity  was  the  location 
of  a  minimum  of  the  voltage  on  the  slotted  line.  The  third 
was  the  location  of  a  minimum  of  the  voltage  when  an  electri¬ 
cal  short  was  placed  at  the  location  of  the  cavity  input. 

The  procedures  for  measuring  this  data  and  the  subsequent 
calculations  of  the  admittance  are  described  very  well  in  a 
number  of  books  and  in  particular.  Ref. [25]  which  was  con¬ 
sulted  during  this  investigation.  The  results  presented  in 
the  next  section  for  each  frequency  are  the  averages  of  cal¬ 
culations  made  using  all  the  possible  combinations  of 
measured  standing  wave  minima.  At  almost  all  frequencies, 
this  total  was  nine  separate  calculations.  All  measured 
values  are  normalized  by  20  millimhos. 

The  input  admittance  was  measured  for  three  cavity  con¬ 
figurations: 

1.  Simple  cavity,  no  partition 

2.  Partitioned  cavity,  aperture  diameter  =  2.540  cm 

3.  Partitioned  cavity,  aperture  diameter  =  8.006  cm. 
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The  dimensions  of  the  various  cavity  sections  are  given  in 
Appendix  E.  During  the  measurements  only  an  air  dielectric 
was  used  in  the  cavities.  The  output  waveguide,  which  is  a 
section  of  50-ohm  air  line,  was  terminated  in  a  50-ohm  load 
at  a  point  13.472  cm  beyond  the  output  port  of  the  cavity. 

The  admittance  was  measured  at  50  MHz  intervals  over  a  band¬ 
width  which  differed  slightly  for  each  cavity.  The  fre¬ 
quency  intervals  measured  for  each  cavity  configuration  are 
given  below. 

Configuration  No.  Frequency  Range 

1  300  -  1700  MHz 

2  400  -  2000  MHz 

3  400  -  1700  MHz 

The  cutoff  frequency  for  the  first  higher-order  mode  (TMQ1) 
in  the  larger  section  of  waveguide  is  1325  MHz,  a  point  which 
is  covered  by  all  three  measurement  intervals. 

C.  Presentation  of  Results 

The  real  and  imaginary  parts  of  the  measured  admittance 
for  each  configuration  are  graphically  presented  in  Figures 
28-30.  Also  plotted  on  each  graph  are  the  real  and  imaginary 
parts  of  the  admittance  calculated  from  solutions  to  the 
appropriate  integral  equations.  In  Figures  31-33,  the  ad¬ 
mittances  computed  from  transmission  line  theory  with 
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susceptance  corrections  are  plotted  together  with  the  pre¬ 
viously  plotted  curves.  The  effects  of  the  size  of  the  ap¬ 
erture  in  the  partition  are  illustrated  in  Figure  34  where 
the  calculated  values  of  admittance  for  the  two  cases  are 
presented  together. 

The  transmission  line  theory  values  agree  quite  well 
with  the  other  results  for  frequencies  well  below  the 
cutoff  in  the  larger  cavity  section;  however,  the  results  of 
this  procedure  begin  to  display  significant  errors  as  the 
cutoff  frequency  is  approached.  The  admittances  computed 
from  the  integral  equation  solutions  display  excellent  agree¬ 
ment  with  measured  results  at  almost  all  sampled  frequencies. 
It  is  very  difficult  to  distinguish  any  differences  in  the 
two  for  frequencies  below  1900  MHz.  It  is  probable  that 
measurement  error  was  becoming  significant  at  frequencies 
above  1900  MHz.  This  contention  is  supported  by  the  results 
shown  in  Tables  15,  16  and  17  where  the  mean  and  standard  de¬ 
viation  for  the  measurements  at  each  frequency  are  presented 
along  with  the  admittance  computed  from  solutions  to  integral 
equations.  The  standard  deviations  of  the  real  and  imaginary 
parts  of  the  admittance  are  calculated  from  the  expression 


where 


y„  =  set  of  N  measurement  results 
n 

y  =  mean  of  N  measurement  results 

The  standard  deviations  for  both  the  real  and  imaginary  parts 
of  the  measured  values  grow  substantially  larger  as  the  fre¬ 
quency  increases  beyond  1900  MHz.  The  standard  deviations 
for  both  parts  of  the  admittance  are  also  large  near  peaks 
in  the  admittance  which  correspond  to  axial  resonances  of 
the  particular  cavity  (length  equal  to  odd  multiple  of 
one-half  wavelength) .  Three  of  the  four  peaks  on  the  ad¬ 
mittance  curve  for  the  simple  two-port  cavity  can  be 
identified  easily.  The  axial  length  of  the  cavity  is  equal 
to  one-half  the  free-space  wavelength  at  1000  MHz  which  is 
the  location  of  the  first  peak.  The  second  peak  occurs  at 
1325  MHz  which  is  the  first  transverse  resonance  frequency. 
The  peak  near  1650  MHz  is  also  immediately  identifiable, 
since  it  represents  the  axial  resonance  of  the  first  higher- 
order  mode.  The  fourth  peak  occurs  at  2000  MHz  where  the 
axial  length  of  the  cavity  is  equal  to  one  free-space  wave¬ 
length.  All  of  these  resonances  are  characterized  by  a 
zero  crossing  of  the  imaginary  part  of  the  admittance  with 
a  steep  negative  slope. 

The  behavior  of  the  admittance  curves  for  the  parti¬ 
tioned  cavities  becomes  more  difficult  to  interpret.  The 
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one-half  wavelength,  resonance  peak  for  the  total  axial 
length  of  the  cavity  occurs  near  600  MHz  for  both  cavities. 

The  nature  of  this  peak  is  verified  by  its  occurance  in 
Figure  34  where  the  admittance  of  an  unpartitioned  cavity 
of  the  same  length  is  presented.  Except  for  the  transverse 
resonance  peak  which  occurs  at  1325  MHz  for  all  structures, 
the  introduction  of  the  partitions  significantly  changes  the 
admittance  response  of  the  two-port  cavity  above  the  pri¬ 
mary  TEM  axial  resonance.  The  most  prominent  feature  of  the 
admittance  curves  for  the  partitioned  cavities  is  a  very 
large,  sharp  peak  between  1100  and  1200  MHz.  The  peak 
appears  to  be  directly  related  to  the  existence  of  the 
partition  and  its  aperture.  The  value  of  the  peak  and  its 
position  vary  as  the  size  of  the  aperture  is  changed,  while 
the  other  features  of  the  admittance  curve  change  very  little. 
The  exact  physical  origin  of  this  resonance-like  peak  has 
not  been  explained  but  it  is  almost  certainly  an  axial 
resonance  since  it  is  also  described  by  the  transmission  line 
theory  calculations.  It  is  also  interesting  to  note  in 
Figure  35  the  almost  precise  agreement  of  the  admittance  for 
the  two  partitioned  cavities  between  the  transverse  resonance 
frequency  (1325  MHz)  and  1800  MHz.  In  this  region,  the  par¬ 
tition  has  a  definite  effect,  but  the  size  of  the  aperture 
appears  to  have  little  influence. 
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TABLE  15.  CALCULATED  AND  MEASURED  ADMITTANCE  OF  A 
TWO-PORT  COAXIAL  CAVITY 


Freq 

MHz 

Integral  Eq. 
Calculation 

Mean  of  Measured 
Results 

Standard 

Deviation 

Real 

Imag 

Real 

Imag 

Real 

Imag 

300 

0.0835 

-0.1102 

0.0589 

-0.1164 

0.0000 

0.0000 

350 

0.0696 

-0.0563 

0.0656 

-0.0630 

0.0000 

0.0000 

400 

0.0614 

-0.0076 

0.0578 

-0.0188 

0.0000 

0.0041 

450 

0.0569 

0.0379 

0.0513 

0.0315 

0.0000 

0.0029 

500 

0.0552 

0.0821 

0.0497 

0.0712 

0.0022 

0.0044 

550 

0.0560 

0.1265 

0.0546 

0.0887 

0.0041 

0.0088 

600 

0.0596 

0.1726 

0.0550 

0.1349 

0.0007 

0.0075 

650 

0.0665 

0.2222 

0.0628 

0.1917 

0.0001 

0.0037 

700 

0.0785 

0.2774 

0.0677 

0.2501 

0.0009 

0.0089 

750 

0.0986 

0.3415 

0.0816 

0.3253 

0.0009 

0.0122 

800 

0.1339 

0.4189 

0.1133 

0.3993 

0.0023 

0.0080 

850 

0.2003 

0.5150 

0.2047 

0.4961 

0.0059 

0.0120 

900 

0.3382 

0.6278 

0.3266 

0.5670 

0.0091 

0.0250 

950 

0.6370 

0.6872 

0.6197 

0.6744 

0.0076 

0.0097 

1000 

1.0325 

0.3635 

1.Q218 

0.3036 

0.0096 

0.0065 

1050 

0.8430 

-0.1719 

0.7745 

-0.2245 

0.0109 

0.0104 

1100 

0.4723 

-0.2242 

0.4431 

-0.2276 

0.0012 

0.0049 

1150 

0.2889 

-0.1084 

0.2706 

-0.1172 

0.0004 

0.0055 

1200 

0.2164 

0.0223 

0.1924 

0.0186 

0.0001 

0.0083 

1250 

0.2195 

0.1676 

0.2201 

-0.0348 

0.0138 

0.2584 

1300 

0.4694 

0.3747 

0.5048 

0.3545 

0.0029 

0.0069 

1350 

0.2015 

-0.4217 

0.1410 

-0.4259 

0.0010 

0.0097 

1400 

0.0069 

0.0073 

0.0088 

0.0275 

0.0009 

0.0061 

1450 

0.0014 

0.1750 

0.0053 

0.1892 

0.0005 

0.0068 

1500 

0.0021 

0.3093 

0.0065 

0.2467 

0.0009 

0.00062 

1550 

0.0104 

0.4704 

0.0137 

0.3974 

0.0008 

0.0070 

1600 

0.0747 

0.7450 

0.0688 

0.6659 

0.0005 

0.0072 

1650 

0.7550 

1.1941 

0.6050 

1.1892 

0.0089 

0.0121 

1700 

0.9587 

-0.0166 

0.8626 

-0.0561 

0.0021 

0.0070 
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TABLE  16.  CALCULATED  AND  MEASURED  ADMITTANCE  OF  A  PARTITIONED, 

TWO-PORT  CAVITY  -  SMALL  APERTURE 


Freq 

Integral  Eq. 

Mean  of 

Measured 

Standard 

Deviation 

MHz 

Calculation 

Results 

Real 

Xmag 

Real 

Xmag 

Real 

Xmag 

400 

0.1214 

0.2097 

0.1080 

0.2163 

0.0003 

0.0061 

4S0 

0.2177 

0.3186 

0.1937 

0.3445 

0.0047 

0.0384 

500 

0.5189 

0.3565 

0.4668 

0.3815 

0.0011 

0.0120 

550 

0.6174 

-0.1778 

0.5174 

-0.2553 

0.0026 

0.0080 

600 

0.2199 

-0.2130 

0.1608 

-0.2097 

0.0003 

0.0052 

650 

0.0937 

-0.0790 

0.0704 

-0.0833 

0.0005 

0.0053 

700 

0.0527 

0.0228 

0.0435 

0.0173 

0.0000 

0.0149 

750 

0.0360 

0.1046 

0.0287 

0.1186 

0.0012 

0.0183 

800 

0.0284 

0.1794 

0.0257 

0.1559 

0.0005 

0.0735 

850 

0.0253 

0.2565 

0.0227 

0.2346 

0.0008 

0.0066 

900 

0.0252 

0.3454 

0.0227 

0.3387 

0.0013 

0.0108 

950 

0.0287 

0.4624 

0.0283 

0.4648 

0.0007 

0.0100 

1000 

0.0395 

0.6466 

0.0422 

0.6737 

0.0021 

0.0130 

1050 

0.0793 

1.0451 

0.0908 

1.0666 

0.0080 

0.0855 

1100 

0.6530 

3.2300 

1.3069 

4.8789 

0.0729 

0.2019 

1150 

0.1574 

-1.6821 

0.1239 

-1.3361 

0.0126 

0.1010 

1200 

0.0111 

-0.4278 

0.0134 

-0.3652 

0.0007 

0.0090 

1250 

0.0026 

-0.0541 

0.0065 

-0.0148 

0.0012 

0.0061 

1300 

0.0320 

0.3348 

0.0408 

0.3958 

0.0008 

0.0066 

1350 

0.1488 

-0.5198 

0.1056 

-0.4640 

0.0010 

0.0117 

1400 

0.0104 

0.0087 

0 . 0111 

0.0060 

0.0010 

0.0059 

1450 

0.0034 

0.1760 

0.0061 

0.1735 

0.0006 

0.0067 

1500 

0.0038 

0.3101 

0.0067 

0.3042 

0.0005 

0.0064 

1550 

0.0094 

0.4669 

0.0109 

0.4602 

0.0008 

0.0076 

1600 

0.0462 

0.7258 

0.0308 

0.7138 

0.0007 

0.0080 

1650 

0.5148 

1.3673 

0.4281 

1.4438 

0.0072 

0.0169 

1700 

1.1483 

-0.1735 

0.9656 

-0.2007 

0.0076 

0.0053 

1750 

0.4774 

0.0913 

0.3951 

0.1521 

0.0008 

0.0060 

1800 

0.5003 

0.3702 

0.4458 

0.4568 

0.0052 

0.0129 

1850 

0.9779 

0.4549 

1.1128 

0.4568 

0.0122 

0.0071 

1900 

0.8825 

-0.5633 

0.4745 

-0.3327 

0.0095 

0.0271 

1950 

0.1796 

-0.3126 

0.5640 

-0.4516 

0.0592 

0.1067 

2000 

0.0407 

0.0151 

0.1912 

-0.3896 

0.0232 

0.1763 
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TABLE  17  CALCULATED  AND  MEASURED  ADMITTANCE  OF  A  PARTITIONED, 

TWO-PORT  CAVITY  -  LARGE  APERTURE 


Freq 

Integral  Eq. 

Mean  of 

Measured 

Standard 

Deviation 

MHZ 

Calculation 

Results 

Real 

Imag 

Real 

Imag 

Real 

Imag 

400 

0.0828 

0.1786 

0.0631 

0.2179 

0.0002 

0.0062 

450 

0.1194 

0.2703 

0.0759 

0.1292 

0.0007 

0.03S9 

500 

0.2178 

0.3925 

0.1719 

0.4055 

0.0006 

0.0048 

550 

0.5350 

0.4829 

0.4575 

0.4419 

0.0057 

0.0012 

600 

0.8247 

-0.1037 

0.7436 

-0.2050 

0.0013 

0.0012 

650 

0.3197 

-0.2802 

0.2075 

-0.2756 

0.0011 

0.0103 

700 

0.1348 

-0.1311 

0.1154 

-0.1219 

0.0005 

0.0165 

750 

0.0757 

-0.0148 

0.0439 

0.0259 

0.0013 

0.0225 

800 

0.0521 

0.0753 

0.0388 

0.1015 

0.0006 

0.0770 

850 

0.0415 

0.1542 

0.0358 

0.1423 

0.0007 

0.0064 

900 

0.0375 

0.2326 

0.0291 

0.2326 

0.0003 

0.0126 

950 

0.0381 

0.3199 

0.0305 

0.3121 

0.0002 

0.0091 

1000 

0.0444 

0.4307 

0.0269 

0.4367 

0.0010 

0.0097 

1050 

0.0624 

0.5986 

0.0505 

0.5757 

0.0026 

0.0542 

1100 

0.1254 

0.9411 

0.1503 

0.9552 

0.0037 

0.0126 

1150 

0.8239 

2.4444 

0.9953 

2.6990 

0.2599 

0.3148 

1200 

0.3237 

-1.5965 

0.2454 

-1.3653 

0.0065 

0.0391 

1250 

0.0133 

-0.3077 

0.0088 

-0.2577 

0.0013 

0.0074 

1300 

0.0215 

0.2384 

0.0206 

0.2754 

0.0003 

0.0076 

1350 

0.1716 

-0.5995 

0.0909 

-0.5370 

0.0019 

0.0110 

1400 

0.0083 

-0.0344 

0.0094 

-0.0370 

0.0002 

0.0049 

1450 

0.0026 

0.1455 

0 . 0  "  ~  2 

0.4004 

0.0026 

0.3960 

1500 

0.0026 

0.2815 

0.0.55 

0.2691 

0.0006 

0.0053 

1550 

0.0039 

0.4244 

0.0089 

0.4189 

0.0005 

0.0084 

1600 

0.0241 

0.6780 

0.0196 

0.6566 

0.0004 

0.0091 

1650 

0.3426 

1.3212 

0.2680 

1.3458 

0.0049 

0.0158 

1700 

1.2635 

-0.2282 

0.9497 

-0.3228 

0.0153 

0.0066 
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CHAPTER  VIII 


INTERNAL  FIELD  MEASUREMENTS  IN  COAXIAL  CAVITIES 

A.  Introduction 

In  Chapter  II  it  is  demonstrated  that  the  rotational 
symmetry  of  the  coaxial  cavities  under  investigation  limits 
the  number  of  field  components  excited  within  the  cavities. 
With  only  the  Ep,  Ez  and  Hp  components  present  inside  a 
cavity,  it  is  possible  to  construct  probes  which  do  not 
couple  to  multiple  field  components  and  to  provide  access 
to  the  cavity  interiors  without  seriously  perturbing  the 
field  within  the  cavity.  These  factors  have  made  possible 
an  extensive  program  of  interior  field  measurements  on  all 
possible  configurations  of  the  two  sets  of  coaxial  cavities 
described  in  Appendix  E.  The  field  distributions  measured 
within  the  cavities  are  compared  to  distributions  computed 
from  solutions  to  the  appropriate  integral  equations  for 
the  same  geometries.  The  measured  results  serve,  therefore, 
both  as  a  direct  observation  of  the  electromagnetic  response 
of  a  given  cavity  and  as  a  rigorous  test  of  the  validity  of 
the  integral  equation  technique. 


B.  Procedures  for  Measurements  Inside  Cavity  End  Plates 

It  is  desirable  to  perform  field  measurements  in  such 
a  manner  that  the  transmission  line  from  the  probe  to  the 
sensing  instrumentation  does  not  interact  with  the  subject 
field.  This  may  be  effectively  accomplished  by  placing 
the  probe  on  the  surface  of  a  highly-conducting  boundary 
and  allowing  the  transmission  line  to  exit  through  and  be 
shielded  by  the  conductor.  Since  the  cavities  under  investi¬ 
gation  were  constructed  entirely  of  brass,  several  options 
for  providing  shielded  access  to  the  cavity  interior  were 
available.  The  options  were  to  measure  at  the  inner  surfaces 
of  the  end  plates  or  on  the  walls  of  the  inner  or  outer 
coaxial  conductors.  Boundary  conditions  constrain  Ep  to  be 
zero  on  the  end  plate  surfaces  and  Ez  to  be  zero  on  the 
walls  of  the  coaxial  conductors,  therefore,  it  was  de¬ 
sirable  to  choose  a  combination  of  measurement  locations 
which  would  sample  all  field  components. 

The  center  conductor  for  the  smaller  cavity  set  is 
less  than  one  centimeter  in  diameter  which  renders  it  an 
impractical  probe  location.  Difficulty  of  machining  argued 
against  cutting  slots  in  the  outer  coaxial  wall.  Con¬ 
versely,  it  was  relatively  simple  to  cut  slots  in  the  cavity 
end  plates  for  probe  access.  These  slots  were  required  to 
be  along  a  radial  line  in  order  to  provide  no  interference 
with  the  axially-directed  currents  flowing  through  the 


129 


cavity.  These  radial  slots  allowed  the  measurement  of  both 
the  Ez  and  the  components  on  the  front  and  rear  end 
plates  of  the  cavities  with  the  front  of  the  cavity  being 
defined  as  the  end  which  contains  the  port  being  excited. 

Radial  slots  were  also  cut  in  the  end  plates  of  the 
larger  set  of  coaxial  cavities.  In  addition,  the  center 
conductor  is  large  enough  to  accomodate  a  moveable  probe, 
so  it  was  possible  to  cut  an  axial  slot  in  this  conductor 
which  extended  the  entire  length  of  the  cavity  and  beyond. 
The  axial  slot  in  the  center  conductor  allowed  the  measure¬ 
ment  of  the  E  and  KL  components. 

P  $ 

The  field  probes  were  designed  for  operation  over  a 
conducting  screen  and  their  size  was  chosen  to  achieve  an 
acceptable  balance  between  providing  a  measureable  signal 
and  imposing  minimal  pertubation  of  the  field  being 
measured.  The  construction  of  the  probes  is  briefly  de¬ 
scribed  in  Appendix  F .  The  lead  from  the  probe  in  use  was 
connected  to  one  input  port  (test)  of  a  network  analyzer. 
The  reference  input  port  of  the  network  analyzer  was 
connected  to  the  RF  generator  which  provided  the  excita¬ 
tion  for  the  cavity.  The  network  analyzer,  configured  in 
this  manner,  performed  as  a  vector  voltmeter  displaying 
relative  amplitude  and  phase  as  a  function  of  position  as 
the  probe  was  moved  along  the  slot.  A  schematic  diagram 
of  the  measurement  instrumentation  is  presented  in  Figure  36 


Field  component  measurements  were  made  for  six  dif¬ 
ferent  configurations  of  the  smaller  set  of  cavities  and 
for  three  configurations  of  the  larger  set.  As  indicated 
earlier,  only  two  field  components  were  measured  in  the 
smaller  set  of  cavities  while  all  three  components  were 
measured  in  the  larger  set.  In  the  following  sections  of 
this  chapter,  graphical  presentations  are  given  for  se¬ 
lected  measured  field  distributions  for  each  of  the  nine 
cavity  configurations.  In  each  graph,  the  field  distri¬ 
bution  calculated  from  an  integral  equation  solution  is 
given  for  comparison.  The  calculated  distributions  have 
been  obtained  assuming  an  excitation  resulting  from  the 
application  of  a  time-varying,  one  volt  potential  be¬ 
tween  the  coaxial  conductors  (See  Eq.30).  The  amplitude 
distributions  in  each  graph  have  been  normalized  twice; 
the  first  normalization  is  by  a  factor  composed  of  geo¬ 
metrical  and  excitation  constants  which  cancels  the  units 
and  the  second  is  the  division  of  all  the  calculated 
values  by  the  largest  value.  The  peak  value  of  the  measured 
amplitudes  is  then  set  equal  to  the  corresponding  calcula¬ 
ted  value  at  the  same  radial  or  axial  location.  The  re¬ 
mainder  of  the  measured  values  are  then  divided  by  this 
factor.  The  calculated  phase  is  plotted  relative  to  an 
assumed  phase  of  zero  degrees  at  the  first  point.  The 
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Generator  Frequency  Counter 


36-  Schematic  Diagram  of  Instrumentation  for  Interior  Field  Measurements 


measured  phase  is  set  equal  to  the  calculated  phase  at  the 
same  point  where  the  amplitudes  are  made  to  be  equal  and 
the  other  measured  phase  points  are  plotted  relative  to 
this  value. 

C.  Field  Interior  to  Simple  One-Port  Cavity 

Measurements  of  the  interior  field  components  were 

performed  for  two  sizes  of  one-port,  coaxial  cavities. 

The  E„  and  H  components  were  measured  on  both  cavities 
z  $ 

and  the  component  was  also  measured  on  the  larger 
cavity.  The  measurements  on  the  smaller  cavity  were  per¬ 
formed  for  the  frequency  band  400  -  2000  MHz  at  100  MHz 
intervals.  For  the  larger  cavity,  the  frequency  band  was 
400  -  1200  MHz  at  100  MHz  intervals.  There  was  great 
diffici lty  in  obtaining  useful  measurements  of  Ez  on  the 
rear  plate  of  the  smaller  cavity.  This  was  the  only 
measurement  set  during  the  entire  investigation  which  pro¬ 
duced  such  difficulty.  The  problem  was  based  on  the  ex¬ 
treme  mismatch  between  the  input  transmission  line  and 
the  cavity.  Essentially  all  the  incident  energy  was  re¬ 
flected  back  toward  the  generator  and  it  was  not  possible 
to  safely  produce  enough  power  from  the  generator  to  yield 
meaningful  field  distributions  on  the  rear  plate.  The 
measurements  made  on  the  front  end-plates  showed  far  more 
variation  with  frequency  than  the  values  measured  on  the 
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rear  plates.  The  latter  values  were  almost  frequency  inde¬ 
pendent.  Representative  distributions  of  E  and  H  on  the 

z  $ 

front  and  rear  end  plates  of  a  simple  one-port,  coaxial 
cavity  are  presented  in  Figures  37  -  40.  The  two  letters, 

LC ,  under  the  small  picture  of  the  cavity  indicate  that  the 
particular  measurement  was  performed  on  the  larger  cavity 
set.  Such  a  distinction,  with  SC  indicating  the  smaller 
cavity  set,  is  made  for  the  one-port  cavities.  All  two-port 
cavity  measurements  were  made  on  the  smaller  cavity  set. 

D.  Field  Interior  to  Simple  Two-Port  Cavity 

A  two-port  cavity  configuration  could  only  be  con¬ 
structed  from  the  smaller  cavity  set,  therefore  only  one  set 
of  such  measurements  was  made.  The  output  waveguide  section 
was  terminated  with  a  50-ohm  load  in  the  same  manner  as  for 
the  admittance  measurements  of  Chapter  vil.  The  frequency 
range  of  these  measurements  was  400  to  2000  MHz  in  100  MHz 
steps.  There  was  no  problem  in  feeding  energy  into  this 
cavity  configuration  and  the  measurements  were  made  with  no 
difficulty.  Distributions  of  E  and  H.  on  the  front  and  rear 
end  plates  of  a  simple,  two-port  coaxial  cavity  are  presented 
in  Figures  41  -  44.  The  magentic  field  measurements  show  some 
interesting  variations  which  are  predicted  quite  accurately  by 
the  integral  equation  calculations 


E.  Field  Interior  to  Partitioned,  One-Port  Cavity 

Partitioned,  one-port  cavities  were  constructed  from  both 
sets  of  cavities.  Two  partitions  for  each  set  with  different 
size  apertures  were  made,  thus  it  was  possible  to  form  two 
variations  of  each  one-port  cavity.  The  frequency  range  for 
the  smaller  set  of  cavities  was  400  -  2000  MHz  at  100  MHz 
intervals.  The  field  components  in  the  larger  set  were 
measured  from  400  -  1200  MHz  at  100  MHz  intervals  with  a  few 
additional  measurements  at  1300  and  1400  MHz.  Distribtuions 
of  E  and  H  on  the  front  and  rear  end  plates  of  several 

Z  $ 

partitioned,  one-port  cavity  structures  are  presented  in 
Figures  45  -  52.  The  distributions  do  not  differ  markedly 
from  those  of  the  simple,  one-port  cavity. 

F.  Field  Interior  to  Partitioned,  Two-Port  Cavity 

It  was  only  possible  to  construct  two-port,  partitioned 
cavities  from  the  smaller  set  of  cavities,  but  by  making  two 
partitions  with  different  size  apertures,  two  variations  of 
such  a  cavity  were  formed.  The  frequency  range  of  these 
measurements  was  400  -  2000  MHz  at  100  MHz  intervals.  Dis¬ 
tributions  of  Ez  and  H^  on  the  end  plates  of  both  cavity  con¬ 
figurations  are  presented  in  Figures  53  -  68.  The  amplitude 
and  phase  curves  in  many  of  these  figures  have  a  wide  variety 
of  shapes  that  are  predicted  very  accurately  by  the  integral 
equation  solutions. 


G.  Measurements  Along  Center  Conductor  of  One-Port  Cavity 
The  most  interesting  measurements  in  terms  of  variety 


of  waveforms  were  made  at  the  surface  of  the  inner  conductor 
of  the  larger  cavity  set.  The  transverse  field  in  a  coaxial 
waveguide  has  its  greatest  intensity  at  the  surface  of  the 
inner  conductor  and  this  fact  appeared  to  enhance  the  quali¬ 
ty  of  the  measurements  made  there.  The  spatial  range  over 
which  the  measurements  were  made  was  greater  than  for  the 
end-plate  measurements  which  added  to  the  possibility  of  in¬ 
teresting  field  variations.  The  measurements  were  made  for 
the  simple,  one-port  cavity  and  .the  partitioned,  one-port 
cavity  with  two  sizes  of  apertures.  The  size  of  the  aper¬ 
ture  in  a  particular  figure  is  designated  by  the  letters  SA 
for  small  and  LA  for  large  immediately  following  the  LC  des¬ 
ignation  for  the  larger  cavity  set.  Distributions  of  and 
are  presented  in  Figures  69  -  78  for  all  three  cavity 
configurations.  The  integral  equation  solutions  do  an  excel¬ 
lent  job  of  predicting  the  amplitude  and  phase  in  all  cases. 

H.  Evaluation  of  Measured  and  Calculated  Results 

The  primary  objective  of  the  field-measurement  program 
was  to  provide  verification  for  the  accuracy  and  demonstrate 
the  utility  of  the  integral  equation  approach.  Sufficient 
measurements  were  made  to  provide  a  very  good  description  of 
the  field  in  the  various  cavities  in  the  frequency  range 
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which  extends  from  below  to  above  the  cutoff  frequency  for 
the  first  higher-order  mode  in  the  cavity  section.  Space 
restrictions  prevent  the  presentation  of  more  than  a  rep¬ 
resentative  sample  of  the  entire  body  of  measurements. 

A  small  number  of  the  measurements  display  evidence  of 
excessive  noise  and  the  strong  probability  of  unwanted  fre¬ 
quency  harmonics.  It  was  discovered  during  the  investiga¬ 
tion  that  the  generator  plug-in  used  in  the  400  -  1000  MHz 
band  produced  substantial  frequency  harmonics  until  it  had 
warmed  up  for  an  extensive  period.  Most  of  the  measurements 
where  this  effect  was  suspected  were  repeated,  but  several 
of  the  measurement  sets  have  unexplained  variations  from  the 
calculated  values  below  1000  MHz  but  have  very  good  agree¬ 
ment  above  1000  MHz. 

The  great  majority  of  the  measurements  display  very 
good  to  excellent  agreement  with  the  calculated  field  dis¬ 
tributions  throughout  the  measured  frequency  band.  Many 
times  the  agreement  follows  rather  complex  variations  of 
both  amplitude  and  phase.  The  overall  impact  of  the  com¬ 
parisons  is  a  very  strong  verification  that  the  integral 
equation  approach  provides  accurate  solutions  for  the  aper¬ 
ture  fields  and,  consequently,  all  other  field  components 
in  the  cavities  that  have  been  treated  during  this  investi¬ 
gation  . 
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riGURE  38  MEASURED  AND  CALCULATED  Z  COMPONENT  OT  ELECTRIC  HELD  800  nhZ 

1  38 


a  b  RflOIUS  NORMALIZED  BT  c  e 

FIGURE  39  MEASURED  ANO  CALCULATED  «  COMPONENT  OF  MAGNETIC  FIELD  BOO  MHZ 
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FIGURE  «  MEASURED  ANO  CALCULATED  •  COMPONENT  OF  MAGNETIC  FIELD  900  MHZ 
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FIGURE  52  MEASURED  ANO  CALCULATED  •  COMPONENT  OF  MAGNETIC  FIELD  600  MHZ 
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RELATIVE  PHASE  ANGLE  IN  DEGREES  J  RELATIVE  PHASE  ANGLE  IN  OCGREES 


LD  NORMAL  I Z 


PEAK  AMP  -  .057  a  -  .0121  >  o  •  .0286  .  c  *  .4246  . 


r I GORE  57  MEASURED  ANO  CALCULATED  2  COMPONENT  0E  ELECTRIC  TIELD  1200  I 


PEAK  AMP  -  .067  a  -  .0134  i  o  -  .0310  >.  c  •  .4600  k 


El  GORE  58  MEASURED  ANO  CALCULATED  Z  COMPONENT  OT  ELECTRIC  HELO  1300  MHZ 

7  48 

t 


RELATIVE  PHASE  ANGLE  IN  OEGREES  i  RELATIVE  PHASE  ANGLE  IN  DEGREES 


AMPLITUDE  Of  E-flELO  NORMALIZED  BT  V/o«LOGIb/ol  ->  AMPLITUDE  Of  E-flELD  NORMALIZED  8T  V/o-LOGlb/ol 


AMPLITUDE  Gf  H-EIELO  NORMALIZED  ai  V/Z-o-LOGIb/al  ^  AMPLITUDE  0 T  H-EIELD  NORMALIZED  Bl  V/Z-o-LDGlb/ol 


PEAK  AMP  -  .SO  a  -  .01 14  v  b  -  .0262  >  c  -  .3893  > 


1GURE  61  MEASURED  ANO  CALCULATED  •  COMPONENT  0T  MAGNETIC  HELD  1100  mmZ 
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riGURE  76  MEASURED  ANO  CALCULATED  •  COMPONENT  OT  MAGNETIC  HELD  1100  MHZ 
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RELATIVE  PHASE  ANGLE  IN  DEGREES  3  RELATIVE  PHASC  ANGLE  IN  QEGRECS 
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FIGURE  78  MEASURED  ANO  CflLCULflTED  «  COMPONENT  OF  MA6NETIC  HELD  1200  MHZ 
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CHAPTER  IX 


CONCLUSIONS 


The  primary  objective  of  this  investigation  has  been 
to  develop  a  proven  and  practical  technique  for  determin¬ 
ing  the  electromagnetic  properties  of  a  set  of  planar  dis¬ 
continuities  in  coaxial  waveguides.  It  is  demonstrated 
early  in  this  report  that  a  knowledge  of  the  transverse 
electric  field  distribution  in  the  plane  of  the  discon¬ 
tinuity  (or  discontinuities)  is  sufficient  to  characterize 
all  electromagnetic  properties  of  the  structure.  According¬ 
ly,  the  bulk  of  the  analysis  has  dealt  with  the  derivation 
of  tractable  integral  equations  which  can  be  solved  to  ob¬ 
tain  these  desired  aperture  field  distributions.  The 
technique  for  deriving  such  equations  has  proved  to  be 
quite  general  and  has  been  easily  extended  to  the  treat¬ 
ment  of  coupled  discontinuities.  The  general  approach  has 
also  been  applied  to  the  related  problem  of  the  discontinuity 
occuring  at  the  junction  of  a  coaxial  and  circular  waveguide 
where  the  coaxial  center  conductor  is  truncated. 

A  significant  outcome  of  this  effort  is  that  any  of  the 
coaxial  junctions,  whether  single  or  one  of  several  forming  a 
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complex  structure,  can  be  treated  through  the  application  of 
one  general  equation  form.  Use  of  such  a  general  form  is 
frequently  cumbersome  and  experience  has  shown  that  three 
somewhat  less  general  equation  forms  are  more  practical. 

The  first  of  these  general  forms  treats  isolated  discontinu¬ 
ities,  the  second  treats  discontinuities  which  are  coupled 
to  one  other  discontinuity,  and  the  third  treats  discontinu¬ 
ities  which  are  coupled  to  two  other  discontinuities,  one  on 
either  side.  The  general  equation  form  for  the  truncated 
coaxial  center  conductor  may  be  similarly  structured  with  an 
obvious  limitation  on  the  types  of  discontinuities  to  which 
this  discontinuity  may  be  simultaneously  coupled. 

The  integral  equations  developed  during  this  investi¬ 
gation  have  proved  to  be  quite  amenable  to  a  straightforward 
numerical  solution  technique.  The  choice  of  an  expansion 
function  set  for  the  aperture  fields  which  allows  analytic 
evaluation  of  the  integrals  significantly  reduces  the  number 
of  machine  operations  required  and  thus  reduces  the  effect 
of  this  unavoidable  source  of  error.  Another  very  favorable 
factor  is  that  the  primary  electromagnetic  quantities  of 
interest,  e.g.,  admittance  and  internal  field  components, 
are  computed  from  the  integrals  of  the  aperture  fields,  and 
integration  is  an  operation  which  tends  to  reduce  the  effects 
of  many  distribution  errors  rather  than  enhance  them. 

A  significant  amount  of  effort  during  this  investigation 


was  devoted  to  experimental  measurements  designed  to  verify 
the  accuracy  of  the  integral  equation  formulation  and  the 
numerical  solution  techniques.  The  comparisons  between 
measured  and  calculated  admittance  values  for  three  cavity 
structures  demonstrated  extremely  good  agreement.  Signi¬ 
ficant  differences  between  measured  and  calculated  values 
appeared  only  at  the  highest  frequencies  measured  and  it  is 
very  probable  that  these  differences  are  the  result  of 
experimental  error.  It  is  quite  interesting  to  observe  the 
exceptional  agreement  of  the  admittance  calculated  from 
transmission  line  theory  with  susceptance  corrections  from 
the  Waveguide  Handbook  at  lower  frequencies,  and  then  to  ob¬ 
serve  the  fairly  rapid  breakdown  of  this  method  as  the  fre¬ 
quency  approaches  the  cutoff  for  the  first  higher-order  mode 
The  measurements  of  internal  field  components  for  nine 
different  cavity  structures  and  subsequent  comparisons  with 
calculated  distributions  has  provided  very  impressive  veri¬ 
fication  of  the  ability  of  the  analytical  technique  to  per¬ 
form  such  calculations.  There  is  not  exact  agreement  for 
every  measurement,  but  there  is  strong  evidence  that  all 
such  disagreements  result  from  either  instrumentation  or 
experimental  judgement  errors.  There  are  no  cavity  con¬ 
figurations  where  there  is  not  exceptional  agreement  between 
measured  and  calculated  field  distributions  for  a  majority 
of  the  frequencies  sampled.  The  field  distributions  them- 
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selves  provide  an  interesting  study.  There  is  generally 
little  variation  in  the  field  distributions  except  near  the 
transverse  resonance  frequencies.  There  is  very  little 
variation  even  at  this  frequency  for  measurements  made  at 
the  rear  end-plates  of  the  one-port  cavities.  Measurements 
made  on  the  end-plates  of  the  two-port  cavities  generally 
display  a  gradual  phase  variation  while  the  only  phase  varia¬ 
tions  observed  for  one-port  cavity  measurements  are  abrupt 
180  degree  shifts.  If  the  primary  goal  of  this  effort  had 
been  to  perform  a  rigorous  chacteri2ation  of  several  cavity 
structures,  a  substantial  amount  of  additional  data  could 
have  been  generated  and  interesting  comparisons  drawn  between 
field  and  admittance  measurements  over  the  frequency  band. 
However,  the  amount  of  data  obtained  is  more  than  adequate 
to  verify  the  capability  of  the  integral  equation  approach. 

The  descriptive  term  "practical",  as  applied  in  the 
first  sentence  of  this  chapter,  can  be  interpreted  to  mean 
either  useful  or  n^t  excessively  complicated  to  use,  or, 
more  likely,  a  given  individual's  combination  of  the  two. 

A  few  suggested  applications  for  the  results  of  this  investi¬ 
gation  are  given  in  Chapter  I  but  this  list  does  not  emphasize 
possible  microwave  circuit  applications.  The  demonstrations 
of  the  accuracy  of  the  technique  indicate  possible  appli¬ 
cations  in  a  large  area  of  waveguide  design  problems  where 
the  expressions  in  the  Waveguide  Handbook  have  been  used  in 
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the  past.  The  integral  equation  technique  has  given  evidence 
of  being  more  accurate  for  susceptance  calculations  than  the 
expressions  in  the  Waveguide  Handbook  for  apertures  that  are 
small  compared  to  the  largest  radial  dimension.  Examples 
presented  in  this  report  indicate  that  the  accuracy  of  the 
aperture  field  solutions  and,  thus,  the  admittance  calcula¬ 
tions  can  be  improved  by  increasing  the  number  of  expansion 
pulses  used  up  to  an,  as  yet,  undetermined  number.  In¬ 
creasing  the  number  of  series  terms  alone  appears  to  have  a 
less  direct  effect  on  convergence,  but  a  larger  number  of 
pulses  does  require  a  larger  number  of  series  terms  to  achieve 
convergence.  There  is  very  likely  a  set  of  coaxial  wave¬ 
guide  problems  where  either  the  requirements  for  extreme 
accuracy  or  extended  bandwidth  make  the  subject  approach 
extremely  appropriate.  None  of  the  publications  listed  in 
the  Bibliography  offer  the  means  for  calculating  the  admit¬ 
tance  of  coaxial  discontinuities  when  higher-order  modes  are 
capable  of  propagation.  The  method  of  Wexler,  Ref. [22], 
which  can  treat  higher-order  modes,  could  be  modified  to 
treat  coaxial  geometries;  however,  it  is  very  cumbersome  to 
treat  more  than  a  few  such  modes.  It  is  the  hope  of  the 
authors  that  the  existence  of  this  technique  will  play  a  role 
in  developing  useful  applications  for  which  it  is  well  suited. 


1  63 


APPENDIX  A 


COMPLETENESS  OF  FIELD  DERIVED  FROM  A  MAGNETIC  VECTOR  POTENTIAL 
FOR  A  ROT AT I ONALLY  SYMMETRIC  WAVEGUIDE 

All  planar  discontinuities  in  the  coaxial  waveguides  con¬ 
sidered  in  the  present  investigation  are  surfaces  which  sup¬ 
port  a  transversely-directed  electric  current.  The  existence 
of  such  current  components  implies  that  a  field  description 
derived  exclusively  from  a  magnetic  vector  potential  with 
only  an  axial  component  might  not  be  complete.  It  is  known 
that  a  field  description  which  consists  of  the  linear  com¬ 
bination  of  all  possible  TEM,  TM,  and  TE  modes  in  a  cylindri¬ 
cal  waveguide  is  complete.  In  the  present  investigation, 
the  TEM  and  TM  modes  are  included  so  it  remains  to  be  deter¬ 
mined  whether  or  not  there  is  possibly  a  necessary  TE  mode 
contribution  to  the  field. 

Within  a  cylindrical  waveguide,  all  TE-mode  field  com¬ 
ponents  can  be  derived  from  the  axial  component  of  an  elec¬ 
tric  vector  potential.  The  electric  vector  potential  must 
satisfy  the  wave  equation  and  for  a  coaxial  region  the  general 
form  is  (note  that  the  m  =  0  term  is  not  included) 

oo  QO 

F  ( p , <J>  ,z )  =  l  l  am„  Z  ( z )  eDmtl)  E  ( p )  (A-l) 

z  _  t  mn  mn  mn 

m=-  n=i 
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where  the  axially-varying  function  Z(z)  is  a  solution  of  the 
harmonic  equation 


and  5mn  represents  a  set  of  solutions  of  Bessel's  equation  of 
order  m.  The  form  of  Emn  is  dependent  upon  the  requirement 
that  the  tangential  component  of  the  transverse  electric  field 
vanish  at  the  walls  of  the  waveguide.  For  a  rotationally- 
symmetric  waveguide  with  symmetric  excitation,  m  must  be  zero 
and  Eq. (A-l)  becomes 


Fz(p,z)  =  l  an  Zn(z)  E n ( p )  .  (A-3 ) 

n=l 

The  field  components  derived  from  this  electric  vector 
potential  are 


For  the  field  composed  of  TEM  and  TM  modes  to  be  incom¬ 
plete,  the  p-directed  current  on  the  transverse  surface  must 
generate  a  field  component  which  could  be  derived  only  from 
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an  axially  directed  electric  vector  potential.  A  $-inde- 
pendent,  p-directed  electric  current  generates  only  a  p 
component  of  magnetic  vector  potential,  A,,  and  the  field 
derived  from  A  is  given  by 


E  =  - 


j  [*2PA  +  7(7*  pA  ) ] 

]c  °  P 


•  U)  r  «  4*  /n  _ 

-  -D  -j  [*  PA, 
k 


1  3 


(A-6 ) 


and 


H  = 


-  7xgA 
U  p 


1  3A„ 

-  -  $  — jr-2. 

p  ^  3z 


(A-7 ) 


It  should  be  noted  that  the  scalar  quantity  is  not  a 
solution  to  the  wave  equation  and,  therefore,  is  not  useful 
in  obtaining  a  solution  for  the  field  inside  a  waveguide. 

It  is  useful  to  note  that  that  the  field  components  described 
by  Eqs . ( A-6 )  and  (A-7)  are  precisely  the  same  as  those  ob¬ 
tained  from  a  ^-independent ,  axially  directed  component  of 
magnetic  vector  potential,  Az .  Furthermore,  neither  of  the 
field  components  described  by  Eqs.(A-4)  and  (A-5)  has  a  vec¬ 
tor  component  that  is  common  to  the  same  field  component  of 
Eqs. (A-6)  and  (A-7).  Therefore,  there  is  no  contribution  to 
the  field  by  a  postulated  Fz,  so  none  is  required  for  a  com¬ 
plete  field  description. 


APPENDIX  B 


PROCEDURE  TO  IMPROVE  CONVERGENCE  RATE  OF  BESSEL  FUNCTION  SERIES 

An  often  productive  procedure  to  increase  the  convergence 
rate  of  a  given  series  consists  of  the  termwise  subtraction 
of  a  second  series  having  a  known  sum  and  whose  higher-order 
terms  approach  and  eventually  become  equal  to  the  corres¬ 
ponding  terms  of  the  first  series.  The  sum  of  the  initial 
series  is  then  computed  by  adding  the  known  sum  of  the  second 
series  to  the  series  of  difference  terms.  To  illustrate  the 
procedure,  assume  that  the  sum  F  of  the  following  series  is  de¬ 
sired  : 

oo 

F  =  y  f  .  (B-l) 

i,  n 
n=l 

Next  assume  that  there  exists  a  second  series  with  a  known 
sum  G 

00 

G  =  l  g  (B-2) 

n=l 

whose  individual  terms  have  the  following  property 

lim  (fn  -  gfi)  =  0  .  (B-3 ) 

n-*-® 

If  both  series  are  convergent,  then  their  sum  or  difference 


is  convergent,  and  Eq. (B-l)  may  be  expressed 

00 

F-  l  (f  -  g  )  +  G  .  (B-4) 

n=l 

The  difference  series  converges  much  more  rapidly  than  either 
of  the  individual  series.  The  rate  of  convergence  is  changed 
from  a  dependence  upon  how  rapidly  the  magnitude  of  f 
approaches  zero  for  increasing  n  to  a  function  of  how  rapidly 
the  individual  f  approaches  the  value  of  the  corresponding 


APPENDIX  C 


ASYMPTOTIC  FORMS  FOR  EIGENFUNCTION  EXPANSIONS  OF  FIELD 
COMPONENTS  IN  COAXIAL  AND  CIRCULAR  WAVEGUIDES 


A.  Introduction 

In  order  to  apply  the  procedure  described  in  Appendix  B 
to  the  Bessel  function  series  appearing  in  the  definitions  of 
the  field  components  in  coaxial  and  circular  waveguides,  the 
behavior  of  the  series  elements  must  be  determined  for  in¬ 
creasing  n.  The  roots  of  Bessel  functions  of  the  first  and 
second  kinds  and  of  the  function 


Vp)  =  N0(Yna)J0(Ynp)  _  J0(Yna)N0(7np) 


(C-l) 


which  is  frequently  obtained  for  coaxial  geometries  (with 
inner  radius  a),  are  increasing  functions  of  n;  i.e., 

Yn+^  >  Yn  for  n+1  >  n,  where  y denotes  the  nth  root  of  the 
particular  function.  Since  the  value  of  Yn  increases  rapidly 
with  n,  the  large  argument  forms  of  the  Bessel  functions 
become  applicable.  The  following  forms  are  used: 


J 

N 


0 

0 


(Y„0> 

<V” 


/_2_ 

n+co  J  tty  P 
v  n 


cos(YnP  -  J) 
sin(YnP  -  J) 


(C-2a) 


(C-2b) 


169 


(C-2c) 


•VV  TZS-fiki  cos(-<n»  -  r-> 

▼  n 

N.  (yj)  - — ►  /— —  sin  (y  p  -  Js-) 
1  n  n-H»  tty  p  1  nK  4 

^  n 


(C-2d) 


B.  Asymptotic  Form  for  in  Coaxial  Waveguide 

With  the  assumption  of  no  variation  with  4> ,  the  trans¬ 
verse  component  of  magnetic  field  in  a  coaxial  waveguide  is 
given  by  a  positive  and  a  negative  traveling  TEM-mode  wave 
and  an  infinite  sum  of  higher-order  mode  terms: 


B 


+  e 


-jkz 


HP 


-  B 


-  e 


jkz 


np 


+  $  A  °'lv  -  e3“nZ)  £ 

n=l  n  dp 


where  $n  and  an  are  defined  by  Eqs. (20)  and  (21),  respec¬ 
tively.  In  this  exercise,  interest  is  limited  to  the  higher- 
order  mode  terms  which  make  up  the  infinite  series.  For  a 
coaxial  region  between  two  discontinuities  located  at  z  =  s 
and  t,  t  >  s,  the  constants  are  defined 


+ 


.  £  ja  t 
±e  J  n 


j2sina  (t-s)  M„ 
n  n 


P  d<1>n  dp 

ar* 


(C-3) 


where  the  aperture  electric  field,  E°,  is  defined  by  Eq. (31) 
2 

and  Mn  is  defined  by  Eq. (36) .  The  expression  for  the  total 
TM  contribution  to  the  field  is  given  by 


30 

V 

L 

n»l 


cosa  (t-z) 
n 


jsina  (t-s)  all 
n  n  n 


E°  P  d$n  dp 
h  p  dp 


d* 

n 

d T 


(C-4 ) 
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The  aperture  electric  field  may  be  expanded  in  a  set  of  N 
pulses  as  defined  by  Eqs.{54)  and  (55).  The  integral  is 
then  evaluated  to  obtain 


rTM  _  k  j 


cosa  (t-z) 
n 


n=l  jsina  (t-s)  a  M 
n  n  n 


N 


[  E„  [*  <r+)  -  *  (tjl 


qil  9  '  "  9 


n  q' 


d$ 

_ i 

dp 


(05) 


where  r^  and  r^  are  defined  by  Eq.  (61) .  The  expression  may 
be  rearranged  to  the  form 


H™  -  *  l  E 
*  n  q=l  ^ 


N 


00  cosa  (t-z)  4>  (r  )  -  $  (r  )  d$ 

z  _ n _  n  q _ n  q  n 


n=l  jsinan(t-s) 


a  M2 
n  n 


Attention  will  now  be  directed  toward  the  expression  within 
the  braces  which  is  represented  by  Rn(p): 


Rn(P)  - 


4>_  (r*) 
_ 3_ 


-  $  (r  )  d<? 
n  _3 _ Q. 


2 

a  M 
n  n 


dp 


(06) 


It  is  the  expression  given  by  Eq. (C-6)  which  determines  the 

asymptotic  behavior  of  the  terms  in  the  infinite  series  of 

Eq.  (C-5) .  The  elements  of  Eq.(C-6)  are  expanded  in  Bessel 

2 

functions  through  the  use  of  the  definition  for  and  the 
following 


dO  (P) 

Sq—  “  -WYna)Jl{Ynp)  "  J0(Yna)Nl(Ynp)] 


(07) 
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Vrq>  *  Vrq>  * 

'  tVV’WV  ‘  VYi^Wq’l.  (C-8) 

Substitution  of  the  Bessel  function  forms  into  Eq. (C-6) 
yields  the  following  expression: 

R_<p>  -  -VV^'VV11  *  J0<Yn1,1,l(V,> 

n  2  2  ~  ~  - - 

tW)Jl(V»  -  J0lYna)N1,’'nb,J2  ‘  ^ 

7T  > 

f[N0(Yna)J0(v„r*)  -  W’Wq’l 

■  ■W’VVV  '  W’Wq’lt  (C-9) 


The  asymptotic  forms  of  Eq.(C-2)  are  now  substituted  into 


Eq.  (09)  to  yield 


VP)  = 


- <JL=Z  [sin  (y  a-  J)cos(y  p~)  -  cos  (y  a-J)  sin  (y  p~)  ] 

a  tty  /ao 

n  n 

4b^y^ 

— 2_[sin(Yna-J)cos(Ynb-^)  -  cos  (Y„a-J)  sin  (Ynb-|l)  ]  2  -  % 

2abTf  y  n 

n 


21  lr.,  TT  .  .  +  TT ,  ,  TT .  .  ,  +  TT  .  , 

WIT!  7?[Sln(Yna'4>CC'S(1,nVT)  *  =os(Yna-T)sin(Ynr  -j)] 
n  r 

v  q 


7=f[sin(yna-?)cos(ynr"-j)  -  cos  ( Yna-J)  sin  (Ynr"-J)  ]  .  (OlO) 
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The  trigonometric  identity 


sin4>  cos©  -  cos0  sin9  =  sin  ($-6) 


(C-ll) 


is  employed  to  reduce  Eq.  (C-10)  to  the  form 


-4sin[vn (a-o) +j] 


^-^rsin[Yn(a-rg)  }  -  _L_Sin[Yn  (a-r‘)  ] 


Vp)  - 


u2a/p  Ynan  h  I  sin2[Yn(a-b)-J]  -  1 

7T  v 


The  expression  may  be  further  reduced  by  noting  that 


sin($  +  4-)  =  cos  <p 


which  yields 


-2cos [ Yn (a-p ) ] 


Vo)  - 


f_i_sin[Yn(a-rJ)]  -  -^sin[Yn  (a-r’)  ]  j 

rq  rq  J 


a/p  Ynan  l[b 


cos  [Yn (a-b) ]  -  a 


As  n  becomes  larger,  an  approaches  the  value  -jy  ,  yieldina 


-2cos [Yn (a-p ) ] 


Rn(P)  « 


1  +1  - 
— sin[Yn(a-rg) ]  -  — sin [yr <a-rq > ] 

- 


-jyn  ^ 


b  cos  [y  (a-b) ]  -  a 


As  n  becomes  larger,  y  approaches  the  value  The  form 

n  “  b-a 

of  the  expression  then  becomes 
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rniT 


-j2c°s[j^j(a-c)  ] 


Rn(P)  - 


^sin[£_(a-rg)]  - 

r 

q 


-tsintgfjfa-^)] 


2 

/o  [r— I  ^b  cos2[£-(a-b>  ]  -  al  . 
[b-aj  lb-a  J  j 


This  expression  may  be  simplified  through  the  use  of  two 
trigonometric  identities: 


and 


cos  (-nTr)  =  1 


cose  sin6  =  *5  sin(C+5)  -  sin(f-r)  . 


The  expression  now  becomes 


...  i  [  sin  (nS ,  )  sin(n£,)'i 

v-->  *  -h  \—r - 

i  ►  P  ►  +  (  n  n  ) 


fsinfnS^)  sinfnC^)] 

2  1 
n  J 


v  -  i  n 

rq 


(C-12 ) 


where 


91  *  b^(2a  “  *  "  rq> 


_  IT  .  + 

0,  =  vTT’(rr,  "C) 
2  b-a  q 


9 3  "  b^I(2a  -  D  -  V 


-  II  .  —  » 

9 4  =  b^(rq  '  0)  * 


This  exercise  has  demonstrated  that  the  asymptotic  (large 
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argument)  form  of  the  Bessel  function  expansion  for  H ' 
consists  of  four  series  which  have  the  general  form: 


SO)  =  l  5i£* 
n=l  n 


sin (n9) 


(C-13) 


In  order  to  apply  the  technique  described  in  Appendix  B, 
another  series  is  required  which  has  a  known  sum  and  whose 
individual  elements  have  a  form  similar  to  those  of  Eq.(C-13). 
In  Ref. [26]/  expression  604  is  given  as 


j  sin (nS) 
n=2  n2-l 


-  sinG  ln^sinj}  0  <  l  <  2~  .  (C-14) 


The  procedure  of  Appendix  B  may  now  be  employed  to  obtain  a 
modified  form  of  Eq.  (C-12)  which  converges  more  rapidly. 

The  resulting  series  expansion  for  is 


„TM  .  k  8  ?  eos°n(t-2)  fVrq>  -  Vrc>  dV-M 

4  n  q-1  “3  nil  3sinan(t-s)  0  „2 

{  n  n  j 


j(b-a)  1  fsinUSi)  -  sin(nS2) 


sinCnO^)  -  sin(n£<})'l 


tt2/p  n2-l 


.  2  (b-a)  _L_ 

ir2/p  /~+ 

rq 


1  ® i  i  ®2  ) 

sinO^  ( j-ln{  2sin^— } )  -  sin©2  ( “ln{  2sin^— } )  j 


if  i  0-5  -l  0 . 

-  sin0,  ( j-ln{  2siny-} )  -  sin0 ,  ( j-ln{  2sin-=— } ) 


(C-15 ) 


Because  the  series  of  Eq.(C-14)  begins  with  n  =  2,  the  asymp¬ 
totic  terms  are  not  evaluated  for  n  =  1.  Care  must  be  taken 
that  the  arguments  of  the  sine  functions  in  Eq.  (C-15)  be  in 
the  range  specified  by  Eq.(C-14).  If  the  argument  falls  out¬ 
side  the  range,  the  argument  can  be  brought  into  the  proper 
range  by  adding  or  subtracting  multiples  of  2~  or  by  employ¬ 
ing  appropriate  trigonometric  relationships. 

The  sine  and  cosine  functions  in  Eq.  (C-5)  are  not  in¬ 
cluded  in  the  derivation  of  the  asymptotic  form  for  two  rea¬ 
sons:  the  first  is  that  these  terms  are  not  always  present 

and  the  second  is  that  when  the  terms  are  present,  they  are 
easily  treated  because  of  the  behavior  of  hyperbolic  func¬ 
tions.  The  behavior  of  the  quotient  cosa^ (t-z) /jsinan  (t-s) 
as  n  increases  is  primarily  dependent  on  the  value  of  z. 

Since  an  approaches  -jy  for  large  n,  the  trigonometric 
functions  become  hyperbolic  functions  with  the  following 
behavior : 


For  z  =  s : 


lim  =°=an<t-s) 


n-*-<» 
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coshy  (t-s) 
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-jsinhYn (t-s) 


(C-16 ) 


For  z*;:  s<C<t 
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An  inspection  of  Eq. (C-5)  will  verify  that  the  asymptotic 
form  of  Eq. (C-16)  normally  occurs  in  conjunction  with  a  y 
factor.  Since  the  product  of  j  and  i-  is  unity,  the  limit 
of  the  quotient  cosa^ (t-s) /jsinan (t-s)  as  n  increases  with¬ 
out  bound  is  unity  and,  therefore,  it  does  not  enter  into 
the  asymptotic  form  for  .  An  even  more  advantageous  re¬ 
sult  is  obtained  when  the  quotient  is  of  the  form  given  by 
Eq. (C-17) .  These  factors  approach  zero  so  rapidly  with  in¬ 
creasing  n  that  it  is  unnecessary  to  employ  the  procedure 
of  Appendix  to  obtain  rapid  convergence. 


C.  Asymptotic  Form  for  Ep  in  Coaxial  Waveguide 

An  expression  for  the  higher-order  mode  expansion  of  E. 
in  the  same  coaxial  region  for  which  Eq .(C-4)  is  defined  is 
given  by 


-  I 


sina  (t-z) 
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n=l  sina  (t-s)  M 
n  n 


fb_  o  d*  .  d$ 

j  E  p  _ n  dp  _ n 

;  Jh  p  dp  dp 


(C-18 ) 


The  only  difference  between  Eq.  (C-5)  and  Eq. (C-18)  which 

would  affect  the  asymptotic  form  of  the  latter  is  the  lack 

of  the  1/a  term  found  in  the  former.  The  use  of  this  fact 
n 

allows  the  derivation  of  the  asymptotic  form  Sn  ( p )  equiva- 

TM 

lent  to  Rn(p)  for  H  by  multiplying  Eq. (C-12)  by  the  asymp¬ 
totic  form  of  a  : 
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lim  a 


-  jmr 

b-a 
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t 


The  resulting  expression  for  Sn(p)  is  given  by 


(P) 


-1 

’T  /p 


f  sin (n-S^) 


f sin (n©^ ) 


sin (n©2) 
n 

sin (n©^ ) 
n 


(C-19) 


The  definitions  of  the  ©n  are  the  same  as  those  given  for 
Eq. (C-12) .  The  convergence  improving  procedure  of  Appendix 
B  can  now  be  applied  by  employing  expression  number  508  of 
Ref . [26]  : 

oo 

l  -S--nj— =  *5  (tt-9  )  ,  0  <  0  <  2*  .  (020) 

n*l  n 


D.  Asymptotic  Form  for  E  in  Coaxial  Waveguide 

z 

The  z-directed  component  of  electric  field  has  no  TEM 
component  and  is  defined  entirely  by  its  higher-order  mode 
expansion.  An  expression  for  Ez  which  is  valid  in  the  same 
coaxial  region  as  that  for  Eqs.(C-4)  and  (C-18)  is  given  by 


Ez  "  l. 

n=l 


cosa  (t-z) 
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jsina  (t-s)  <*nMn 


E°  p  ^n  dp 
h  0  dp 


(021) 


The  process  for  obtaining  the  asymptotic  form  is  begun  by 
substituting  the  large  argument  forms  of  the  Bessel  functions 
into  Eq.  (021)  and  results  in  the  expression: 


178 


¥ 


.  ’nfo'V^o'V1  -W'W’J 

Tn (p)  =  22 

IW»Jl<'rnb>  -  J0(''na)VYnbl>2  - 


n 
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(C-22 ) 


After  the  steps  outlined  in  Eqs.(C-7)  through  (C-12)  are 


followed,  an  asymptotic  expression  for  Tn(p)  is  obtained: 
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(C-23 ) 


The  0  are  the  same  as  those  defined  for  Eg. (C-12) .  The 
n 


convergence  improving  procedure  of  Appendix  B  may  be  ap¬ 
plied  with  the  use  of  expression  number  503  of  Ref. [26]: 
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n=l 


cos (n0) 
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-ln{ 2sinj} 


0  <  0  <  2tt 


(C-24 ) 


E.  Asymptotic  Form  for  in  Circular  Waveguide 

The  field  in  a  circular  waveguide  consists  only  of 
higher-order  modes.  The  field  expansion  is  different  from 
that  of  a  coaxial  region  and,  consequently,  the  asymptotic 


9 


forms  for  the  individual  field  components  are  different.  In 
Chapter  VI,  it  is  shown  that  the  transverse  component  of  mag¬ 
netic  field  is  proportional  to  the  following  expression: 
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(C-25) 


After  E°  is  expanded  in  a  set  of  weighted  pulses  defined  by 
Eq.  (31) ,  the  integral  is  evaluated  and  a  function  Qn(p)  is 
extracted: 


Qn(p)  = 


-2Jl(Xnp) 
Xnb\  Jl(Xnb> 


lJ0(XnV  -  Wq>]  * 


(C-26 : 


Substitution  of  the  asymptotic  forms  for  the  Bessel  functions 


into  Eq.  (C-26)  yields 
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The  trigonometric  identity 

cos4>  cos9  =  h[cos(4>+9)  +  cos(4>-9)]  -  (C-28) 


is  now  employed  to  yield 
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Three  additional  identies: 
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COS($-1T)  =  -COS$  , 


cos ($~j)  =  sin$  , 


COS^i})  =  *5[cOS2({)  +  1]  , 


(C-2  9a ) 
(C-29b) 
(C-29c) 


are  used  to  obtain 


-2 


Q  (p)  a  - — - - -  f-~  {“COS[X  (p+r^)] 

a  A  b  /p  [-sin(2A  b)  +  1]  /~+  q 

n  n  n  rg 


+  sin[Xn(p-rg)  ] }  -  { -cos [  Xn  ( p+r^)  ]  +  sin[  Xn ( p-rq)  ]  } 

rq 


For  increasing  n: 


n 


-jX 


n 


and 


(n-h)  tt 


(C-30a) 


(C-3  Ob) 


The  asymptotic  form  of  Qn(p)  is  now  given  by 
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-2 


X 


Q„(P)  - 


-jb  ln~‘i’21'2  S5 


k  (sin[i£^(p-rJ)J  -  cos[i^iii(p+r*)]( 


~  {sin[  .(n-~%)7T(p-r~)l  -  cos[---n~ig)  77  (p+r~)  ] } 

rq 


Note  that: 


sin  [  (2n-*s)iT]  =  -1 


The  asymptotic  expression  for  Q  (  )  may  be  presented  more 
compactly  by 


Qn(P) 


— ^ - ^ -i—  {  sin  (n->s)  -  cos(n-%)uJ 

r~+  x  z 

rq 

-  — {  sin  <n-*s)  $3  -  cos(n-*s)$4} 


(n-3s)2-rr2  /p 


t  (C 


where 


*  71  /  +x 

$1  =  E(p-rg)  , 


4,  = 


tt  ,  .  +  , 

T-(P+r  , 
b  q 


and 


4,  = 


E1P"rq>  ' 


4,  = 


b(P+rq>  * 


-31) 
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The  use  of  the  trigonometric  identies : 


and 


sin(4>-0)  =  sin$  cos8  -  cosiji  sinQ 
cos{(J>-0)  =  cos<j)  cos8  +  sin<j)  sin0  , 


(C-32 ) 
(C-33) 


results  in  the  expression 


Qn(p)  " 


-lib 


(n-^s)  2tt2/p 


If  <J>  (t) 

-  {cos  1  sin(n4>.)  -  sin  1  cos(n<f,) 

/“+  4~  1  T  1 


'  rq 


-  cos^2  cos(n<f_)  -  sin^2  sin(n$-)}  -  — {cos^3  sin(n<t>,) 
4^4^  S~^  4  J 

rq 


-  sin<t>3  cos(n<t.,)  -  cos^4  cos(n$  .)  -  sin^4  sin(n<I>  )} 
4  J  4  4  4 


A. 


.  (C-34) 


The  procedure  of  Appendix  B  may  now  be  applied  by  employing 
Eq. (C-14)  for  the  sine  series  and  one  of  the  following  ex¬ 
pressions  from  Ref. [26]:  (The  second  expression  is  #605.) 


00 


l 

n=l 


cos  (n<t ) 

2 

n 


0  <  4  <  2tt 


( C—  3  5 ) 


l  COy(-n<i>)  =  H  +  hcost  -  *5  ( — <t> )  sin  ,  0  <  <J>  <  2tt  (C-36) 

n=2  n  -1 

2 

If  the  series  with  n  -1  in  the  denominator  are  used,  the 
n  =  1  terms  are  neglected. 
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F.  Asymptotic  Form  for  Ed  in  Circular  Waveguide 

The  asymptotic  form  for  E  can  be  obtained  from  that 
for  by  multiplying  the  latter  by  For  very  large  n, 

on  approaches  -  j  (n-*s)  Ti/b ,  and  the  subsequent  multiplication 
of  Eq.(C-34)  by  this  factor  yields  the  expression 


vn(P) 


(n~*5)  * 


/p 


-  cos{[_2  cos  (ni  )  - 
4 


-  sin  3  cos (n4 ,)  - 
4  J 


14  ■* 

-  cos  1  sin(n4.)  -  sinTl  cosin':.) 

—  1  —  1 

q 


*  1C 

sin\2  sin(n4-)}  -  -  {cos  3  sin(n4  ) 

4  2  4 


4  4  it 

cos_4  coslr.:,)  -  sm_4  sin(n4.);j  .  ( C— 3 7 ) 

A  **  A  ** 


The  convergence  improving  procedure  of  Appendix  B  may  be 
applied  through  the  use  of  the  series  in  Eqs.(C-20)  and 


APPENDIX  D 


PROPERTIES  OF  SOLUTIONS  TO  BESSEL'S  EQUATION 

1.  Introduction 

In  the  present  analysis,  frequent  use  of  various  Bessel 
functions  (solutions  to  Bessel's  equation)  is  made  and,  thus, 
for  convenience,  pertinent  properties  are  developed  and  sum¬ 
marized  here. 


2.  Orthogonality 

Bessel ' s  equation  is 

fn(Xx)  =  0 


I  dr  SL  - 

x  dx*  dx 


VAx)} 


,2  n 

A  -  2 

X 


where  f  is  any  solution.  Consider  two  solutions  B  (A  x)  and 
n  n  q 


Z  (y  x)  of 
m  p 


-  ~-{x  B  (A  x)  }  + 
x  dx 1  dx  n  q  1 


x2  -  £L_ 
q  2 
X 


B  ( Ax)  =  0 
n  q 


(D-la) 


and 


x  5c  Zm(Ypx)  *  + 


v2  .  5L_ 

YP  2 
x 


Z  ( Y_x)  =  0 
m  p 


(D-lb) 


Multiply  Eq.  (D-la)  by  x  zm(YpX)  and  Eq.  (D-lb)  by  x  B^A^x) 
to  obtain 
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r 


W>  fe 


x  3—  B  (Ax) 
dx  n  q 


+  (Xq  *  ^2)x  W*  VXqx)  =  0 


(D-2a) 


ana 


W>  fe 


,  X  -r—  Z  (y  x) 
[  dx  m  ' p 


+  (Y?  -  **>x  Z  ( Y  x)  B  ( A  x)  =  0 


P  2  m  p  n  q 


(D-2b) 


Subtract  Eq. (D-2b)  from  Eq. (D-2a)  and  integrate  over  (a,b) : 

rb. 


f  [Zm  IT'*  37  V  -  Bn  37(l‘  37  V]  dx 

-  fbf(V^>  +  iy(n2-m2)]x  2  B  dx 
J  L  p  q  2  J  m  n 


x=a 


(D-3) 


x=a 


The  left  side  of  Eq. (D-3) ,  after  one  integration  by  parts, 
reduces  to 
rb , 


f  [zm  37(’1  37  Bn>  -  Bn  37<x  37  V]  dx 
=a 

■  [x  Z  B  -  xB  Z  ] 

I  m  dx  n  n  dx  ml 

u  J  X=cl 


[*<37  Bn)l37  z*>  -  x(37  zm>(37  B„» 


x=a 


xZ  §-  B  -  x  B  |-Z 
m  dx  n  n  dx  m 


dx 


b 

x=a 
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so  that  one  may  write 


2  2  fk  2  2  fk  1 

(Y  -X  )  x  Z  B  dx  +  (n  -m  )  -  Z  B  dx 

'p  q  J  m  n  x  m  n 


=  [x  Z  4—  B  -  B  Z 

m  dx  n  n  dx 


a-z  V 


m  J 
J  x=a 


Eq. (D-4)  serves  as  the  basis  for  determination  of  various 
useful  properties  of  solutions  to  Bessel's  equation.  These 
properties  are  developed  in  turn  below.  When  m=n,  Eq.(D-4) 
reduces  to 

<''p-'q>  f-  VV  VV1  d* 

x=a 


*  *  W°  h  W>  *  VV>  Ic  W> 


Observe  that: 


x  Z  ( yx )  B(X  x)  dx  =  0 
n  p  n  q 


2  2 

if  Y  i*  x  and  if  the  following  is  true: 
P  ^ 


X  Z  (Y  X)  B  ( X  x)  -  B  ( X  x)  Z  (y  x)  =  0 

in  p  dx  n  q  n  q  n  p 


Z  and  B  can  be  the  same  solution  to  Bessel's  equation  or 
n  n 

different  solutions;  e.g.,  J  ,  N  ,  or  combinations  of  these. 

n  n 

If  the  term  within  the  braces  of  Eq. (D-7)  is  zero  at  both  a 
and  b,  or  one  is  zero  at  a  and  the  other  at  b,  the  condition 
is  satisfied. 

Another  condition  which  leads  to  similar  results  is 


x  T-  Z  +  Z  =  0 
dx  n  n 


(D-8a) 


x  ~  Bn  +  Bn  -  0 
dx  n  n 


(D-8b) 


at  x=a  and  at  x=b  as  can  be  seen  from  the  following  simplifi¬ 
cation  of  the  left-hand  side  of  Eq.(D-7)  subject  to  Eq. (D-8) : 


LlW)[-VXqb)]  ’  VXqb)[Zn(V)] 


LZn(Ypa)[-Bn(Aqa)]  -  Bn(Aqa)[Zn(Ypa)] 


=  0  ( D—  9 ) 


It  is  worth  noting  that  Eq.(D-9)  follows  from  Eq. (D-8)  also 

in  the  special  case  that  Z  =  B  (the  eigenvalues  are  main- 

n  n 

tained  distinct) . 

If  Eq. (D-la)  is  replaced  by  the  zero  eigenvalue  equation 
below: 


1  d_ 
x  dx 


d_ 

dx 


Rn(x) 


4  v*> 


0 


and  the  preceding  procedure  is  repeated,  one  arrives  at 
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X  Z  (Y„x)  R  (x)  dx 
n  p  n 


x=a 


iW1  E  Rn<x'  '  Rn(x)  E  W°  J 


ib 


(D-10) 


x=a 


which  shows  that  the  solution  Rn(x)  for  the  zero  eigenvalue 
exhibits  the  same  property  as  those  associated  with  non-zero 
eigenvalues . 


3 .  Normalization 


2  d 


It  Bessel's  equation  is  multiplied  by  2x  ^  fn(Ax), 


one  obtains 


2x  E  fn  E<x  E  fn>  +  2<x2x2-n2»E  fn  £n  ‘  0 


or 


d_ 

dx 


tx  E  fJ2 


-  "2  e([£i>I2)  +  x2(2x2  E  £n  £n)  ■  0  <D-ni 


By  direct  differentiation,  one  obtains 

J— [x2f 2]  =  2x2f  §-  f  +  2x  f2 
dx  n J  n  dx  n  n 


or 


2x2f  —  f  =  —  [x2f2l  -  2x  f2 
2x  fn  dx  n  dxlx  rnJ  2X  rn 


which  enables  one  to  reduce  Eq. (D-ll)  to 


2X2x  f2  =  — f  fx  —  f  l2 
2K  x  fn  dx[lx  dx  nJ 


+  -  "2  h[il]  (D'12) 
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Integration  of  both  sides  of  Eq. (d-12)  yields  the  following 
result: 


|bx  f2(Xx) 
x=a 


fn(Xx) 


2  2  2  2  2  2  1 
+  X  x^f„(Xx)  -  n^f’C-x) 
n  n  i 


(D-13) 


Another  form  which  is  useful  in  applications  follows  from  use 
of  the  identity  (valid  for  any  solution  to  Bessel's  equation) 

x  fn(Ax)  =  n  fn(Xx)  -  Xx  fn+1(^x) 

in  Eq. (D-13)  . 

f*5  2 

x  r  ( Xx)  dx  = 
n 

x=a 

for  n^O.  When  n=0,  fQ(Xx)  =»  -X  f-^fXx)  and  Eq.  (D-13)  is 
reduced  to 


■ff2 

n 


(Xx) 


-  fn+1ttx) 


fn-l(Ax) 


ib 


x=a 


(D-14a ) 


f  Q  (  X  x ) 


dx  =  j 


x2f2 (Xx) 


+  x2f2(Xx) 


-ib 


x=a 


x=a 


(D-14b) 


4  .  Orthogonality  of 

The  function  4>  ,  commonly  found  in  problems  dealing 
P 

with  coaxial  geometries,  is  defined  by 


♦  p(x) 


V  V> 


vy>  - 


W> 


v  V1 


(D-15) 
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where  a  and  b  are  the  inner  and  outer  radii,  respectively,  of 
the  coaxial  structure  and  the  eigenvalues  Xp  are  the  roots  of 

vy  w1  ■  vy  vy  • 0  1D‘16) 


Since  p(x)  is  a  solution  to  Bessel's  equation,  one  can 
write  from  Eq. (D-6) 


t'2  \2\ 

(a  -  A  ) 

p  q 


b 

x  $p(x)  dx 


$p(x) 


k  $q(X) 


vx) 


d_ 

dx 


4p(x) 


b 

x=a 


(D-17) 


Notice  that  the  right-hand  side  of  Eq.  (D-17)  is  zero  at  the 
upper  limit  because  $  (b)  and  $  (b)  are  zero  due  to  the  con- 

p  q 

ditions  that  X  and  X  are  the  appropriate  associated  roots. 

p  q 

At  the  lower  limit,  each  $  and  0  is  identically  zero  re- 

tr  '•4 

gardless  of  the  values  of  the  roots.  Hence,  clearly  and 
are  orthogonal. 


5.  Norm  of  $> 


From  Eq. (D-14b)  we  see  that  for  (order)  n=0: 

rb 


x  $  (x)  dx  = 
P 


x=a 


x*  N0U  a)  Jx(Xpx)  Jq  ( Xpa)  N^x) 


+  x2*p(x) 


b 

J  x=a 


191 


a 


which,  because  of  the  boundary  conditions  on  $p(x)  at  x  • 
and  x  ■  b,  reduces  to 


jbx  $2(x)  dx  -  l[x2[N0(Xpa)J1(Xpx)  -  jQfXpaJN^XpX)]2 


x»a 

or 

fb 

x*a 


i  b 
x»a 


x  *p(x)  dx  *  |[b2[NQ  (XpaJ^tXpb)  -  Jq  (Xpa)NL  (Xpb)  | 


-  a 


lN0(Xpa)J1(Xpa)  -  J0(Xpa)N 


i'V>)a] 


2 

Since  the  coefficient  of  a  above  is  in  the  form  of  a 
Wronskian  squared,  one  has 


rb 


x=»a 


•pW 


d!t  *  ir(w>Jlupb)  ' 


Another  form  is 


"p  -  [bx*2<x,dx  -  If  y<’)f  •  ID-19) 

x=a  p 

6.  Orthogonality  of  $p(x) 

The  eigenfunction  ^  $p(x) 

lx  VX)  -  ^ptNo<V)JlUpX)  "  Jo(V)Nl(V)]  (D“20) 
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where  the  Xp  are  the  roots  of  Eq. (D-17) .  Notice  that  the 
boundary  conditions  on  $p(x)  are  given  by 


and 


S  l«*pj 

1  *  0 

at  x  ■  a,b 

Hx(* 

d$P|  + 
dx 

fx2  -  y*  » 

IP  xj  P 

(D-21J 


(D-22) 


from  which  it  may  be  concluded  that 


dx^dx  $p(x) ^ 


+  h  *p(x)  *  0  at  x  *  a'b  •  (D— 23) 


Hence,  the  procedure  leading  to  Eq. (D-9)  ensures  the 
orthogonality  of  4p (x) . 

If  one  considers  the  zero  eigenvalue  equation  and  the 
corresponding  eigenfunction,  the  orthogonality  of  -  and 

X  X 

lx  $p^x^'  subject  to  the  weight  function  x,  can  be  demon¬ 
strated  by  writing 


fb‘(D( 


dx  *  «p{b)  -  <tp  (a) 


x«a 


I®  VX)J 

which  reduces  to,  in  view  of  the  boundary  conditions, 

b 


(D-24 ) 


fb  d 

x(3x  VX) 


2 

2a*  *  f-fe4p<x>j 


x*a 


x*a 


2fx2  2l b 

kp[T-fNo(Xpa)Ji(XPx)  '  Jo(XPa)Ni(V)]  j  *  (D~25) 
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7.  Norm  of  - 


The  eigenfunction  corresponding  to  the  zero  eigenvalue 
has  the  following  norm,  with  respect  to  the  weight  function 
x. 


(D-26) 


x*a 


8 .  Summary 

For  convenience,  the  orthogonality  properties  and  norms 
developed  in  this  Appendix  are  summarized  below. 
Orthogonality 


|  x  $p(x)  «g(x)  dx  -  0  ,  p  +  q 

x«a 


(D-27) 


rb  i 

X 

■a 

[dx 

*p(x)l 

life  V*’) 

dx 

*  o  ,  p  q 

(D-28) 

rb  , 

X 

d 

dx 

Vx\ 

(1)  -  ■ 

0 

'  xl  *  o 

(D-29) 

x=*a 


Norms 


n: 


rb  2  h2 

x  $„(x)  dx  ■  — s 

'  p  2X‘ 


Vx) 


x»a 


,2.  2 


£ 

2 


»o(XpllJi(V”  '  W’W’)2  ‘  77  (D"30) 


1  94 


I 


APPENDIX  E 


DETAILS  OF  COAXIAL  CAVITY  CONSTRUCTION 

1 .  Introduction 

Two  sets  of  coaxial  cavity  structures  were  fabricated  in 
order  to  perform  the  experimental  investigations  described  in 
Chapters  VTI  and  VTII.  The  two  sets  of  cavities  are  function¬ 
ally  quite  similar  but  differ  greatly  in  construction  detail 
and  materials.  The  primary  reason  for  this  difference  is  that 
one  set  is  almost  twice  the  diameter  of  the  other  and  the 
smaller  set  was  constructed  from  the  largest,  cylindrical, 
brass  pipe  available.  Each  set  consists  of  two  cavity  (wave¬ 
guide)  sections  having  identical  cross  sections,  but  differ¬ 
ent  lengths,  plus  interchangeable  end  plates  which  give  the 
cavities  one-  or  two-port  characteristics.  In  addition,  each 
set  includes  two  very  thin  partitions  with  centered,  circular 
apertures  (of  different  sizes)  which  are  used  to  divide  the 
cavities  into  multiple  sections.  All  cavity  components  are 
made  of  brass  in  order  to  achieve  boundary  conditions  of  high 
electrical  conductivity.  Following  are  short  physical  de¬ 
scriptions  of  each  set  of  cavities  and  appropriate  construc¬ 
tion  details. 
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2.  Smaller  Coaxial  Cavity 


The  geometry  and  dimensions  of  the  smaller  of  the  two 
sets  of  coaxial  cavities  are  given  in  Figure  E-l  and  Table 
E-l.  The  cavity  sections  were  cut  from  a  section  of  approxi¬ 
mately  eight  inch  inner  diameter,  red  brass  pipe  and  were 
bored  on  a  lathe  to  their  final  inner  diameters  after  the  1/4 
inch  thick  flanges  were  soldered  on  each  end.  As  indicated 
in  Table  E-l,  the  shorter  cavity  section  required  a  slightly 
larger  bore  to  achieve  a  circular  cross  section  which  left  it 
0.04  cm  larger  in  inner  diameter  than  the  other  section.  In 
order  to  position  the  inner  conductor  as  near  to  perfect  cen¬ 
ter  as  possible,  the  end  plates  and  cavity  sections  were 
aligned  by  inserting  steel  pins  through  precisely  positioned 
holes  drilled  in  flanges  and  the  end  plates.  After  in¬ 
sertion  of  the  pins,  nuts  and  bolts  were  used  to  securely 
fasten  the  adjacent  parts  together.  Use  of  the  available 
pieces  allowed  four  different  cavity  configurations  to  be 
realized: 


1. 

One-section, 

one-port 

cavity 

2. 

One-section, 

two-port 

cavity 

3. 

Two-section, 

one-port 

cavity 

4. 

Two-section, 

two-port 

cavity 

It  should  be  noted  that  configurations  3  and  4  have  the  op¬ 
tions  of  no  partition  between  the  sections  or  a  partition  with 
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TABLE  £-1 .  DIMENSIONS 

OF  COMPONENTS 

OF  SMALLER  SET  I 

OF  CAVITIES 

Cavity 

l 

5 

V 

5 

a 

T 

5 

Component 

[cm] 

[cm] 

[cm] 

[cm] 

[cm] 

[cm] 

[cm] 

Long  Cavity 
Section 

IS. 095 

30.48 

21.232 

0.635 

Short  Cavity 
Section 

8.637 

30.48 

21.273 

0.635 

Front  End 

Plate 

29.34 

11.684 

0.313 

Rear  End 

Plate 

29.84 

13.472 

0.318 

Partition  with 
Large  Aperture 

29.80 

8.006 

0.051 

Partition  with 
Small  Aperture 

29.80 

2.540 

0.051 

One-Port  Rod 

26.799 

Two-Port  Rod 

40.300 

Rod  Used  with 
Partition 

48.939 

TABLE  E-2.  COMPARISONS  OF  ADMITTANCES  OF  TRANSMISSION  LINES 
WITH  DIFFERENT  RODS  FOR  COAXIAL  CENTER  CONDUCTOR 


Transmission 

OO  of  Inner 
Conductor 

ID  of  Outer 
Conductor 

Admittance 

Line  Section 

[cm] 

[cm] 

[mhos] 

GenRad  50 -Ohm  Feed  Line 

GenRad  Center  Rod 

0.6205 

1.428 

0.020000 

Nominal  5*- inch  Rod 

0.6350 

1.428 

0.020565 

%-inch  Rod  -  0.005" 

0.6223 

1.428 

0.020065 

!*-inch  Rod  +  0.005" 

0.6477 

1.428 

0.021080 

Long  Cavity  Section 

GenRad  Canter  Rod 

0.6205 

21.232 

0.004718 

Nominal  %-inch  Rod 

0.6350 

21.232 

0.004749 

%-inch  Rod  -  0.005" 

0.6223 

21.232 

0.004722 

%-inch  Rod  +  0.005" 

0.6477 

21.232 

0.004770 

Short  Cavity  Section 

GenRad  Center  Rod 

0.6205 

21.273 

0.004715 

Nominal  %-inch  Rod 

0.6350 

21.273 

0.004746 

%-inch  Rod  -  0.005" 

0.6223 

21.273 

0.004719 

%-inch  Rod  +  0.005" 

0.6477 

21.273 

0.004773 
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one  of  two  aperture  sizes.  Thus,  there  are  eight  possible 
cavity  variations  available  from  the  four  basic  configurations. 

The  coaxial  outer  conductors  which  were  soldered  to  the 
end  plates  to  form  feed  lines  to  the  cavity  are  GenRad  Type 
874,  50-ohm,  air  line  sections.  The  companion,  50-ohm  inner 
conductors  were  used  for  cavity  configurations  1  and  2.  These 
rods  have  a  diameter  of  0.6205  cm  (0.2443  inches).  The  inner 
conductors  of  the  two-section  cavities  were  constructed  of 
circular,  brass  rod  with  a  diameter  of  0.635  ±0.013  cm  (0.250 
±0.005  inches).  The  possible  variations  in  the  characteris¬ 
tic  impedances  of  the  feed  lines  and  the  cavity  sections  is 
presented  in  Table  E-2.  This  change  in  inner  conductor  diam¬ 
eter  is  included  in  the  appropriate  calculations.  The  only 
impact  would  be  in  matching  the  50-ohm  load  for  the  two-port 
cavity  configuration.  Comparisons  of  measured  and  calculated 
results  for  these  cases  show  no  effects. 

In  order  to  allow  the  insertion  and  movement  of  probes 
designed  to  sample  the  field  inside  the  cavity,  radial  slots, 
0.3175  cm  (1/8  inch)  wide,  were  cut  in  each  of  the  end  plates. 
These  slots  were  designed  to  be  parallel  to  the  direction  of 
current  flow  on  the  end  plates  and,  thus,  provide  minimum  dis¬ 
turbance  to  the  field  within  the  cavity.  The  orientation  of 
the  slots  also  served  to  inhibit  the  excitation  of  any  higher- 
order  modes  which  might  be  associated  with  ^-directed  cur¬ 
rents. 
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The  use  of  the  brass  pipe,  the  1/4  inch  brass  flanges, 
and  the  1/8  inch  brass  end  plates  made  these  cavity  configu¬ 
rations  very  durable  and  heavy.  It  was  possible  to  tighten 
the  flange  bolts  very  firmly  in  order  to  enhance  electrical 
conductivity  at  the  interfaces.  The  excellent  agreement  be¬ 
tween  measured  and  calculated  field  and  admittance  values  for 
the  smaller  cavities  demonstrates  that  the  cavities  contain 
no  electrical  abnormalities.  Photographs  of  cavity  compo¬ 
nents  are  presented  in  Figure  E-2. 

3 .  Larger  Coaxial  Cavity 

The  geometry  and  dimensions  of  the  larger  set  of  coaxial 
cavities  are  given  in  Figure  E-3  and  Table  E-3 .  The  cavity 
sections  were  formed  from  0.005  cm  (0.020  inch)  thick  brass 
sheets  which  were  rolled  into  circular  cylinders  and  soldered 
to  1/8  inch  thick  brass  flanges  which  serve  to  hold  the  cylin¬ 
ders  in  shape.  The  axial  seams  where  the  sheet  brass  ends 
join  were  filled  with  solder  and  sanded  smooth.  Because  both 
cavity  sections  were  constructed  simultaneously,  the  inner 
diameter  (average  of  four  measurements)  of  one  section  agrees 
to  within  0.1  cm  of  the  other.  The  cavity  sections  and  end 
plates  are  joined  by  the  use  of  alignment  pins  and  bolts  pas¬ 
sing  through  holes  in  the  flanges  as  was  done  for  the  smaller 
set  of  cavities.  The  inner  conductor  was  constructed  from  a 
3.195  cm  (1.258  inch)  diameter  brass  pipe.  The  outer  conduc- 
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tor  of  the  feed-line  section  was  constructed  from  a  10.424  cm 
(4.104  inch)  inner  diameter,  brass  pipe. 

The  feed  point  on  the  feed  line  was  a  GenRad  50-ohm  con¬ 
nector  mounted  on  the  sedi  of  the  outer  conductor  with  the 
center  pin  of  the  connector  extending  down  to  the  inner  con¬ 
ductor  of  the  feed  line.  The  end  of  the  feed  line  not  at¬ 
tached  to  the  cavity  was  blocked  off  with  a  moveable  aluminum 
plug  which  served  as  a  shorting  plate.  Because  of  the  diffi¬ 
culty  of  fabricating  a  matched  load  for  the  non-standard  coax¬ 
ial  line,  only  one-port  versions  of  the  large  cavities  were 
constructed.  This  resulted  in  two  possible  cavity  configura¬ 
tions  : 

1.  One-section,  one-port  cavity 

2.  Two-section,  one-port  cavity 

Two  partitions  with  different  aperture  sizes  were  constructed 
for  U3e  with  the  two-section  cavity.  The  relevant  dimensions 
of  the  larger  cavity  set  are  given  in  Table  E-3  and  admittan¬ 
ces  of  different  waveguide  segments  are  given  in  Table  E-4. 

Slots  for  field  probe  insertion  were  cut  into  the  cavity 
end  plates  as  was  done  for  the  smaller  cavity  set.  In  addi¬ 
tion,  an  axial  slot  was  cut  in  the  inner  conductor  to  allow 
sampling  of  the  interior  field  along  the  axial  coordinate. 

This  slot  extended  well  into  the  feed-line  section  and  would 
serve  quite  well  for  a  slotted  line  for  admittance  measure¬ 
ments.  The  extremely  large  standing  wave  ratios  caused  by 
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Pin  Holes 


Figure  E-3.  Illustration  of  Components  of  Larger  Cavity  Set 


TABLE  E-3.  DIMENSIONS  OF  COMPONENTS  OF  LARGER  SET  OF  CAVITIES 


Cavity 

l 

5 

V 

6 

a 

T 

Component 

[cm] 

[cm] 

[cm] 

[cm] 

[cm] 

[cm] 

Long  Cavity 
Section 

27.318 

41.28 

35.631 

0.159 

Short  Cavity 
Section 

13.459 

41.28 

35.551 

0.159 

End  Plate  with 
Inner  Conductor 

41.28 

177.01 

0.159 

End  Plate  with 
Outer  Conductor 

41.28 

105.41 

0.159 

Partition  with 
Large  Aperture 

41.00 

14.813 

0.051 

Partition  with 
Small  Aperture 

41.00 

5.121 

0.051 

Truncated 
Coaxial  Center 
Conductor 

TABLE  E-4.  ADMITTANCE  OF  EACH  SECTION  OF  LARGER 
SET  OF  CAVITIES 


Transmission 
Lina  Section 

Feed  Line  Section 
Long  Cavity  Section 
Short  Cavity  Section 


OD  of  Inner 
Conductor 
[cm] 

3.195 

3.195 

3.195 


ID  of  Outer 
Conductor 
[cm] 

10.424 

35.631 

35.551 


Admittance 

[mhos] 

0.014094 

0.006911 

0.006917 
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the  one-port  cavity  configuration  made  such  measurements  im¬ 
practical  . 

There  was  some  difficulty  in  obtaining  extremely  tight 
interfaces  between  the  equally-spaced  bolts  which  hold  the 
flanges  together.  During  the  measurements,  C-clamps  were 
used  to  tighten  these  joints  although  tests  made  at  several 
frequencies  did  not  show  any  significant  effect  due  to  adding 
the  clamps.  The  agreement  between  measured  and  calculated 
results  on  the  end  plates  of  the  larger  cavities,  while 
generally  good  to  very  good,  did  not  attain  the  almost  com¬ 
plete  consistency  of  the  end-plate  measurements  for  the 
smaller  cavity  set.  The  greater  difficulty  in  obtaining  very 
tight  flange  junctions  is  one  likely  cause  for  such  disagree¬ 
ments.  The  very  good  agreement  between  measured  and  calcu¬ 
lated  data  at  the  surface  of  the  coaxial  inner  conductor  of 
the  large  cavity  contributes  strong  support  to  the  idea  that 
the  impedance  mismatch  of  the  one-port  cavities  makes  it 
extremely  difficult  to  couple  sufficient  energy  into  the 
cavity  to  perform  accurate  measurements. 


APPENDIX  F 


DESCRIPTIONS  OF  ELECTRIC  AND  MAGNETIC  FIELD  PROBES 

1 .  Introduction 

Two  types  of  probes  were  fabricated  in  order  to  sample 
the  interior  field  of  the  coaxial  cavities.  The  electric 
field  or  "charge"  probe  consists  of  a  short,  thin  electrical 
conductor  (in  this  case,  the  inner  conductor  of  a  short  piece 
of  0.035  inch  semi-rigid  coaxial  cable)  which  is  oriented  nor¬ 
mal  to  the  perfectly-conducting  wall  of  the  cavity.  The  nor¬ 
mal  component  of  the  electric  field  induces  a  small  current 
on  the  conductor  which,  in  turn,  creates  a  potential  differ¬ 
ence  between  the  conductor  and  the  cavity  wall  which  is  de¬ 
tectable  by  a  network  analyzer  operating  in  a  vector  voltmeter 
mode.  The  magnetic  field  or  "current"  probe  is  composed  of  a 
semi-circular  loop  of  0.035  inch  coax  with  the  outer  conductor 
grounded  to  the  probe  transport  at  both  ends.  A  short  section 
of  the  outer  conductor  is  cut  away  at  the  center  of  the  loop 
to  expose  a  short  piece  of  the  inner  conductor.  The  component 
of  magnetic  field  tangential  to  the  cavity  wall  induces  a  cur¬ 
rent  on  the  loop  and,  in  turn,  creates  a  potential  difference 
between  the  outer  and  inner  conductors  of  the  coaxial  loop. 
This  potential  difference  is  again  detected  by  a  network  ana- 


lyzer.  The  basic  probe  designs  have  been  in  use  by  Professor 
R.  W.  P.  King  and  others  for  many  years.  The  rest  of  this  ap¬ 
pendix  describes  the  particular  probe  transport  designs  used 
for  the  end  plates  and  the  center  conductor  of  the  cavities. 

2.  End-Plate  Probes 

Front  and  side  views  of  the  electric  and  magnetic  field 
probes  used  to  sample  the  cavity  field  at  the  inner  surface  of 
the  cavity  end  plates  are  presented  in  Figure  F-l.  The  probes 
are  mounted  on  brass  transports  which  have  a  T-shaped  cross 
section  for  insertion  into  and  sliding  back  and  forth  in  a  1/8 
inch  (0.3175  cm)  wide  slot.  The  bottoms  of  the  probe  trans¬ 
ports  are  designed  to  be  flush  with  the  inside  surface  of  the 
end  plates.  Good  electrical  contact  between  the  probe  trans¬ 
port  and  the  cavity  end  plate  is  very  important  in  order  to 
obtain  a  valid  reading  from  the  network  analyzer.  It  is  also 
important  that  the  probe  transport  slide  easily  in  the  slot  to 
facilitate  making  measurements.  These  two  contradicting  re¬ 
quirements  add  an  element  of  art  to  the  science  of  probe  manu¬ 
facture  and  use.  Adequate  pressure  on  the  top  of  the  probe 
transport  during  measurements  can  ensure  a  good  electrical 
contact  with  the  cavity  end  plate.  It  is  important  that  the 
bottom  of  the  probe  transport  be  flush  with  the  inner  surface 
of  the  end  plate  in  order  that  the  probe  see  no  sharp  corners 
or  other  discontinuities  in  its  neighborhood. 
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All  dimensions 
are  in  centimeters 


Since  the  end  plates  on  the  small  cavities  are  1/8  inch 
thick  and  those  of  the  large  cavities  sure  1/16  inch  thick,  a 
small,  flat,  1/16  inch  thick  clip  is  used  to  adjust  the  depth 
of  the  probe  transport  to  the  thickness  of  the  thinner  end 
plates.  A  photograph  of  the  two  probes  (with  a  quarter  for 
size  comparison)  is  shown  in  Figure  F-2. 

3.  Center  Conductor  Probes 

Drawings  of  the  electric  and  magnetic  field  probes  used 
to  sample  the  cavity  field  at  the  outer  surface  of  the  coaxial 
inner  conductor  are  presented  in  Figure  F-3 .  The  probes  are 
mounted  on  1/8  inch-thick,  brass  guides  which  sure  in  turn 
mounted  on  thin,  brass  plates.  The  plates  have  been  contoured 
to  fit  the  outside  of  a  brass,  cylindrical  probe  transport  and 
make  a  snug  fit  with  the  inner  surface  of  the  brass  pipe  which 
forms  the  coaxial  center  conductor.  The  thin  plates  are  se¬ 
cured  to  the  sides  of  the  cylindrical  probe  transports  by  min¬ 
iature  screws.  The  coaxial  lead  from  each  probe  passes  into  a 
hole  in  the  side  of  the  cylinder  and  emerges  from  the  hole  at 
one  of  the  flat  ends.  The  probe  transport  slides  easily  in¬ 
side  the  brass  pipe  because  of  thick.  Teflon  disks  which  are 
fastened  to  each  end  of  the  brass  cylinder.  The  disks  are 
slightly  larger  than  the  cylinders  and  make  the  sliding  con¬ 
tact  with  the  inner  walls  of  the  pipe.  The  thin  brass  plates 
make  the  primary  electrical  contact  with  the  brass  pipe.  The 


tops  of  the  brass  guides  are  designed  to  be  flush  with  the 
outer  surface  of  the  brass  pipe. 


The  probe  transports  are  positioned  by  moving  a  long, 
thin,  brass  tube  whose  end  has  been  threaded  to  receive  a 
bolt  which  passes  through  the  center  of  the  brass  cylinder. 
The  flexible  probe  lead  enters  a  hole  in  the  side  of  the  thin 
tube  (just  behind  the  threaded  section) ,  passes  through  the 
inside  of  the  tube,  and  exits  at  the  far  end  of  the  tube.  A 
marker  was  placed  near  the  fair  end  of  the  tube  for  use  in  re¬ 
cording  the  position  of  the  probe  inside  the  cavity.  Two 
photographs  of  the  magnetic  and  electric  field  probes  are 
presented  in  Figure  F-2. 
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